gr-qc/0311052 



Dimensional regular izat ion of the third post-Newtonian dynamics 
of point particles in harmonic coordinates 

Luc Blanchet/'B Thibault Damour,^'Q and Gilles Esposito-Farese^'H 

^gReCO, FRE 2435-CNRS, Institut d'Astrophysique de Paris, 
ggbis iioulevard Arago, F-75014 Paris, France 
^Institut des Hautes Etudes Scientifiques, 
35 route de Chartres, F-91440 Bures-sur-Yvette, France 
(Dated: November 17, 2003) 

Abstract 

Dimensional regularization is used to derive the equations of motion of two point masses in 
harmonic coordinates. At the third post-Newtonian (3PN) approximation, it is found that the 
dimensionally regularized equations of motion contain a pole part [proportional to {d — 3)~^] which 
diverges as the space dimension d tends to 3. It is proven that the pole part can be renormalized 
away by introducing suitable shifts of the two world-lines representing the point masses, and that 
the same shifts renormalize away the pole part of the "bulk" metric tensor g^y[x^). The ensuing, 
finite renormalized equations of motion are then found to belong to the general parametric equations 
of motion derived by an extended Hadamard regularization method, and to uniquely determine the 
3PN ambiguity parameter A to be: A = —1987/3080. This value is fully consistent with the recent 
determination of the equivalent 3PN "static ambiguity" parameter, = 0, by a dimensional- 
regularization derivation of the Hamiltonian in Arnowitt-Deser-Misner coordinates. Our work 
provides a new, powerful check of the consistency of the dimensional regularization method within 
the context of the classical gravitational interaction of point particles. 
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I. INTRODUCTION 



A. Relation to previous works 

The problem of motion, one of the cardinal problems of Einstein's gravitation theory, 
has received continuous attention over the years. The early, classic works of Lorentz-Droste, 
Eddington- Clark, Einstein-Infeld-Hoffmann, Fock, Papapetrou and others led to a good 
understanding of the equations of motion of bodies at the first post-Newtonian (IPN) 
approximation^ (see, e.g., P for a general review of the problem of motion). In the 1970's, 
an important series of works by a Japanese group jl, IE 1^ led to a nearly complete control 
of the problem of motion at the second post-Newtonian (2PN) approximation. Then, in the 
early 80 's, motivated by the observation of secular orbital effects in the Hulse- Taylor binary 
pulsar PSR1913-I-16, several groups solved the two-body problem at the 2.5PN level (while 
completing on the way the derivation of the 2PN dynamics) 0, S 0, i, S H Q Q (fo^ 
more recent work on the 2.5PN dynamics see 3 3 3)- 



In the late 90 's, motivated by the aim of deriving high-accuracy templates for the data 
analysis of the upcoming international network of interferometric gravitational-wave detec- 
tors, two groups embarked on the derivation of the equations of motion at the third post- 
Newtonian (3PN) level. One group used the Arnowitt-Deser-Misner (ADM) Hamiltonian 
approach [ly, [13, [H, 113 worked in a corresponding ADM-type coordinate system. 
Another group used a direct post-Newtonian iteration of the equations of motion in har- 
monic coordinates 21 , 2^, 2^ 2^, IE ipl. The end results of these two approaches have been 



proved to be physically equivalent |20l.l25|. However, both approaches, even after exploiting 
all symmetries and pushing their Hadamard-regularization-based methods to the maximum 
of their possibilities, left undetermined one and only one dimensionless parameter: in the 
ADM approach and A in the harmonic-coordinates one. The unknown parameters in both 
approaches are related by 

X 3 1987 

A = Us , 1.1 

11 3080 ^ ^ 

as was deduced from the cornparison between the invariant energy functions for circular 
orbits in the two approaches j2l|, and from two independent proofs of the equivalence be- 
tween the two formalisms for general orbits [l^, l25| |. The appearance of one (and only 
one) unknown parameter in the equations of motion is quite striking; it is related with the 
choice of the regularization method used to cure the self-field divergencies of point particles. 
Both lines of works [H, 113, E, dS, M and |2l|, 22^ 2A, 25., M\ regularized the self-field 



divergencies by some version of the Hadamard regularization method. The second line of 
work defined an extended version of the Hadamard regularization (isl . 23 |. which permitted 



a self-consistent derivation of the 3PN equations of motion, but its use still allowed for the 
presence of arbitrary parameters in the final equations. On the other hand, the Hadamard 
regularization also yielded some arbitrary parameters in the gravitational radiation field of 
point-particle binaries at the 3PN order^he most important of which being the parameter 
9 entering the binary's energy fiux (j?! . |2^ . 

Let us notice that the regularization (when dealing with point particles) and the renor- 
malization (needed when dealing either with point particles or with extended bodies^) of 



^ As usual the nPN order refers to the terms of order 1 / c^" in the equations of motion. 

^ In the case of extended compact bodies, the gravitational self-energy (divergent when the radii of the 
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self-field effects has recurrently plagued the general relativistic problem of motion. Even at 
the IPN level, early works often contained incorrect treatments of self- field effects (see, e.g., 
Section 6.14 of for a review). At the 2PN level, the self- field divergencies are more severe 
than at the IPN level. For instance, they caused Ref. [3] to incorrectly evaluate the "static" 
{i.e., velocity-independent) part of the 2PN two-body Hamiltonian. The first correct and 
complete evaluation of the 2PN dynamics has been obtained by using the Riesz analytical 
continuation method [2^ . (See [1, for a detailed discussion of the evaluation of the static 
2PN two-body Hamiltonian.) In brief, the Riesz analytical continuation method consists 
of replacing the problematic delta-function stress-energy tensor of a set of point particles 

2/^(5a), 

T-^(x) = Y.^a^ j ds/jf-'^^^ [-g{ya)r" S^'\x' - y'^is^)) , (1.2) 
[where ds^^ = —g^uiUa) dVady^, g = detgfj,u] by an auxiliary, smoother source 

Tr(a:) = J2'^ac'fdsa^^^^ [-giyaT'^' Zf{x' - yl{s^)) . (1.3) 

a 

[Actually, in the implementation of 0, one works with T'^^(x) = \g{x)\ T^'^(x).] In Eq. ()1.3|1 
the four-dimensional delta function entering Eq. ()1.2j) has been replaced by the Lorentzian ^ 

four- dimensional Riesz kernel Z^\x — y), which depends on the complex number A. When 
the real part of A is large enough the source T'^{x) is an ordinary function of x^, which is 
smooth enough to lead to a well-defined iteration of the harmonically relaxed Einstein field 
equations, involving no divergent integrals linked to the behavior of the integrands when 
x^ — s> y^. One then analytically continues A down to 0, where the kernel Z^^{x — y) tends 
to 5^^\x — y). The important point is that it has been shown that all the integrals 
appearing in the 2PN equations of motion are meromorphic functions of A which admit a 
smooth continuation at A = (without poles). It was also shown there that the formal 
construction based on (jl.3p does generate, at the 2.5PN level, the metric and equations of 
motion of "compact" bodies {i.e., bodies with radii comparable to their Schwarzschild 
radii). 

The Riesz analytic continuation method just sketched works within a normal 4- 
dimensional space-time (as recalled by the superscript (4) in ()1.3|) ). However, it was 
mentioned in |3^ that the same final result (at the 2.5PN level) is obtained by replac- 
ing Z^\x — y) by Zq^ ^\x — y) = 6^'^~"^\x — y), i.e., by formally considering delta-function 
sources in a space-time of complex dimension 4 — A. In other words, at the 2.5PN level, the 
Riesz analytic continuation method is equivalent to the dimensional regularization method.^ 
However, it was also mentioned at the time P| that the generalization of Riesz analytic con- 
tinuation beyond the 2.5PN level did not look straightforward because of the appearance of 
poles, proportional to A~^, at the 3PN level (when using harmonic coordinates). 



bodies formally tend to zero) must be renormalized into the definition of the mass. 
^ Z^"^^-* is the Lorentzian version of the Euclidean kernel Slf' discussed in Appendix lEl 
^ Dimensional regularization was invented as a mean to preserve the gauge symmetry of perturbative quan- 
tum gauge theories ,31„ ,,32„ 33». i34] . Our basic problem here is to respect the gauge symmetry associated 
with the diffeomorphism invariance of the classical general relativistic description of interacting point 
masses. 
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Recently, Damour, Jaranowski and Schafer j35[ showed how to use dimensional regular- 
ization within the ADM canonical formalism. They found that the reduced Hamiltonian 
describing the dynamics of two point masses in space-time dimension D = d + 1 was finite 
(no pole part) as — 3. They also found that the unique 3PN Hamiltonian defined by the 
analytic continuation of d towards 3 had two properties: (i) the velocity-dependent terms 
had the unique structure compatible with global Poincare invariance,^ and (ii) the velocity- 
independent ("static") terms led to an unambiguous determination of the unknown ADM 
parameter a;^, namely 

^dim. reg. ADM _ Q _ (1.4) 

The fact that the dimensionally regularized 3PN ADM Hamiltonian ends up being globally 
Poincare invariant is a confirmation of the consistency of dimensional regularization, because 
this symmetry is not at all manifest within the ADM approach which uses a space-plus-time 
split from the start. By contrast, the global Poincare symmetry is manifest in harmonic 
coordinates, and indeed the 3PN harmonic-coordinates equations of motion derived in pli 
2^ were found to be manifestly Poincare invariant. 



B. Method and main results 



In the present paper, we shall show how to implement dimensional regularization (hence- 
forth often abbreviated as "dim. reg." or even "dr") in the computation of the equations of 
motion in harmonic coordinates, i.e., following the same iterative post-Newtonian formal- 
ism as in Refs. |l3l . 21, i^]. Similarly to the ADM calculation of Ref. jl^, our strategy 



will essentially consist of computing the difference between the (i-dimensional result and 
the 3-dimensional one j^, 23| corresponding to Hadamard regularization. This difference is 



computed in the form of a Laurent expansion in e = d — 3, where d denotes the spatial di- 
mension. The main reason for computing the e-expansion of the difference is that it depends 
only on the singular behavior of various metric coefficients in the vicinity of the point par- 
ticles, so that the functions involved in the delicate divergent integrals can all be computed 
in d dimensions in the form of local expansions in powers of ri or r2 (where = |x — y^l; 
Ya, a = 1,2, denoting the locations of the two point masses). Dimensional regularization as 
we use it here can then be seen as a powerful argument for completing the 3-dimensional 
Hadamard-regularization results of 0, i^l and fixing the value of the unknown parameter. 



We leave to future work the task of an exact calculation of the (i-dimensional equations of 
motion, instead of the calculation of the first few terms in a Laurent expansion in e around 
d = 3, as done here. The first step towards such a calculation is taken in Appendix 
where we give the explicit expression of the basic quadratically non-linear Green function 
^■(x, yi, y2) in d dimensions. 

The detailed way of computing the difference between dim. reg. and Hadamard's reg. will 
turn out to be significantly more intricate than in the ADM case. This added complexity has 
several sources. A first source of complexity is that the harmonic-gauge d-dimensional cal- 
culation will be seen to contain (as anticipated long ago 0) poles, proportional to {d — 3)'^, 

^ Thus the "kinetic ambiguity" parameter uJk, originaUy introduced in the ADM approach , takes 



the unique value Wfc = This value was obtained in |21| using the result for the binary energy function 
in the case of circular orbits, as calculated in the harmonic-coordinates formalism, and also directly from 
the requirement of Poincare invariance in the ADM formalism [loj . 
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by contrast to the ADM calculation which is finite as c? 3. A second source of complex- 
ity is that the end results 2^1 for the 3-dimensional 3PN equations of motion have been 



derived using systematically an extended version of the Hadamard regularization method, 
incorporating both a generalized theor y o f singular pseudo-functions and their associated 
(generalized) distributional derivatives [2^, and an improved definition of the finite part as 
X — >• yi, say [F]i, of a singular function F(x, yi,y2), designed so as to respect the global 
Poincare symmetry of the problem 2^1 . We shall then find it technically convenient to sub- 



tract the various contributions to the end results of [2^ which arose because of the specific 
use of the extended Hadamard regularization methods of ji^ 2^1 before considering the dif- 



ference with the (i-dimensional result. A third source of added complexity (with respect to 
the ADM case^) comes from the presence in the harmonic-gauge integrals we shall evaluate 
of "hidden-distributional" terms in the integrands. By hidden distributional terms we mean 
terms proportional to the second spatial derivatives of the Poisson kernel di^^ oc or 
to the fourth spatial derivatives of the iterated Poisson kernel A^^ 6^^^ ocrf"'^. Such terms, 
dij rl~^ or dijki r^~'^, considered as Schwartz distributional derivatives [36| , contain pieces 

proportional to the delta function S g^^ , which need to be treated with care. The gen eralized 
distributional derivative defined in p^, and used to compute the end results of i^l, led to 



an improved way, compared to the normal Schwartz distributional derivative, of evaluating 
contributions coming from the product of a singular function and a derivative of the type 
dijr~^ or dijkiTa, and more generally of any derivatives of singular functions in a certain 
class. We shall find it convenient to subtract these additional non-Schwartzian contributions 
to the 3PN equations of motion before applying dimensional regularization. However, we 
shall note at the end that dim. reg. automatically incorporates all of these non-Schwartzian 
contributions. 

A fourth, but minor, source of complexity concerns the dependence of the end results of 
2^ 1 for the 3PN acceleration of the first particle (label a = 1), say a^^, on two arbitrary 



length scales r[ and rg, and on the "ambiguity" parameter A. Explicitly, we define 

af^ [A; r'l, r^] = R.H.S. of Eq. (7.16) in Ref. . (1.5) 



Here the acceleration is considered as a function of the two masses mi and m2, the relative 
distance yi — y2 = ri2ni2 (where ni2 is the unit vector directed from particle 2 to particle 
1), the two coordinate velocities vi and V2, and also, as emphasized in ()1.5|) . the parameter 
A as well as two regularization length scales r[ and r!^. The latter length scales enter the 
equations of motion at the 3PN level through the logarithms ln(ri2/r^) and In(ri2/r2). They 
come from the regularization as the field point x' tends to yi or y2 of Poisson-type integrals 
(see Section [111 Bl below) . The length scales r[, r!^ are "pure gauge" in the sense that they 
can be removed by the effect induced on the world-lines of a coordinate transformation of 
the bulk metric |22]. On the other hand, the dimensionless parameter A entering the final 



^ The specific form of the 3PN ADM Hamiltonian H derived in [1^ and used (in its rf-dimensional general- 
ization) in was written, on purpose, in a way which does not involve any hidden distributional terms 
(the only delta- function contributions it contains being explicit contact terms F{x)Sa )• This allowed 
one to estimate the difference between the d-dimensional Hamiltonian ij'''^ and the Hadamard-regularized 
3-dimensional one Hr[i?(^^] without having to worry about distributional derivatives. However, as a check 
on the consistency of dim. reg., the authors of |35j did perform another calculation of H based on a 
starting form of the Hamiltonian which involved hidden distributional terms, with the same final result. 
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result p.5|) corresponds to genuine physical effects. It was introduced by requiring that the 
3PN equations of motion admit a conserved energy (and more generally be derivable from 
a Lagrangian). This extra requirement imposed two relations between the two length scales 
r[, and two other length scales Si, S2 entering originally into the formalism, namely the 
constants Si and S2 parametrizing the Hadamard partie finie of an integral as defined by Eq. 
(|3.4|) below. These relations were found to be of the form 

l„(!l).l| + A!!!l±i!^ andl«2, (1.6) 

VS2/ 308 1712 

where the so introduced single dimensionless parameter A has been proved to be a purely 
numerical coefficient (independent of the two masses). When estimating the difference be- 
tween dim. reg. and Hadamard reg. it will be convenient to insert Eq. into (|1.5p 
and to reexpress the acceleration of particle 1 in terms of the original regularization length 
scales entering the Hadamard regularization of ai, which were in fact r'l and S2- Thus we 
can consider alternatively 

af^[r;,S2] = af^[A;r;,r^(s2,A)] and 1 ^ 2 , (1.7) 

where the regularization constants are subject to the constraints p.6|) [we will then check 
that the A-dependence on the R.H.S. of ()1.7|) disappears when using Eq. ()1.6p to replace r2 
as a function of S2 and A]. 

Our strategy will consist of two steps. The first step consists of subtracting all the 
extra contributions to Eq. ()1.5p . or equivalently Eq. (jl.TJ L which were specific consequences 
of the extended Hadamard regularization defined in j23l 2^ 1 . As we shall detail below, 



there are seven such extra contributions ^'^ai, A = 1, ■ ■ ■ ,7. Essentially, subtracting these 
contributions boils down to estimating the value of ai that would be obtained by using a 
"pure" Hadamard regularization, together with Schwartz distributional derivatives. Such a 
"pure Hadamard-Schwartz" (pHS) acceleration was in fact essentially the result of the first 
stage of the calculation of ai, as reported in the (unpublished) thesis ji^]. It is given by 

ar[ri,.2] =af K,S2]-X:5V. (1.8) 

A=l 

The second step of our method consists of evaluating the Laurent expansion, in powers of 
e = d — 3, of the difference between the d-dimensional and the pure Hadamard-Schwartz (3- 
dimensional) computations of the acceleration ai. We shall see that this difference makes a 
contribution only when a term generates a pole ~ 1/e, in which case dim. reg. adds an extra 
contribution, made of the pole and the finite part associated with the pole [we consistently 
neglect all terms 0{e)]. One must then be especially wary of combinations of terms whose 
pole parts finally cancel ( "cancelled poles" ) but whose dimensionally regularized finite parts 
generally do not, and must be evaluated with care. We denote the above defined difference 

Pai = Dai [e, 4; r[, S2] = Vsl^ [e, 4; A; r[, r^] . (1.9) 

It depends both on the Hadamard regularization scales r[ and S2 (or equivalently on A and 
r[, Tg) and on the regularizing parameters of dimensional regularization, namely e and the 
characteristic length associated with dim. reg. and introduced in Eq. ()2.4|) below. We 
shall explain in detail below the techniques we have used to compute Va.i (see Section HVj). 
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Finally, our main result will be the explicit computation of the ^-expansion of the dim. reg. 
acceleration as 

af [e, io] = af '[ri, + Vai,[e, io; r[, S2] . (1.10) 
With this result in hands, we shall prove (in Section rVT|) two theorems. 

Theorem 1 The pole part oc 1/e of the dimensionally-regularized acceleration m.l(J\} . as 
well as of the metric field g^^{x) outside the particles, can be re-absorbed (i.e., renormalized 
away) into some shifts of the two "bare" world-lines: Ya ^ Ya + $,a, with, say, oc l/e 
("minimal subtraction"; MS), so that the result, expressed in terms of the "dressed" quanti- 
ties, is finite when e 0. 

The situation in harmonic coordinates is to be contrasted with the calculation in ADM-type 
coordinates within the Hamiltonian formalism [35j, where it was shown that all pole parts 
directly cancel out in the total 3PN Hamiltonian (no shifts of the world- lines were needed). 
The central result of the paper is then as follows. 

Theorem 2 The "renormalized" (finite) dimensionally-regularized acceleration is physi- 
cally equivalent to the extended- Hadamard-regularized acceleration ( end result of Ref. \23il ). 
in the sense that there exist some shift vectors ^i{B,io;r[) and £2(£^, ^0; ^^2)? such that 

[A, r[,r',] = lim [af [e, io] + %A;ri,r^) ai] (1.11) 

( where ai denotes the effect of the shifts on the acceleration ^ if and only if the heretofore 
unknown parameter X entering the harmonic- coordinates equations of motion takes the value 

1Q87 

^dim. reg. harmonic /j^ "1^2^ 

3080 ■ ^ ■ ^ 

The precise shifts ^a(e) needed in Theorem 2 involve not only a pole contribution oc l/e, 
which defines the "minimal" (MS) shifts considered in Theorem 1, but also a finite contri- 
bution when e ^ 0. Their explicit expressions read: 



11 Glmj n _ 2 In /^i^^^^ ^ ^27 



c6 



e \ in 1540 



aivi and 1^2, (1-13) 



where Gn is the usual Newton's constant [see Eq. ()2.4|) below], a^vi denotes the acceleration 
of the particle 1 (in d dimensions) at the Newtonian level, and q = Ane'-' depends on the 
Euler constant C = 0.577 ■ ■ ■ . [The detailed proofs of Theorems 1 and 2 will consist of 
our investigations expounded in the successive sections of the paper, and will be completed 
at the end of Sections IVIBI and I VI (.1 respectively, taking also into account the results of 
Section IVTDI ] 

Notice that an alternative way of presenting our central result is to say that, in fact, 
each choice of specific renormalization prescription (within dim. reg.), such as "minimal 



^ When working at the level of the equations of motion (not considering the metric outside the world-lines), 
the effect of shifts can be seen as being induced by a coordinate transformation of the bulk metric as in 
Ref. (we comment on this point in Section IVl Bl belowl . 
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subtraction" as assumed in Theorem 1 for conceptual simplicity,^ leads to renormalized 
equations of motion which depend only on the dim. reg. characteristic length scale £o 
through the logarithm \n{ri2/io), and that any of these renormalized equations of motion 
are physically equivalent to the final results of 2^1 . In particular, this means, as we shall see 
below, that each choice of renormalization prescription within dim. reg. determines the two 
regularization length scales r[, r'^ entering Eq. (jl.5|) . Of course, what is important is not the 
particular values these constants can take in a particular renormalization scheme [indeed 
r[ and are simply "gauge" constants which can anyway be removed by a coordinate 
transformation], but the fact that the different renormalization prescriptions yield equations 
of motion falling into the "parametric" class {i.e., parametrized by r[ and rg) of equations 
of motion obtained in 23| . 



An alternative way to phrase the result p.lH) - ()1.12j) . is to combine Eqs. fll.Sp and (jl.lOp 
in order to arrive at 



lim 

e^O L 



Pai[£,4;-ii;r;,r^] + aJ = J] 5 % . (1.14) 



A=l 



Under this form one sees that the sum of the additional terms 6^ai differs by a mere shift, 
when and only when A takes the value ()1.12|) . from the specific contribution Vai we shall 
evaluate in this paper, which comes directly from dimensional regularization. Therefore 
one can say that, when A = — IH^, the extended-Hadamard regularization jl^ [i^ is in 
fact (physically) equivalent to dimensional regularization. However the extended-Hadamard 
regularization is incomplete, both because it is unable to determine A, and also because it 
necessitates some "external" requirements such as the imposition of the link ()1.6|) in order 
to ensure the existence of a conserved energy — and in fact of the ten first integrals linked to 
the Poincare group. By contrast dim. reg. succeeds automatically (without extra inputs) in 
guaranteeing the existence of the ten conserved integrals of the Poincare group, as already 
found in Ref. [35^ . 

In view of the necessary link (jl.lj) provided by the equivalence between the ADM- 
Hamiltonian and the harmonic-coordinates equations of motion, our result ()1.12j) is in perfect 
agreement with the previous result (|1.4j) obtained in |35|i].^ Our result is also in agreement 
with the recent finding of Itoh and Futamase [ssl. l39| (see also jl3]); who derived the 3PN 
equations of motion in harmonic gauge using a "surface-integral" approach, aimed at describ- 
ing extended relativistic compact binary systems in the strong-field point particle limit. The 
surface-integral approach of Refs. ^] is interesting because, like the matching method 
used at 2.5PN order in it is based on the physical notion of extended compact bodies. 



However, for technical simplicity we shall prefer in Section [VI CI below to use a modified minimal subtrac- 
tion that we shall denote 

One may wonder why the value of A is a complicated rational fraction while tOs is so simple. This is because 
LOs was introduced precisely to measure the amount of ambiguities of certain integrals, and that the ADM 
Hamiltonian reported in IS] was put in a minimally ambiguous form, already in three dimensions, for 
which an a posteriori look at the "ambiguities" discussed in the Appendix A of [l8| already showed that 
LOg = 0. By contrast, A has been introduced as the only possible unknown constant in the link between 
the four arbitrary scales r'i,r2, si, S2 (which has a priori nothing to do with ambiguities of integrals), in 
a framework where the use of the extended Hadamard regularization makes in fact the calculation to be 
unambiguous. 
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In this respect, we recall that the matching method used in Q showed that the internal 
structure (Love numbers) of the constituent bodies would start influencing the equations of 
motion of (non-spinning) compact bodies only at the 5PN level. This effacement property 
strongly suggests that it is possible to model, in a physically preferred manner, two compact 
bodies as being two point-like particles, described by two masses and two world-lines, up to 
the 4.5PN level included. It remains, however, to prove that the dimensional regularization 
of delta-function sources does yield the physicall;^unique equations of motion of two com- 

proved it at the 2.5 PN level, and the 
m indicates that this is also true at the 
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pact bodies up to the 4.5 PN order. The work 
agreement of the present results with those of jSS 
3PN level. 

Besides the independent confirmation of the value of uts or A, let us also mention that 
our work provides a confirmation of the consistency of dim. reg., because our explicit 
calculations [which involved combinations of hundreds of Laurent expansions of the form 
a_i + aQ + 0{e)] are entirely different from the ones of js^: We use harmonic coordinates 
(instead of ADM- type ones), we work at the level of the equations of motion (instead of the 
Hamiltonian), we use a different form of Einstein's field equations and we solve them by a 
different iteration scheme. 

Finally, from a practical point of view our confirmation of the value of Ug orA allows one 
to use the full 3PN accuracy in the analytical computation of the dynamics of the last orbits 
of binary orbits [i^, 41 1. It remains, however, the task of computing, using dimensional 



regularization, the parameter 6 entering the 3.5PN gravitational energy flux j27„ ^] to be 
able to have full 3.5PN accuracy in the computation of the gravitational waveforms emitted 
by inspiralling compact binaries (see, e.g., |42| and references therein). 

The organization of the paper is as follows. In Section IH] we derive our basic 3PN solution 
of the field equations for general fluid sources in d spatial dimensions, using (i-dimensional 
generalizations of the elementary potentials introduced in Ref. ji^]. Section UTTl collects all 
the additional terms included in j22] which are due specifically to the extended Hadamard 
regularization, and derives the pure Hadamard-Schwartz (pHS) contribution to the equations 
of motion. The differences between the dimensional and pHS regularizations for all the 
potentials and their gradients are computed in Section IIVI Then the dim. reg. equations 
of motion are obtained in Section |3 where we comment also on their interpretation in 
terms of space-time diagrams. Section IVIl is devoted to the renormalization of the dim. reg. 
equations by means of suitable shifts of the particles' world-lines, and to the equivalence 
with the end results of ji^] when Eq. (jl.l2j) holds. At this stage, the proofs of Theorems 1 
and 2 stated above are finally completed. 

We end the paper with some conclusions fSection fVII|) and three appendices. Appendix 1X1 
provides further material on the (i-dimensional metric and geodesic equation. Appendix ^ 
gives a compendium of useful formulae for working in d dimensions, and Appendix gen- 
eralizes the well-known quadratic-order elementary kernel = ln(ri + r2 + ri2) to d 
dimensions. The latter calculation of the (i-dimensional kernel is not directly employed 
in the present paper, but represents a first step in obtaining the equations of motion in any 
dimension d (not necessarily of the form 3 + e). 
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II. FIELD EQUATIONS IN d + 1 SPACE-TIME DIMENSIONS 



This section is devoted to the field equations of general relativity in ci + 1 space-time 
dimensions, and to the geodesic equation describing the motion of point particles. We use 
the sign conventions of Ref. ji^, and in particular our metric signature is mostly +. Space- 
time indices are denoted by greek letters, and spatial indices by latin letters ■ ■ ■ run 
from 1 to d). A summation is understood for any pair of repeated indices. We work in the 
harmonic gauge, which is such that 

r^ = ^?°^r^^ = o. (2.1) 

As usual, g°'^ denotes the inverse metric and F^^ the Christoffel symbols. Using this gauge 
condition, one can easily prove that the Ricci tensor reads in any dimension 

~^9fJ,a,'y 9/35,1/ ~^ 9va,'^f (2.2) 

where a comma denotes partial derivation. Note that the spatial dimension d does not 
appear explicitly in this expression, whereas some (i-dependent coefficients do appear when 
expressing the Ricci tensor in terms of the so-called "gothic" metric g^"^ = y/—g g^'^ [see 
Eq. (|X9| in Appendix \M below] . 

In any dimension, the Einstein field equations read 

_\ ^,u^^^]^rJ.^,^ (2.3) 
2 V y 

where T^^^ denotes the matter stress-energy tensor, given by the functional derivative 
^yzr^rptiu ^ 2cSSm/Sg^iy of the matter action Sm with respect to the metric tensor. By 
definition, G denotes the constant involved in Eq. (|2.3|) . which shows that its dimension is 
such that 

G = G^<-3, (2.4) 

where Gjv is the usual Newton constant (in 3 spatial dimensions) and io is an arbitrary 
length scale. This scale will be involved in our dimensionally regularized results below, but 
we will finally show that the physical observables do not depend on it. As is well known, 
the combination of Eq. ()2.3|) with its trace allows us to rewrite it as 

8nG 1 



R,^ = - ^ 9^.n\ , (2.5) 



in which the spatial dimension d now appears explicitly. 

We wish to expand in powers of 1/c the field equations resulting from ()2.2j] and ()2.5j) . The 
basic idea is to introduce a sequence of "elementary potentials", V , Vi, Wij, . . .which allow 
one to parametrize conveniently the successive post-Minkowskian contributions to the metric 
gniy{x). For instance, at the first post-Minkowskian order it is convenient to parametrize the 
metric as 

goo = -l + 2V/c' + 0{G'), go^ = -^/c^ + 0{G^) , (2.6) 
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where the so-introduced elementary potentials V and Vi satisfy equations of the form 



nV = -AnGa , nVi = -AnGai 



(2.7) 



where n = df — {l/(?)d^ denotes the flat d'Alembertian and where, hy definition, the sources 
a and cXi are linear combinations of the contravariant components T^^ of the stress-energy 
tensor of the matter. Let us underline that the factor —Ai^G in these equations is a choice. 
We could of course introduce here a functional dependence on the spatial dimension d, for 
instance by replacing the factor Att by the surface of the unit [d — l)-dimensional sphere 
[see Eq. flB3|) in Appendix^, but this would only complicate the intermediate expressions 
without changing our final result. The matter sources a and cxj defined by Eqs. (|2.6p . (|2.7|) 
read (in d spatial dimensions): 



a = 



2 {d- 2)ro° + T 



d 



a 



(2.8) 



The definition for cxjj has been added for future use. Note that cxj and aij take the same forms 



as usual in 3 dimensions (see Eqs. (3.9) of Ref. |22]), but that the definition of a involves 
an explicit dependence on d. Conversely, the first and third of these equations allow us to 
express in terms of the above matter sources: = ac^ — cTjj) /{d — 2). A simple 
consequence of the expression of a is that the (i-dimensional Newtonian potential generated 
by a mass rua located at ja reads explicitly 



Uai^) = K(X) + 0{- 



2 ^ ~ ^ k^^"" 



d 



(2.9) 



where the factor 2{d — 2)/{d — 1) comes from a {i.e., from Einstein's equations), while the 
factor k = r(^^) /vr^" comes from the expression of the Green function of the Laplacian 
in d dimensions (see Eq. fl4.12|l below and Appendix lB|) . 

We give below the simplest forms of the metric and of the potential equations that we 
could obtain. We will explain afterwards which rules we followed to simplify them. Let us 
first define the useful combination 



V = V 



d-3 
d-2 



AX 16T 



Then the metric components can be written in a rather compact form: 



(2.10) 



goo 



.g (d-2)c^ 



1 - 



32 A K- 



(^-3)^ AVi 



2v I ^ 4 16 

^2 



^6 



1 + 



+ 



1 /d-1 V 



2\d-2 c2 



8Ri 16 
H H 



2{d-2] 



SijVkVk 



o 



+ 



(2.11a) 



(2.11b) 
(2.11c) 



The various elementary potentials V, Vi, K, Wij, Ri, X, Zij, Yi and T introduced in these 
definitions are ci- dimensional analogues of those used in Eqs. (3.24) of Ref. 0. Actually, 
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an extra potential is needed for (i 7^ 3, and it has been denoted K in Eq. ()2.1Up above. 
We give in Appendix ^ the exphcit expansion of this metric in powers of 1/c, as well as 
its inverse g^^ and its determinant g, which can be useful for future works. Note that the 
first post-Newtonian order of the spatial metric, gij = 5ij[l + pz|y^] + C'(l/c''), explicitly 
depends on d contrary to our choice ()2.6|1 for goQ. This dissymmetry between goo and gij is 
imposed by the field equations (j2.5p . 

The successive post-Newtonian truncations of the field equations ()2.2|) - ()2.5|1 give us the 
sources for these various potentials. The equations for DV and DVi have already been 
written in Eqs. ()2.7|) above. We get for the remaining potentials: 



aw, 



OR,, 



nx 



BY,. 



-AnGaV , 
-AttG ( cxij - 6ij 



4:7iG 



d-2 
d-1 
'd-2 

-AttG 



Va,, 



d 



dkVdiVk 



Va, 



+ 2 



2 

d{d-l) 
4(rf-2)2 

rf- 3 
d- 1 



1 /d-1 

2 V^^- 2 

dtVdiV . 



d^Vd.V, 



+W^jd^jV + 2VidtdiV + 



1 /d 



2 \d-2 

+^^2(9tVr-2d.V,d,V,, 
4:TtG _ / . cTkh \ d-1 



d-3 
d-2 

Vd^V 



d 



V OiA - 



d-2 



+dkVidkV,-2dkV(,d,)Vk 



5. 



d-2 



dkVm {dkVm — dmVk 



{d - 2)3 

1 fd-1 



-A-kG 



+ 



d-2 
1 



2{d-2) 



crRi 
CkkVi 



2{d 

{5-d){d 



o-VVi + - akWik 



d-3 



5-d 



K 



+WuiduV,-- 



1 fd-1 



d-2 

d{d 



2'^ikVk 



'J^2^'^^'^'~ A{d-2f 
2 



{d-2f 



VdfVdiV 



d-1 



d-2 

{d-l){d-3) 



vdkVdkV, + 



1 d-1 



2\d 



VdlV, + 2VkdkdtV., 



+ ' ' d,Kd,V, + 3) ^^^^^^^^ ^ d.VdtK) 



{d-2y 



4{d-2y 



(2.12a) 
(2.12b) 



(2.12c) 



(2.12d) 



(2.12e) 



(2.12f) 
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□T 



-47rG 



.2(^-1) 



5 — d 1 
(^ijW^j + — — — V^au + — r (rV^Vi 



4(rf-2)2 



1 /d-3 

2 W- 2 



1 



12 

d-3 ^ 



d-l 

-2diV, djRi - diV, dtWi, + 





1 

" 2 


d- 


3 


2{d- 


2)2 




/d 


^ij + 2 


\d 



-2 

- 1 



-^^+ 2(d-l)(rf-2V -^'^ 



d-2 



A{d-2y 



d~ 1 



8 \d-2 



2 



^idV? (d-l){d-3) rf(rf-l)(rf-3) 
"2^^*^^^ ~ 4(d-2)2 4(d-2)3 



(c;-i)(rf-3) d-3 lfd-3 



(2.12g) 



In Eq. (I2.12ep . parentheses around indices mean their symmetrization, i.e., a^ij) = \ {aij + 
aji). ^OY d = 3, the above set of equations ()2.12|) reduces to Eqs. (3.26) and (3.27) of 



Ref. |22]. The order of the terms and their writing has been chosen to be as close as possible 
to this reference. 

The harmonic gauge conditions (j2.ip impose the following differential identities between 
the potentials: 







dt 



1 fd 



2 \d-2 

'd-l 



V + 



2c2 



W + 



d 



d-2 



2{d-l){d-3) 



K 



2 

+ ^ 



\d-2 



X + ^ Z - ^ VkVu + ]-VW 
d-l d-l 2 



+ 



6{d-2) 



+dAVi + 



+ 



Ri + 



1 /d-l 



d-2 
1 



{d- l){d-3) 
{d-2f 



VK 



Y^--W,,V, + -WV + 



1 fd-1 



dt{V + - 



{d-l){d-3) 
2{d-2f 

4 



KV, 

VV 
1 



2 \d-2 



O 



VR^ + 



1 fd-1 



4 \ d-2 



Vv. 



(2.13a) 



% " 2 ^^''i 



(2.13b) 



where W = Wkk and Z = Zj^^ denote the traces of potentials Wij and Zij. For d = 3, these 
identities reduce to Eqs. (3.28) of Ref. jl^]. In this paper we shall check (see Sections HvXl 
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and IVI A|) that all the dimensionally-regularized potentials we use obey, at the indicated 
accuracy, the differential identities ()2.13p equivalent to the harmonic gauge conditions. 

In order to simplify as much as possible the above equations ()2.12|) for the potentials, we 
used the following rules: 

(i) We used the harmonic gauge condition ()2.13ap to replace everywhere diVi in terms 
of dtV and higher post-Newtonian order terms, and the gauge condition ()2.13b|) to 
replace djWij in terms of diW and dtVi [our knowledge of the higher order terms 
0{l/c^) in Eq. ()2.13b|) was actually not necessary for the simplification of Eqs. ()2.12|) ]. 
We also used the lowest order terms of Eqs. ()2.13|) to simplify their own higher order 
contributions. 

(ii) If the source of a potential P contained a double (contracted) gradient of the form 
□P = dkAdkB + (other terms), where A and B were two lower-order potentials, we 
got rid of the double gradient by defining another potential P' = P — \AB. We could 
then write its equation as UP' = -^{nA)B - lA{DB) + ^ dtAdtB + (other terms), 
in which \3A and \3B were replaced by their own explicit sources. The contribution 
proportional to was then transferred into the source of a higher order potential. 

(iii) At order C(l/c^), equation (j2.2j) - (j2.5|) for i?oo (^-c, for 0(700) contains the term 
WijdijV, that we introduced in the source of potential X, Eq. ()2.12d|) . In all other 
equations involving the same source Wij dijV, we used OX to eliminate it, instead of 
reintroducing it in the sources of other potentials. This is the reason why X is involved 
in the spatial metric Qij too at order 0{l/c^) [via the exponential of V in Eq. ()2.11c|) ]. 
and why VX appears again in ^fQo at order 0(l/c^). See the expanded form of the 
metric ()A1|) in Appendix [X] 

(iv) In the equation for Rqq at order 0{l/c^), we chose to eliminate a source propor- 
tional to VdiVj diVj, by including an all- integrated term VViVi in the definition of (700, 
Eq. ()2.11a|) . On the other hand, we could not eliminate at the same time the source 
term proportional to VidjVidjV in □(700, although it involves a double (contracted) 
gradient too. This is the reason why such a term appears in Eq. ( |2.12g| ). 

The above simplification rules have been applied systematically with a single excep- 
tion. Indeed, Eq. ()2.12ep for DZy involves double (contracted) gradients dkVidkVj and 
^ijdkVmdkVm- Therefore, the application of rule (ii) would have yielded another potential 

% = 4- - \Vry, + ^^^^ 5,,VuVk , (2.14) 

such that no double gradient appears in its source (but extra compact sources (T(iVj) and 
^ijCTfcVfc would have been involved). Although this modified potential Zij actually simplifies 
slightly some equations (but not all of them), we have chosen to use Zj ^ which is the direct d- 
dimensional analogue of the potential written in Eq. (3.27c) of Ref. |22|. Indeed, as explained 
in the following sections, the 3-dimensional results of this reference will be necessary for our 
(i-dimensional calculations, and it is more convenient to keep the same notation. 

Notice also that after the above simplifications, the resulting metric involves only poten- 
tials which are at most cubically non-linear (like for the term Wij dijV in the potential X — 
using the terminology of Section IV CI below) . There is no need to introduce any quartically 
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non-linear elementary potential because it turns out that it is possible to "integrate directly" 
all of them (at the 3PN level) in terms of other potentials. The only quartic contributions 
are the terms composed of l^^ and VX in the metric component goo [see Eq. ()Ala|) in Ap- 
pendix^. The fact that there are no intrinsically quartic potentials at the 3PN order made 
the closed-form calculation in 21, i^l possible. We shall comment more on this interesting 



fact in Section HV Al 

Let us now apply the general potential parametrization of the metric defined above to 
the specific case of (monopolar) point particles, i.e., to the action 



S = I ^-^V^Y^Ri9)-J2'^aC J ^-g,M)dyUy^^. (2.15) 



The stress-energy tensor T^'^(x) = [2c/ ^y—g{x)]6 Sraattcr/^ 9iJ.u{x) deduced from this action 
reads 

T^'^lx) = J2mac' jdsa^^^^ [-9{yaT'^ S^'^'\x' - y',{s,)) , (2.16) 

where dsa = \/ —gp^v (z/a) dya dy^ is (c times) the proper time along the world-line of the 
a^^ particle and where is the Dirac density in ci -|- 1 dimensions (J d^~^^x 5^'^'^^Hx) = 



1). Here, we take advantage of the fact (emphasized in [2^) that dim. reg. respects 
the basic properties of the algebraic and differential calculus: associativity, commutativity 
and distributivity of point-wise addition and multiplication, Leibniz's rule, Schwarz's rule 
{dfj, dy f = dy dfj, /), integration by parts, etc. In addition, the post-Newtonian expansion of 
9fiu{x) yields "(i-dimensional functions" which are formally as smooth as wished (by taking 
the real part of d small enough) in the vicinity of the world-lines: see for instance Eq. ()2.9|) . 
This allows us to work with self-gravitating point particles essentially as if they were test 
particles. For instance, we can use F[g^y{x)] S^'^^ (x — y^) = F[g^y(ya)] 5^'^^ (x — y^). In 
particular, the ?/a-evaluated determinant factor [—g{ya)]^^ in (!2.16|) came from the field- 
point dependent factor [—g{x)]~2 in the definition of T^^{x). Similarly, the usual derivation 
of the equations of motion of a test particle formally generalizes to the case of self-gravitating 
point particles in d dimensions. One then finds that the equations of motion of point particles 
can equivalently be written as 

V^T^'^(x)=0, (2.17) 

or as the usual geodesic equations. The latter can be written either in covariant form 
< = « = dy^/dsa), I.e., explicitly 

^^TUMiy.)]'£f^^o, (2,18) 

where F^j, = \g^'^{d^gy^ + d^g^^ — df^g^y) as usual, or in the explicit form corresponding 
to using the coordinate time t = ya/c as parameter along the world-lines, which is easily 
derived from the covariant expression with a lower index, Ma^i''^= d^g^pU^^ / dsa = 
\(^ii9vp'Ki'^a- Like in 3 dimensions, cf. Eqs. (3.32)-(3.33) of Ref. [2j|, it can thus be put in 
the form 

- = F-. (2,19) 
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where 



~ 2 yZ-QpaVpy'/c^ 



= dx^/dt 



, , . _ ^ ^ , (2.20) 

(c, v) denoting the coordinate velocity. Let us emphasize again that in d 



dimensions, all the non-linear functions of gfiuiVa) and OxgiMuiUa) that will enter our calcu- 
lation of (j2.19j) - (j2.2m) can be treated as in the x ^ evaluation of smooth functions of x. 
For instance, denoting for simplicity / = 2{d — 2)/{d — 1), the Newtonian approximation, 
say f/*^'^^(x) = t/(x), of the basic scalar potential V{'x), reads, in the vicinity of x = yi. 



f/(x) = fkGmirl-'^ + U2ix) 



(2.21) 



where ?72(x) = f kG 1712 is (in any (£) an indefinitely different iable function of x near yi. 
Analytically continuing d to sufficiently "low" (and even with negative real part, if needed) 
values, we see not only that [f/(x)]x=yi = ?72(yi), but that [t/"(x)]x=yi = (?72(yi))", and, 
e.g., [f/^'(x)9,f/(x)]x=y, = (f/2(yi))^9,f/2(yi), etc. 

Although the expressions (j2.2(jp do not depend explicitly on the dimension ci, the met- 
ric (j2.1H) does depend on it, and therefore the post-Newtonian expansions of Eqs. (j2.2(jp 
involve many rf-dependent coefficients. We give their full expressions in Appendix 1X1 
Eqs. (IA11|) - ()A12|) . but we quote below only their Newtonian orders and the very few terms 
which will contribute to the poles oc l/(c? — 3) in our dimensionally regularized calculations: 



4 16 " /I 



c" V c° 

The acceleration & = dw/dt can thus be written as 
d (P* - v') 



F' - 



diy + - 



dt 



4 



-t'V ^ dRj 

diX + 2—^ + 
dt 



16 
+ ^ 



9.r + ^ + 

dt 



+ o[- 



(2.22) 



(2.23) 



In Section FlV Al we shall give flesh to the formal expressions written above by explaining by 
which algorithm one can compute, with the required accuracy, the explicit ci-dimensional 
expansions near x = yi [analogous to the simple case (j2.21|) ] of the various elementary 
potentials entering Eq. ()2.23|1 . and notably of the crucial ones X, Ri (to be computed with 
IPN accuracy) and T, (to be computed at Newtonian order only). 



III. HADAMARD SELF-FIELD REGULARIZATIONS IN 3 DIMENSIONS 

The main aim of this Section is to complete the first step of the strategy outlined in the 
Introduction, i.e. to collect a complete list of the additional contributions to the equations 
of motion which are specific consequences of the use of the extended Hadamard regulariza- 
tion methods defined in ji^ |2^. However, to do that we need to start by recalling some 
material concerning the Hadamard regularization in 3 dimensions, and by contrasting it 
with dimensional regularization. Such material is needed for understanding our computa- 
tion based on the "difference" in Section IIVI We shall start by recalling the definition of 
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the "ordinary" Hadamard regularization and complete it by defining what we shall call the 
"pure" Hadamard regularization. Then we shall recall the main new features of the extended 
Hadamard regularization defined in jisl . 2^1, and collect the additional contributions to the 



equations of motion which are specific consequences of the use of the extended Hadamard 
regularization (there are seven such additional contributions). 



A. Ordinary and "pure" Hadamard regularizations 

The phrase "Hadamard regularization" covers two distinct concepts: (i) the regulariza- 
tion of the "limit" limx^yi F(x;yi,y2) where F(x;yi,y2) belongs to a class of singular 
functions (generated by the iteration of Einstein's equations), and (ii) the regularization of 
the 3-dimensional integral J (i^xF(x; yi, y2) of some function F & The class of functions 
T consists of all functions -F(x) on that are smooth except at yi and y2, around which 
they admit Laurent-type expansions in powers of ri or r2 (see Section II. A of [23] for the 
precise definition of JF). When ri = |x — yi| — > {i.e., around singularity 1) we have, 
VA^ G N, 

^(x)= E r?/,(nO + o(rn, (3.1) 

PQ<p<N 

where the Landau o-symbol takes its usual meaning, and the i/p(ni)'s denote the coefficients 
of the various powers of ri, which are functions of the positions and velocities of the particles, 
and of the unit direction rii = (x — yi)/ri of approach to singularity 1. The powers of ri 
are relative integers, p G Z, bounded from below by some typically negative po depending 
on the F in question. 

The Hadamard ^^partie finie" of the singular function F at the location of the singular 
point 1 (first meaning of Hadamard regularization) is defined as the angular average of the 
zeroth-order coefficient in the expansion (|3.ip . It is denoted (F)i, so that 

m.-(/.>-/^/o(nO, (3.2) 

where dfl{ni) denotes the usual surface element on the 2-dimensional sphere centered on 1. 
We shall employ systematically a bracket notation Q for the angular average of a function 
of the angles (either ni or 112).^^ A distinctive feature of the Hadamard partie finie ()3.2|) is 
its "non-distributivity" in the sense that 

{FG)i ^ (F)i(G)i in general for F,GeJ^. (3.3) 

The non-distributivity represents a crucial departure away from the simple algebraic prop- 
erties of the analog of {F)i in dim. reg. which is merely F*^'^^(yi). It is an interesting fact 
that in a post-Newtonian expansion the non-distributivity starts playing a role only at the 
3PN order (because the functions there become singular enough). Up to the 2PN order one 
can show that {FG)i = (F)i (G)) for all the functions involved in the equations of motion in 
harmonic coordinates [13| . Several of the problems of the Hadamard self-field regularization 



Since this will always be clear from the context, we do not specify on the brackets if the angular integration 
should be performed around point 1 or 2. We do not indicate either if the integration sphere is two- 
dimensional or {d — l)-dimensional (as we shall see later there can be no confusion about this). 
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(in the "ordinary" sense) when apphed at the 3PN level {e.g. the occurrence of the unknown 
constant A) are related to the latter non-distributivity. 

The second notion of Hadamard partie finie (denoted Pf in the following) is to give 
a meaning to the generally divergent integral J (i^xF(x). In this work we shall have to 
consider only the ultra-violet (UV) divergencies of the integrals, i.e., at the locations of the 
two local singularities yi and y2. All functions involved at the 3PN order are such that 
there are no infra-red (IR) divergencies when |x| — oo (this is true not only in 3 dimensions 
but also for any dimension d = 3 + e in a. neighborhood oi d = 3). The Hadamard partie 
finie of the (UV) divergencies is then defined as 



d3xF(x) 



Pf,,,, / d='xF(x) = lim( / 

+4^ E ^(/p>+4vrlnf-V/-3> + l-2|. (3.4) 

The description of this formula in words is as follows. One first excises two spherical balls 
Bi{s) and B2{s) surrounding the two singularities (each one having the same radius s), and 
one computes the integral on the volume external to these balls, i.e., R'^ \ -Bi(s) U B2{s) — 
first term in Eq. ()3.4|) . That integral tends to infinity when s — 0, but we can subtract from 
it its purely divergent part, which is given by the additional terms in (j3.4|) (which obviously 
are to be duplicated when there are 2 singularities, cf. the symbol 1^2). The limit s ^ 
then exists (by definition) and defines Hadamard's partie finie. 

Notice the crucial dependence of the partie finie on two constants Si and S2 entering 
the log-terms. These constants have the dimension of length. We shall say that si is the 
regularization length scale associated to the Hadamard regularization of the divergencies 
near x = yi (similarly for S2). Note also that the Hadamard partie finie does not depend 
(modulo changing the values of si and S2) on the shape of the regularization volumes Bi 
and i?2, above chosen as simple spherical balls (see the discussion in Ref. jl^]). 

An important consequence of the definition (j3.4j) is that, in general, the integral of a 
gradient diF is not zero, because the surface integrals surrounding the singularities become 
infinite when the surface areas tend to zero, and may possess a finite part. We find (see 
Eq. (3.4) in 0) 

Pf j rf^x9iF(x) = -47r(nl/_2> + 1 ^ 2. (3.5) 

For a general F ^ T the R.H.S. is typically non-zero. This fact shows that the application 
of the ordinary Hadamard regularization in the post-Newtonian iteration has to be supple- 
mented by a notion of distributional derivatives, in order to ensure that the integrals of 
gradients are zero like in the case of regular functions. Notice that the constants s\ and 
disappear from the result ()3.5p . [We shall also see the need, within dim. reg., to consider 
some derivatives in the sense of distribution theory.] 

Let us apply the definition ()3.4|) to the integral of a compact-support or "contact" term, 
i.e., made of the product of some F and a Dirac delta-function at the point 1. Let us 
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formally assume that 

Pf j d3xF(x)5(3)(x-yi) = (3.6) 

which is the most natural way, within Hadamard's regularization, to give a sense to such 
integral. Now the problem with that definition is that if we want to dispose of some lo- 
cal meaning (at any field point x) for the product of F with the delta-function, then as 
a consequence of the non-distributivity we cannot simply equate F{x.)5^^\x — yi) with 
(F)i5(3)(x_y,), ^.e.: 

F(x)5(^)(x - yi) ^ (F)i(5(^)(x - yi) in general for F e (3.7) 

Indeed, if it were true that F6f^^ = {F)i6^^ [for simplicity we denote 6^^^ = — yi)], 

then multiplying by any G we would have FG6l = {F)iG6\ , and by integrating over 

following the rule ()3.6|) this would yield {FG)i = {F)i{G)i in contradiction with the 
violation of distributivity ()3.3|) . Therefore, both the violation of distributivity ()3.3|) and the 
consequence (j3.7p are unescapable in the ordinary Hadamard regularization. 

The previous situation should be contrasted with the d-dimensional case for which the 
distributivity is always satisfied, as we have simply 

(i.WG'W)(y^)=i.(^)(y^)G'W(y^), (3.8) 

and 

F(^)(x) 5W(x - yi) = F('^)(yO 5(<^)(x - y^) . (3.9) 

Finally, taking the d ^ 3 limit, we see that the dim. reg. way of regularizing a three- 
dimensional "contact term", i.e. a term like F{x.)5^^\x. — yi), is by considering it as the 
d — > 3 limit of its rf-dimensional analogue ()3.9p . That is, 

dim. reg. [F(x) _ y^^j ^ (\imF^'^\yi)^6^^\^~yi), (3.10) 

where F^'^'^ is the d-dimensional version of F, as obtained by solving Einstein's equations 
in d dimensions (using the method explained in Section fl V Al below) . There are no poles 
in the calculation of the "contact" terms in any of the potentials at the 3PN order so the 
limit d — i> 3 in Eq. ()3.10p always exists. Once again the dim. reg. prescription ()3.10|) owns 
all the good features one wishes, notably the distributivity as we have emphasized in Eqs. 
(|IHl)-(Ea)- 

In the following it will be convenient, in order to compare the present dim.-reg. calcu- 
lation with the Hadamard-based work 2^\, to introduce the terminology pure Hadamard 



regularization to refer to the following "minimal" version of the Hadamard regularization: 
(a) an integral J (i^xF(x), where F is made of some product of derivatives of the non-linear 
potentials V, Vi, ■ ■ ■ , is regularized by the ordinary Hadamard partie finie prescription ()3.4p . 
without bringing in any distributional contributions (see below for the treatment of these); 



Actually this assumption should be viewed as the definition of a new object we can call Pf(5'^^(x — yi) 
and which takes the property . This is exactly what we do in the context of the extended Hadamard 
regularization. 
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(b) the regularization of a product of potentials V, Vi, Wij, ■ ■ ■ (and their gradients) is 
assumed to be distributive, which means that the value at the singular point yi of some 
polynomial in V, Vi, Wij, ■ ■ ■ and their gradients, say J-'[V, Vi, Wij, diV, ■ ■ ■], is given by the 
replacement rule 

(^[y, V, d,v, ■■■]), — my\, {v\, {d,v)i, ■■■]■ (3.11) 

and (c) a contact term, i.e. of the form F(x)(5^^''(x — yi), appearing in the calculation of 
the sources of the non-linear potentials, is regularized by using the rule 

J^[V, V, ■ ■ ■ ](5(3)(x - yi) T[{y\, iVi)^, {Wi,\, ■ ■ ■ ](5(3)(x - yO , (3.12) 

(there are no gradients of potentials in the contact terms). The rules (j3.1H) - (j3.12|) of the 
pure-Hadamard regularization are formally equivalent to assuming the replacement rules 
{FG)i — > {F)i{G)i together with F(x)(5(3)(x - yi) — ^ (F)i5(3)(x - yi), in the case where 
F and G are made of products of our elementary potentials and their gradients. The rules 
of the pure Hadamard regularization are, however, well defined, and are not submitted (by 
their very definition) to the consequences of the ordinary Hadamard regularization ()3.3p 
and J3.7p. Note also that, as done in previous computations of the 3PN ADM-Hamiltonian 
[l6llT7j and the 3PN binary's energy flux 27|, one can formally use fl3.11|) - ()3.12|) at the price 



of adding a limited number of arbitrary parameters (considered as unknown). 

The definition (j3.12j) of pure-Hadamard regularization for contact terms is useful because 
we have checked that, when using the dim. reg. prescription (j3.1(J|) (in the limit where 
d — >• 3), all the contact terms in the sources of the non-linear potentials V, Vi, Wij, ■■■ 
needed at the 3PN order, agree with the result of the pure-Hadamard regularization. [Of 
course, we would not need to introduce a notion of pure-Hadamard regularization in a direct 
calculation of the equations of motion in d dimensions, i.e., not based on the "difference" 
between Hadamard and dim. reg., because in such a pure dim. reg. approach the contact 
terms would be treated unambiguously from the start using Eq. ()3.10|) .] On the other 
hand, when computing the value at the singular point of the potentials for insertion into 
the geodesic equations, we do find some departure between the dim.-reg. calculation and 
the (ordinary or extended) Hadamard one. Let us illustrate these differences by means of 
the simplest example which does enter our 3PN calculation, namely the regularization of 
{UYdiU where U is the Newtonian potential. In d dimensions U^'^^{x) is given by Eq. (|2.2ip 
[we add here a superscript (d) to indicate the ^-dimensionality of a potential and pose 
U = U^^^. Therefore in dim. reg. the result is simply 

lim([[/(^)(yi)]^9,f/('^)(yi)) = [U2{yi)]%UM (dim. reg.) , (3.13) 

where ?72(yi) = Gm2/ri2 is the value at point 1 of the potential of the other particle. The 
result ()3.13p is the same as when using the pure Hadamard regularization. Indeed, we find 
first that {U)i = f/2(yi) and (9jf/)i = 9jf/2(yi), and then, by using the definition ()3.11|) . 

(U'd.U), ^ [{U)if{diU), = [f/2(yi)]'9,f/2(yi) (pure-Hadamard) . (3.14) 



12 Thus we shall write {VW^J)l — > {V)i{W,j)i or {V^d^V)i — > [{V)i]^{d^V)i, but not, for instance, 
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On the other hand, the latter results contrast with the application of the ordinary Hadamard 
regularization for which we find 

{U^d,U\ = [t/2(yi)]'9,f/2(yi) + ^[f/i(y2)]'f/2(yi)5.f/2(yi) (ordinary-Hadamard) . (3.15) 

The first term is in fact the "pure-Hadamard" result which is in agreement with the dim. 
reg. one. The second term is an example of the non-distributivity of the ordinary Hadamard 
regularization/^ which is also systematically taken into account in the extended-Hadamard 
regularization that we shall describe in Section llll CI 



B. Ordinary Hadamard regularization of three-dimensional Poisson integrals 

Let us give some reminders of the way we apply the considerations of Section Qll Al to the 
computation of Hadamard-regularized potentials having the form of Poisson or Poisson-like 
integrals. Let us first discuss the prescription one has taken in d = 3 to define the "value 
at x' = yi" of a (singular) Poisson potential -P(x'). In ci = 3, the Poisson integral -P(x'), at 
some field point x' G M^, of some singular source function -F(x) in the class is defined in 
the sense of the partie-finie integral ()3.4j) . namely 

^(xO = -i^Pfs... / ^%^(x) , (3.16) 
47r J |x — x'l 

where Si and S2 are the two constants introduced in Eq. (|3.4j) . At first sight we could think 
that a good choice for defining the pure Hadamard value [P(x')]x'=yi is simply to replace 
x' = yi in ()3.16|) . i.e., 

P(yi)--^Pf.,./ ^i^(x). (3.17) 



However, the work on the 3PN equations of motion |21l . |22| suggested that the definition 
()3.17|) is not acceptable: it did not seem to be able to yield equations of motion compatible 
with basic physical properties such as energy conservation. 

The choice adopted in j2l|, i^l is to define the regularized "value at x' = yi" of the 



function P(x') by taking the Hadamard partie finie in the singular limit x' yi. Notice 
first that -P(x') does not belong (in general) to the class T because the Poisson integral will 
generate some logarithms of in its expansion when 0. Thus, we shall have, rather 

than an expansion of type ()3.1|) . 



p(x')= ^1 



IP 



^p(n;) + /ip(n;)lnr; 



+ o(r;^), (3.18) 



^•^ In the ADM-Hamiltonian the analogue of this example is the regularization of C/^, which gives automati- 
cally [/72(yi)]'* in dim. reg. and also (by definition) in the pure-Hadamard reg., while 

The latter example represents in fact the only source of ambiguity present in (the static part of) the 
ADM-Hamiltonian formalism l35l. 
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where the coefficients iQp and ihp depend on the angles n'^, and also on the constants si 
and S2, in such a way that when combining together the terms in ()3.18|) the constant r[ 
always appears in "adimensionalized" form like in ln(r^/si). Then we define the Hadamard 
partie finie at point 1 exactly in the same way as in Eq. (j3.2|) . except that we now include 
a contribution linked to the (divergent) logarithm of r[, which is possibly present into the 
zeroth-order power of r[. More precisely, we define 

{P), = {g,) + {h,)\nr[, (3.19) 

where we introduced a new regularization length scale denoted r^, which can be seen as 
some "small" but finite cut-off length scale [so that In r[ in Eq. ()3.19|) is a finite, but "large" 
cut-off dependent contribution]. We shall see later that the dependence on r'^ disappears 
(as it should) when adding to (P)i the difference DP(1) = -P^'^^(yi) — (-P)i- To compute 
the partie finie one must apply the definition ()3.19p to the Poisson integral ()3.16p . which 
involves evaluating correctly the angular integration therein. The result, proved in Theorem 
3 of pH, is 



;/-2>. (3.20) 



We recover in the first term the value of the potential at the point 1: P(yi), given by 
Eq. ()3.17|) . The supplementary term makes the partie finie to differ from the "naiv e" g uess 
P(yi) in a way which was found to play a significant role in the computations of j2lL 2^ . 



The apparent dependence of the result ()3.2()|1 on the scale si is illusory. The Si-dependence 
of the R.H.S. of Eq. ()3.2()|1 cancels between the first and the second term, so the result 
depends only on the constants r[ and S2, and we have in fact the following simpler rewriting 
of (IT^ . 

1 f rf^x 



(P)i--^Pf<..y —P(x) -(/_,). (3.21) 

Similarly the regularization performed at point 2 will depend on r'2 and Si, so that the 
binary's point-particle dynamics in Hadamard's regularization depends on four (a priori 
independent) length scales r^, S2 and rg, si. The explicit expression of the result ()3.2H) is 
readily obtained from the definition of the partie-finie integral ()3.4|) . We find (see the details 
in Ref. Q) 



(P)i = lim< 



1 



\Bi{s)UB2{s) ^1 



P(x) 



p+2<0 
+00 , y 

-Ol I — 



+ 1 



t=0 



i\ Vri2 



.p+^+3<0 



,?2n 



S2 



(3.22) 



Note that the terms corresponding to singularity 2 involve the multipolar expansion around 
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the point y2 of the factor 1/ri = l/|x — yi| present into the integrand. 

Because we work at the level of the equations of motion, many of the terms we shall need 
in this paper are in the form of the gradient of a Poisson-like potential. For the gradient we 
have a formula analogous to (|3.2ip and given by Eq. (5.17a) of namely 



{n\f-i) 



(3.23) 



where we have taken into account (in the rewriting of the second line) the always correct 
fact that the constant Si cancels out and gets "replaced" by r[. Notice that in ()3.23p there 
is no additional term to the partie finie integral similar to the last term in ()3.21|) . The 
corresponding exphcit expression is 



(9,P)i = lim< 



1 

4:71 



R3\Bi(s)UB2(s) 



rf^x ^F(x) 



„, -diL — 
\ri2 



.p+l+3<0 




Finally we must also treat the more general case of potentials in the form of retarded 
integrals [see Eqs. ()2.12p ]. but because we shall have to consider (in Section |IVB|) only the 
difference between the dimensional and Hadamard regularizations, it will turn out that in 
fact the first-order retardation (IPN relative order) is sufficient for this purpose. Actually, in 
this paper we are not interested in radiation-reaction effects, so we shall use the symmetric 
(half-retarded plus half-advanced) integral. At the IPN order we thus have to evaluate 

i?(x') = P(x') + ^Q(x') +o(^^y (3.25) 

where -P(x') is given by ()3.1(ij) . and where Q{'x.') denotes (two times) the double or "twice- 
iterated" Poisson integral of the second-time derivative, still endowed with a prescription of 



We write the multipole expansion in the form 

employing cm' usual notation where capital letters denote multi-indices: L = iii2 ■ ■ ■ ii, and, for instance, 
= n^2 ' ' ' "-2*- The expansion is symmetric-trace-free (STF) because Si^i^_-^dL{l/ri2) = 9L-2A(l/ri2) = 
0. Here 9L(l/ri2) is a short-hand for £ partial derivatives d/dy\2 of l/ri2 = l/|yi — y2|- The multipole 
expansion in d dimensions (also STF) is given by Eq. H4.23|l below. 
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taking the Hadamard partie finie, namely 

g(x') = -^Pfs„s. j d'^ |x - x'|92F(x) . (3.26) 

In the case of Qjx.') the results concerning the partie finie at point 1 were given by Eqs. (5.16) 
and (5.17b) of 0, 

iQ)i = -^^k,s, j d'^nd^Fi^) + ^{k -4> , (3.27a) 
{^,Q)^ = ^Piri,s,J d^^n\d^F{^) + ^{n\k^s), (3.27b) 

where the ifc^'s denote the analogues of the coefficients ifp, parametrizing the expansion of 
F when r[ 0, but corresponding to the double time-derivative d^F instead of F. [In the 
following we shall not need the explicit forms of the results ()3.27p .] 

Let us clarify an important point concerning the treatment of the repeated time derivative 
d'^F[x.) in Eqs. ()3.27j) . As we are talking here about Hadamard- regularized integrals (which 
excise small balls around both yi and y2), the value of dfF{x) can be simply taken in the 
sense of ordinary functions, i.e., without including eventual "distributional" contributions 
proportional to (5(x — yi) or (5(x — y2) and their derivatives. However, we know that such 
terms are necessary for the consistency of the calculation. This is why we must also include 
somewhere in our formalism the difference between the evaluation of these distributional 
terms in d dimensions, and the specific distributional contributions issued from the general- 
ized framework used in ^2^. This difference will be included in Section Fill Dl below, among 
the complete list of additional contributions specifically related to the use of the extended 
regularization approach we shall now describe. 



C. Extended-Hadamard regularization 

The extended- Hadamard'^ regularization, proposed in Refs. [i^ l^l, tackles the partic- 
ular properties of the ordinary Hadamard regularization, notably the non-distributivity of 
Eqs. (j3.3p and (|3.7|) . and the fact that the integral of a gradient is not zero [Eq. (j3.5p ]. 
These properties are implemented within a theory of pseudo-functions, viz linear forms de- 
fined on the set of singular functions JF. The use of pseudo-functions in this context enables 
one to give a precise meaning to the object F5(x — yi) needed in the computation of the 
contact terms, and which is otherwise ill-defined in distribution theory. Furthermore the 
use of some generalized versions of distributional derivatives permits a systematic treatment 
of integrals and a natural implementation of the property that the integral of a gradient 
is always zero. In this paper we shall content ourselves with recalling the principle of the 
extended-Hadamard regularization, and with presenting its "ready-to-use" consequences. 

To any F G we associate the "partie finie" pseudo-function PfF, which is the linear 
form on JF defined by the duality bracket 

< PfF, G >= Pf y rf^x F(x)G(x) , (3.28) 

which means that the action of PfF on any G G JF is the partie-finie integral, as given by 
f|3.4|) . of the ordinary product. The pseudo-function PfF reduces to a distribution in the 
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ordinary sense of Schwartz j36| when restricted to the usual set T) of smooth functions with 
compact support on M^. The product of pseudo-functions coincides, by definition, with the 
ordinary point- wise product, namely PfF.PfG = Pf(FG). In the class of pseudo-functions 
constructed in Ref. j23|, the "Dirac-delta" pseudo-function ViSi is defined by 

< Pf(5i, F>=Pfj d^x 5i(x)F(x) = , (3.29) 

where {F)i denotes Hadamard's partie finie ()3.2j) . This definition, which obviously yields a 
natural extension of the Dirac function 5i(x) = 5*^'^^(x — yi) in the context of Hadamard's 
regularization, leads also to new objects which have no equivalent in distribution theory, 
the most important one being the pseudo-function Pf{F6i) which played a crucial role in 
21 . 2^ 1 for the calculation of the compact-support parts of potentials as well as the purely 
distributional parts of derivatives. It is given by 

< Pf(F(5i), G >= {FG)i , (3.30) 

where one should be reminded that it is in general not allowed to replace the R.H.S. by the 
product of regularizations: (FG)t^ ^_(F)i{G)i. 

In the actual computation [2l|, |2^ the pseudo-function Ff{F6i) acts always on smooth 
functions with compact support (e T>), in which case it reduces to a distribution in the 
ordinary sense, which was shown to admit the "intrinsic" form 

Pf(F(5i) = J2 ^ri^i f -e)dLSi (when restricted to V) . (3.31) 

Here L = denotes a multi-index composed of C multipolar indices ii, ■ ■ ■ 

Ol = ■ ■ ■ di,, means a product of i partial derivatives di = d/dx^, and raf = n^i ■ ■ -rii 
a product of i unit vectors (we do not write the i summation symbols, from 1 to 3, over the 
indices composing L). Notice that the sum in Eq. 1)3.3111 is finite because F admits some 
maximal order of divergency when ri — >■ 0. Now we discover that the "monopole" term in 
the latter multipolar sum, having £ = 0, is nothing but {F)i6i which is exactly the result 
we would get following the pure-Hadamard regularization rule fl3.12|) . [Indeed, as we are 
considering here only the contact terms entering the source terms for the 3PN-level non- 
linear potentials, the "ordinary" Hadamard regularization {F)i coincides with the "pure" 
Hadamard regularization ()3.12|) .] The sum of the other terms then define what we can call 
some non-distributive contributions because their appearance is the direct consequence of 
the violation of distributivity ()3.3|) . Thus, 

Pf(F5i) = (F)i5i + "non-distributivity" contributions . (3.32) 

In the work care has been taken of all such non-distributivity terms. Consider for 
instance the Poisson integral of a compact-support term Pf(F5i) (say, proportional to the 
matter source densities a, ai or aij). Using (|3.3ip the Poisson integral reads 



(3.33) 
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To apply we assume that x' is distinct from the 2 singularities yi and y2; see |23| for more details 
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[where di should be better written 9^ (1/r'^)]. Evaluating now the partie finie ()3.2|) 

at both singular points [i.e., when — > and r'2 0] we obtain 

-Pf(F(5i) = , (3.34a) 



X — X' 



1 

00 



(/ (;!;) - J f «(^«- f^^-^) 

The result fj3.34ap is in agreement with the pure-Hadamard regularization; however 
Eq. ()3.34b|) does involve some extra terms with respect to the pure-Hadamard calculation, 
since the latter is easily seen to simply yield {F)i/ri2, which is nothing but the "monopolar" 
term £ = of the multipolar sum in the R.H.S. of ()3.34b|) . Therefore we decompose ()3.34b|) 
as 

Pf / -. tF6i ] = + "non-distributivity" contributions . (3.35) 

y |x-x'| J 2 ri2 

The second ingredient of the extended-Hadamard regularization concerns the treatment 
of partial derivatives in some (extended) distributional sense. Essentially, one requires [2^ 
that the derivative reduces to the ordinary derivative in the case of regular functions, and is 
such that one can integrate by parts any integrals. The latter property (valid for the spatial 
derivative) translates into 

< ^^(PfF), G >= - < di(PfG), F > . (3.36) 

This rule contains the standard definition of the distributional derivative |3^] as a particular 
case. It implies the important property that the integral of a divergence is zero. Let us pose 

9,(PfF) = Pf(9,F) + D,[F], (3.37) 

where Pf{diF) denotes the derivative of F viewed as an "ordinary" pseudo-function, and 
Dj[F] represents the purely distributional part of the spatial derivative (with support con- 
centrated on yi or y2). 

Looking for explicit solutions of the basic relation ()3.36|) we have found ji^], with the 
help of Eq. (ITHll : 



D,[F] =47rPfUl 



2 1 Ti 1 

fc>o 1 



61] +1^2. (3.38) 



Notice that Dj[F] depends only on the singular coefficients of F (coefficients of negative 
powers of ri in the expansion of F). The derivative operator defined by Eqs. ()3.37j) - ()3.38|l 
does not represent the unique solution of ()3.36j) . but it has been checked during the calcu- 
lation that using another possible solution results in physically equivalent equations of 



motion at 3PN order {i.e., reducing to each other by a gauge transformation). Concerning 
multiple derivatives we dispose of the general formula 

e 

Dl[F] = J29^^■■■^^-M^^,^.■■■^^F] , (3.39) 
k=l 
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giving the distributional term associated with the £-th spatial derivative, Dl[F] = diPfF — 
PfdiF (where L = 11X2 ■ ■ - ie), in terms of the single derivative Dj[F]. As an example, to 
treat the second-derivative of the Newtonian potential, dijU where U = Grrii/ri + Gm2/r2, 
one uses 

477 15 



1 

ri 



Pf S'^ 



-nv 



5i 



(3.40) 



where fi*/ = — ^Sij. Therefore the extended distributional derivative differs in general 
from the usual Schwartz's derivative [cf. the second term in ()3.40|) ]. [This is unavoidable 
if one wants to respect the basic rule of integration by parts fl3.36|) for general functions 
in the class !F.] Notice also that we do find a distributional term in the case of the first 
derivative: Dj [1/ri] = 27rPf {rin\6i). We recall also (for future use) the case of the partial 
time-derivative, dt{PiF) = Pi{dtF) + Dt[F], whose distributional term is given by (following 
Ref. 0) 

I),[F] = vlD ,[F] + v',D ,[F] , (3.41) 

in terms of the partial derivatives with respect to the source points yi and y2, namely iDjF] 
and 2Di[-F]- The explicit expression reads 



Dt[F] = -AnPi{in,v,) 



0^ f-i 
2 1 ' — ' r 

fc>0 



-r / -2-k 



5i 



(3.42) 



where [n\V\) denotes the ordinary scalar product [notice the over-all sign difference with 
respect to Eq. ()3.38|l ]. Multiple time-derivatives can be treated accordingly to Eq. ()3.39j) . 
For instance, 

D«[F] =D49tF]+9iDi[F]. (3.43) 



Following the regularization [23| all the distributional terms [of type Pf(-F5i) and Vi{F82)\ 
issued from the latter distributional derivatives are to be treated when computing the po- 
tentials according to the extended contact term definitions of Eqs. ()3.34a|) - ()3.34b|) . 

Finally let us turn to the extension of the Hadamard regularization (introduced in [j^ D 
concerning the definition of a new operation of regularization, denoted [-F]i, consisting of 
performing the Hadamard regularization (F)i within the spatial hypersurface that is geo- 
metrically orthogonal (in a Minkowskian sense) to the four- velocity of the particle 1. The 
regularization \F\x differs from (F)i by a series of relativistic corrections calculated in (2^ . 
Together with the other improvements of the extended-Hadamard regularization, it resulted 
in equations of motion in harmonic-coordinates which are manifestly Lorentz invariant at 
the 3PN order Here we give a formula, sufficient for the present purpose, for ex- 

pressing \F\\ in terms of the basic regularization (-F)i, defined by ()3.2j) . at the IPN order: 



dtF + -v\d,F 



(3.44) 



The first term is simply (-F)i, while the other terms define a set of relativistic corrections 
required for ensuring the Lorentz invariance of the final equations in Hadamard's regular- 
ization. Hence, we decompose (j3.44|) into 



'Lorentz" contributions . 



(3.45) 
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D. Contributions due to the extended-Hadamard regularization 



After the reminders of the last subsections, we are now in position to explain the origin 
of all the contributions [included in the final result ()1.5|) ] which were due to the specific 
use of the extended-Hadamard regularization. Actually, we shall list here the contributions 
due to the use of the full prescriptions of [IH 01 with respect to those that would follow 
from using what we shall call a ^^pure Hadamard-Schwartz" (pHS) regularization. By this 
we mean: (1) treating the contact terms of all the non-linear potentials V, Vi, 



as m 



(j3.12p [we have checked that for all the potentials involved this is equivalent to (j3.1(Jj) ]: (2) 
treating the distributional part of an integrand such as Fij dij U in the normal Schwartz's 
distributional way, for instance 



da I — 

ri 



yi; 



(3.46) 



and evaluating the contact term generated by the delta function in the "pure Hadamard" 
way (|3.12j) : (3) regularize any three-dimensional integral by the ordinary Hadamard prescrip- 
tion ()3.4|) : and, finally, (4) using systematically, in the last stage of the calculation where 
one replaces the metric into the geodesic equations, the pure Hadamard replacement rule 
appearing in 1)3.111) [for instance, we write (y^diV)i — > [(y)i]'^{diV)i, creating therefore a 
net difference with respect to the ordinary and/or extended Hadamard regularizations for 
which {V%V), ^ [{V)^]\d,V),]. 

The usefulness of the definition of such a pHS regularization is to "localize" the additional 
contributions brought by dim. reg. to the occurrence of poles (x 1/e (or "cancelled" poles) 
in d dimensions. 

Our complete list of additional contributions contains seven items. First of all there are 
four "non-distributivity" contributions of the type given by Eq. ()3.35p : 

(i) The so-called "self" terms, for which the delta-function in Pf{F6i) comes from the 
purely distributional part of the distributional derivative given by (|3.38|) . The self 
terms were derived in Eq. (6.20) in [22]; they read explicitly 



srself , 



151 G^mim: 



9 



n 



12 



+ 



+ - {ni2V2)v\ 
5 



(3.47) 



where {ni2V2) denotes the usual scalar product between ni2 and V2, and where = v|. 
The expression fl3.47|) can be rewritten in a simpler way as (where we denote for 
simplicity V2 = V2V2) 



self „i 



151 G^mim\ 



9 



n\2 + 



1 G^ 



m: 



12 



30 c6 



■V2"'d, 



ijk 



1 

ri2 



(3.48) 



(ii) The so-called "Leibniz" terms, which are additional contributions due to the extended 
distributional derivative, taking into account the violation of the Leibniz rule when 
performing some simplifications of the non-linear potentials at the 3PN order (see the 
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Notice that the distributional term differs from the extended-Hadamard prescription 
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explanations in Section III B in j22]). The Leibniz terms were written in Eq. (6.19) in 
|22l |. and read 

^ «i - - - g d.^u [j- j . (3.49) 

We emphasize that the contributions ()3.48p and ()3.49|1 represent some additive effects 
of the use of the distributional derivative introduced in Ref . [23] , when compared to the 
effect of the Schwartz derivative in the pHS regularization. Note that both Eqs. ()3.48|1 
and ()3.49|1 depend on the choice of distributional derivative, and we have given them 
here in the case of the "particular" derivative defined by ()3.38|) . 

(iii) A special non-distributivity in the compact-support potential V when it is computed 
at the 3PN order. In this case the Pf{F5i) comes siraply from the compact-support 
point-particle source a of the potential. Eq. (4.17) of Q gives for that term 



'^^1 = 5 — —5 — "-12- (3.50) 



(iv) A contribution coming from the compact-support part of the potential T parametrizing 
the metric at the 3PN order, and derived at the end of Section IV A in IIJ" 

f J _ 1 G^ml^jkr, f ^ \ 



5^al = -^<a,-.^-J. (3.51) 

Besides the non-distributivity of the type ()3.35|1 . we have also the more "direct" non- 
distributivity due to the fact that the pure Hadamard prescription for the regularization 
of the value of an expression "at yi", Eq. p.lip , differs from the ordinary and/or extended 



Hadamard ones [see for instance Eq. ()3.15p ]. It plays a role only in the last stage of the 
computation of the 3PN equations of motion, once we substitute all the potentials computed 
at the right FN order into the geodesic equations. We thus have 

(v) A "direct" non-distributivity contribution, which can be called non-distributivity in 
the equations of motion (EOM), and given by Eq. (6.34) in 2^ |. 
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where v''^ = fi2'^^i2 and = v{ —v{. This term involves some combinations of masses 
different from those in Eqs. (jTIH|) - (jTKT]l . Note that because {FG)i ^ (F)i(G)i the 
non-distributivity in the EOM depends on which prescription has been chosen for 
the stress-energy tensor of point-particles. Eq. ()3.52|) corresponds to the particular 
prescription advocated in Section V of ji^]. However it was checked in ji^l that 
different prescriptions yield physically equivalent equations of motion. 



^"^ Another derivative was introduced and discussed in Q where it is called the "correct" one, but it yields 
physically equivalent 3PN equations of motion. 
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The next correction brought about by the extended-Hadamard regularization is the one 
due to the regularization [F]i, performed in the Lorentzian rest frame of the particle. In 
practice the effect of such "Lorentzian" regularization boils down to applying Eq. ()3.44p . 
It turned out that the only new contribution of this type came from the regularization of 
the potential X at the IPN order [and also when deriving the result for 6^'~'^a\, which is 
Galilean- invariant, in Eq. (|3.52j) ]. leading to 

(vi) The so-called "Lorentz" contribution to the acceleration, given by Eq. (5.35) in 2^ 



as 



C6 



Q •, 1 • ' 
70 ^ 5^2 



d,,k ( ^ ) . (3.53) 



This term has been crucial for ensuring the Lorentz invariance of the final 3PN equa- 
tions of motion in 21. 2^ . 



Finally, one must also take care of one additional contribution (with respect to the "pHS" 
definitions) due to the non-Schwartzian way of treating distributional derivatives. We have 
already mentioned two contributions coming from this origin: (i) and (ii) above. Actually, 
there is a third one with the same origin and which comes from our computation (see Section 
IIVI below) of the "difference" between the dimensional and Hadamard regularizations of 
retarded potentials, namely the crucial potentials X and Ri which must both be expanded 
to IPN fractional accuracy. More precisely, this contribution is due to the repeated time- 
derivative operator coming when expanding the time-symmetric Green function of the 
d'Alembertian as D"^ = A'^ ^c'^A-'^d^ + 0{c-'^). We shall explicitly exhibit in Eqs. 
below the way these derivatives enter our calculation of the difference. For technical reasons 
the time-derivative must be kept inside the integrals, so it has to be considered in 
a distributional sense, and we have therefore to take into account the different ways of 
treating the distributional derivatives in both regularizations. In the extended-Hadamard 
regularization the distributional terms are given by Dti[F] which is shown in Eqs. (I3.41|) - 
(I3.43p . and when they enter the source of some Poisson-type integral they are evaluated 
according to Eqs. ()3.34j) . On the other hand, in dim. reg. one uses the ordinary Schwartz 
derivative (in d dimensions) which is described in Section IIVGI In this case the double 
time-derivative is computed with the help of the Gel'fand-Shilov formulae ()4.34j) - ()4.35|) 
below. When examining the difference between the contact terms in Dtt[-^] and those issued 
from liuid^s F^'''\ we find that only the source for the IPN potential X (or rather for the 
combination 4diX / which enters the equations of motion) contributes. This gives: 

(vii) The following "time- derivative" contribution to the acceleration, 

etime-derivative i _ ^ GSu^m^ , 4 G^m| ^ f) ( ^ \ ('\^A\ 



This term was part of the final result of [22|. However it is not mentioned in [2^ 
because this reference never tried to compare the results of the extended distributional 
derivative with those given by the ordinary Schwartz derivative in 3 dimensions, except 
in those cases, items (i) and (ii) above, for which the Schwartz derivative yielded in fact 
some ill-defined (formally infinite) expressions in 3 dimensions. [The latter expressions 
turn out to be rigorously zero when computed in dimensional regularization.] 
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In summary, there are in all seven different terms, (i)-(vii), which are specifically due to the 
extended version of the Hadamard regularization. The "pure-Hadamard-Schwartz" equa- 
tions of motion are then obtained from the end result of [22], i.e., a^^ given by Eq. (7.16) 
of ji^], by subtracting these terms. Therefore we define (see also Section IV Bl below) 

^pHS ^ ^BF _ j^^selfj^^ ^ ^Lcibniz^^ ^ ^ ^ d^^^S^i + ^^"''"^"ai + ^ti^e-derivative j ^ 

(3.55) 

and idem with 1 2 for the other particle. 



IV. DIMENSIONAL VERSUS HADAMARD REGULARIZATIONS 

In this Section we come to the core of our technique for evaluating the difference between 
the (i-dimensional equations of motion and their pure-Hadamard-Schwartz expressions, de- 
fined above and given in practice by Eq. ()3.55p . 



A. Iteration of Einstein's equations in d dimensions 

Let us start by indicating how we solved (with sufficient accuracy) Einstein's field equa- 
tions in d dimensions. One writes the post-Minkowskian expansion of Einstein's equations in 
the guise of explicit formulae for the elementary potentials V, V^, ■ ■ ■ , T, as given in SectionlTTl 
Note that it is crucial to take into account the explicit (i-dependence of the coefficients en- 
tering these equations. The first step of the formalism is to get sufficiently accurate explicit 
expressions for the basic linear potentials V and Vi. As we do not need to consider here 
radiation reaction effects (which do not mix with the UV divergencies arising at the 3PN 
level) it is enough to solve Eqs. ()2.7|) by means of the PN expansion of the time-symmetric 
Green function. For instance, we have 

V = -AnG n-i^cr = -4vrG (^A" V + 1 A'^ + 1 A'^ df a + A"^ df + O (1) . 

From Eq. (j2.16p we see that the source a reads 

a(x,t) = /ii(t) 5(")[x - yi(t)] + 1^2, (4.2) 

where 

/ , ""^ (■i-2) + v;/c' 

Note the presence of many "contact" evaluations of field quantities in a. Such terms are 
unambiguously defined in dimensional regularization. They are computed by successive 
iterations {e.g. to get /ii(t) to IPN fractional accuracy we need to have already computed 
g^i, to order 0(c~^) included). Those evaluations do not give rise to pole terms in a, up to 
the 3PN accuracy. Hence, as we said above, we can consider that their d — > 3 limits define 
a certain (three-dimensional) way of estimating contact terms, that we have checked to be 
in full agreement with the "pure-Hadamard" prescription defined in the previous Section. 
Coming now to the spatial dependence of the scalar potential V we get from Eq. ()4.1|) 



y(x, t) = GMt) ui + ^ [/ii(t) + ■ ■ ■ + 1 ^ 2 , (4.4) 



31 



where we introduced the following elementary solutions ui = A ^ (^—4tt6[^^^ , f i = A ^ ui, 
etc., whose explicit forms are 

Ml = krl''^, (4.5a) 
where k is related to the usual Eulerian F-function by 

.._r(¥) 



d-2 
TT 2 



(4.6) 



Inserting the explicit expression ()4.4j) of V into, say, the non-linear terms in the R.H.S. of 
Eq. ()2.12b|) . yields a d'Alembert equation for the non-linear potential Wij with a "source 
function" which is the sum of some contact terms S'(x) ^'•'^•'(x — y^) and of an extended 
non-linear source F^'^^'k) which belongs to the c?-dimensional analogue of the class J-", say 
JF^"^). More precisely, at each stage of the iteration we find inhomogeneous wave equations 
of the type 

□ iy('^)(x) = F('^)(x) + J2 ^a(x) S^'\^ - , (4.7) 

a 

where the extended source function F^^\ic) is regular everywhere except at the points yi 
and y2, in the vicinity of which it admits an expansion of the general form (VA^ G N) 

F^'\^)= E ^rVji(ni)+o(rf), (4.8) 

Po<p<N ^ 
qo<q<qi 

where p and q are relative integers (p, q & Z), whose values are limited by some po, go and qi 
as indicated. The expansion ()4.8|) differs from the corresponding expansion in 3 dimensions, 
as given in (|3.1|) . by the appearance of integer powers of where e = d — 3. The coefficients 
i/pfg depend on the unit vector rii in d dimensions, on the positions and coordinate velocities 
of the particles, and also on the characteristic length scale £o of dimensional regularization. 
Because — > F when ci — > 3 we necessarily have the constraint (Vp > Pq) 

E f?l = f^- (4-9) 
— 1 1 

i30<'?<ijl 

The iteration continues by inverting the wave operator by means of the time-symmetric 
expansion 1)4.11) . The basic terms of this expansion which will turn out to be crucial for our 
3PN calculation based on the difference are in fact the first two terms. Focussing on the 
terms generated by the extended source F*^'^-' (x) (rather than the simpler contact terms) we 
can write the d-dimensional analogue of ()3.25|) as 

i?W(x') = □-Lt^^'Hx)] = P^'K^) + ^,Q^'\^) + ^ (^) ' (4-10) 



The constant k adopted here is related hy k — ink to the constant k chosen in [35j. Our present choice 
is motivated by the easy-to-remember fact that hm^^s fc = 1. 
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where the d-dimensional Poisson integral of F^'^'^ reads 



pCOlx') = A- [Fl*(x)] ^ -i- / ^-^^F''H^) . (4.11) 

We have used the fact, already mentioned above, that the d-dimensional elementary solution 
of the Laplacian reads 

A (k\x - xf = -Att5^'^\x - x') , (4.12) 

(see Appendix El for a proof of Eq. ()4.12j) and for other useful formulae valid in d dim.), 
while the IPN term is given by 

Q'^'^\^') = 2A-'[d^F'^'^\^)]=- ^^^^_^J d'^^\^-^r-'d^F'-'\^). (4.13) 

Note the important point that in d dimensions, as in 3 dimensions, the time-derivative 
operator df present in the integrand of ()4.13j) is to be considered in the sense of distributions 
(see further discussion in Section HV CI below) . 

An important technical aspect of the (i-dimensional FN iteration of the elementary po- 
tentials V, ■ ■ ■ , T is the existence of the generalization ()4.12j) of the usual Green's function 
for the Laplace equation, as well as of its higher PN analogues A-"5(^), allowing one to 
explicitly compute the spatial dependence of the linear potentials V, Vi and K for in- 
stance. However, starting with Wij we need to Poisson-integrate non-linear sources, such 
as A~^{di U dj U). In three dimensions, these non-linear contributions are reducible to the 
knowledge of the basic non-linear potential g, such that Ag = r^^ . We have succeeded 
in explicitly computing the (i-dimensional analogue of the g potential, namely 

g^''\^) ^ A-\rl-' rl-') . (4.14) 

Our result is reported in Appendix O As indicated there, if we wished to explicitly compute 
some of the higher PN potentials needed to write the closed form of the non-linear sources 
relevant to the 3PN equations of motion, we should extend the calculation of the potential 
^'^^ to the potentials /'^'^^ and of Appendix O 

Luckily, it is not needed to use a closed-form expression for any of the non-linear poten- 
tials. Indeed, similarly to what was used long ago [8] when discussing the iteration gen- 
erated by Riesz-type sources, Eq. (jl.3|) . one can control the UV singular part of □~^F(x) 
from the knowledge of the UV singular part of its non-linear source F(x).^^ More pre- 
cisely, in the vicinity say of yi, at each iteration stage we can decompose the source in 
F(x) = Singp,(x) + Reg^(x) where the singular part Sing^(x) (with respect to yi) is a sum 
of terms of the form Eq. ()4.8|) . which are not (in the limit d — > 3) smooth functions of x — yi, 
and where the regular part Regp^(x) is a smooth (C°°) function of x — yi. [The simplest ex- 
ample of this decomposition is Eq. ()2.21|) with, near point yi, Sing[;(x) = fkGmirl~'^ 
and Reg^(x) = [/2(x).] If, for concreteness, we then consider P(x) = A"^F(x), the 
above decomposition entails a corresponding decomposition of -P(x), and it is easy to see 



We have checked that for aU the non-compact (extended) potentials involved in this calculation, there are 
no IR divergencies, i.e., the integrals converge at infinity |x| — > oo for any small enough value of e = d — 3. 
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that ASingp(x) = Singj7(x). From this result, we can uniquely determine Singp(x) from 
Singp(x) using, e.g., the formula ()B26c|) in Appendix iBl This local procedure does not 
allow one to compute the regular part of the Poisson potential in d dimensions. Fortunately, 
thanks to particular simplifications that occur in the structure of Einstein's field equations, 
the knowledge of Regp(x) in d dimensions for the complicated non-linear sources is not 
needed. Indeed, one can see on our explicit solution of Einstein's field equations at 3PN 
given in Section IH] that there are no "quartically non-linear" source terms of the form, say, 
diVdjVWij or diWjkdjWki for goo at the 3PN order (see Fig. El below). 

As explained in Section IV CI below, a nice way to understand the origin of the poles oc 
((i — 3)^^ appearing in the 3PN equations of motion is to use a diagrammatic representation. 
A pole can arise in ai only when three propagator lines (including the extra one coming 
from when solving □(7^,^ = non-linear source) can all shrink towards the first world- 
line. If terms of the type above {e.g. diVdjVWij) were present in the source one could 
have a diagram where the three shrinking propagators come from D"^, diVi and djVi. Then 
Reg^'^) [W^y(x)] would remain as an external attachment to this diagram (and would then 
fork into two "feet" on the second world-line). In view of the pole cx (with e = d — 3) 
arising from the triplet of shrinking propagators, one would need to know Reg^'^-*[iyij(x)] up 
to e accuracy, i.e., Reg^'^\Wij{x)] = Reg^^\W^j{x)] +£W^^(x) + 0{e^) [in which Wlj{x) is 
defined by this expansion] . If such a term had been present we would have needed to use the 
full d-dimensional, globally determined g-potential given in Appendix O to determine W^j, 
which would have entered the final, renormalized equations of motion. However, because 
all such terms are absent at the 3PN order, the only external attachments to the dangerous 
shrinking diagrams are simple lines, such for instance as the lines ending on y2{t) in Figs. Eli, 
ISb orHjo presented below. Such lines do need to be evaluated to accuracy e, but this is easy 
because they represent linear potentials such as V or Vi which are known in dimension d via 
Eq. (lOD . 

In conclusion, the algorithm we use to solve, with sufficient accuracy, Einstein's equations 
in d dimensions consists of: (1) starting from the fully (i-dimensional expressions for the 
linear potentials V, Vi (and more generally for the parts of the non-linear potentials with 
delta-function sources); (2) determining the local expansions, near ya, of the singular parts 
of the non-linear potentials by inverting ASingp''(x) = Singp^(x) via formulae ()B26c|) of 
Appendix iB) (3) completing P^'^^x.) by adding to Singp^(x) the limit when ci — > 3 of 
Regp^(x), namely Regp'*(x) which is known from the previous work on the 3PN equations 

of motion in 3-dimensions |22]. Note that we denote by Regp (x) a formal d-dimensional 
function, x G M'^, the explicit expression of which in terms of ri, ni, etc. coincides with 
its 3-dimensional counterpart. For instance Reg^^^(x) denotes the usual regular part of 
^'^\x), obtained by subtracting from g'^^)(x) = ln(ri + r2 + ri2) the two three-dimensional 
locally singular expansions of A^^(rf^r^^) around yi and y2 as given by the d ^ 3 limit 

of Eq. ()C9|) and its 1 ^ 2 analog. After this double subtraction, Reg^^^(x) is considered 
as a function in M'^, and we can use as approximation to ^^^{x) the explicit expression 
S'loc i(^) +S'ioc 2(^) + R-sg^^-* (x) . More generally, in our calculations we use as approximation 
to P^'^^ [which symbolizes here the non-linear potentials Ri and Zij at Newtonian order, and 
Wij at the IPN order] the expression Singp''(x) + Regp''(x). Evidently, the subtraction of 
the singular part needs to be performed only up to some finite order in and r^. We have 
checked the choice we made of N in each calculation by doing two separate calculations 
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for the values and + 1, and checking that the corresponding final results are the 
same. We performed also direct checks of the independence of the final results on the 
precise (i-dimensional extensions of the "regular" part of the non-linear potentials, such as 
Regp^(x) = Regp'*(x) + eP'(x) + 0{e'^) [in which P'(x) is defined by this expansion]. We 
systematically added in all our non-linear potentials Ri, ■ ■ ■ , Wij some smooth contributions 
to Regp''(x) vanishing with e, i.e., some substitutes for the actual -P'(x). These "substitutes" 
were determined in such a way that (i) they are homogeneous solutions of the d'Alembertian 
equation at the required post-Newtonian order, (ii) the differential identities obeyed by the 
potentials in d dimensions, Eqs. fj2.13ap - ()2.13b|) . are indeed satisfied up to the order e, and 
with the required precision N in powers^'' of ri or r2. And we checked that our final results 
are totally insensitive to the introduction of such substitutes for the function -P'(x). 

Finally, when evaluating the equations of motion, as given by Eq. (j2.23|) . we must evaluate 
the value at x' = yi of many terms given either by Poisson integrals of the form (j4.11|) or 
their IPN generalizations (j4.13|) . This is quite easy to do in dim. reg., because the nice 
properties of analytic continuation allow simply to get [P'''^^(x')]x'=yi (say) by replacing x' 
by yi in the explicit integral form ()4.11|) . Finally, we simply have for the values at x' = yi 
of the potentials, 

P^'\yi) = -^l (4.15a) 
Q^'\yi) = -^^^(^r^ / d''^rt'd^F(^\^) , (4.15b) 
as well as for their spatial gradients. 

As said above, the main technical step of our strategy will then consist of computing the 
difference between such (i-dimensional Poisson- type potentials ()4.15j) or ()4.1fj|l . and their 
"pure Hadamard-Schwartz" 3-dimensional counterparts, which were already obtained in 
Section ITTTBI 



k{d-2) 1^,^ n 



4:71 
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d'^xn\rl-'^d^F^'^\^) 



(4.16a) 
(4.16b) 



B. Difference between the dimensional and Hadamard regularizations 

We denote the difference between the prescriptions of dimensional and "pure Hadamard- 
Schwartz" regularizations by means of the script letter V. Given the results (P)i and 
P'^'^^(yi) of the two regularizations [respectively obtained in Eqs. ()3.2H) and ()4.15a|) ] we pose 

PP(l) = p(^)(yi)-(P)i. (4.17) 



Therefore, our verification that the potentials we need do satisfy the harmonicity conditions (|2.13|l has 
been done only in the vicinity of the two particles. 
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That is, 'DP(l) is what we shall have to add to the pure Hadamard-Schwartz result ()3.55|) 
in order to get the correct (i- dimensional result. Note that, in this paper, we shall only 
compute the first two terms, a_i + ao + 0{e), of the Laurent expansion of VP{1) when 
e — i> 0. This is the information we shall need to fix the value of the parameter A. We leave to 
future work an eventual computation of the d-dimensional equations of motion as an exact 
function of the complex number d. 

Similarly to the evaluation of the difference VH = H'^'^^ — Hadamard[if ^^''] in Ref. jsij ]. 
the difference ()4.17|) can be obtained by splitting the dimensional integral ()4.15a|) into 
three volumes, two spherical balls Bf\s) and b'^\s) of radius s and centered on the two 
singularities, and the external volume \ Bf\s) U B2^\s). When d 3 (with fixed s), 
b[^\s) and B2^\s) tend to the regularization volumes -Bi(s) and B2{s) we introduced in 
Eq. (|3.22|) . Consider first, for a given value s > 0, the external integral, over \ b['^\s) U 
B2^\s). [If wished, two balls with different radii could be used, with the same result.] Since 
the integrand is regular on this domain, it is clear that the external integral reduces in the 
limit e — s> to the one in 3 dimensions that is part of the Hadamard regularization ()3.22|) . 
So we can write (for any s > 0) 
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\B[''\s)uB!f\s) r 



'^F''\-) = -i[ ^F(x) + (9(.), (4.18) 

1 47r Jm.3\Bi(s)UB2(s) ^1 



and we see that when computing the difference VP{1) the exterior contributions will cancel 
out modulo 0{e). Thus we obtain, after this prehminary step [following Eq. (j3.22|) ]. 
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(4.19) 



See Section IV of Ref. ^35] for a careful justification of the formal interversions of limits 
s ^ and e — > that we shall do here. The point is that in order to obtain the difference 
T)P{1) we do not need the expression of F^'^^ for an arbitrary source point x G M'^ but only in 



the vicinity of the two singularities: indeed the two local integrals over B^f{s) and B^^'is) 
in Eq. ()4.19|1 can be computed by replacing F^'^^ by its expansions when ri — and r2 ^ 
respectively. We substitute the ri-expansion Eq. ()4.8j) into the local integral over b[^\s), 
and integrate that expansion term by term. This readily leads to 
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where we still use the bracket notation to denote the angular average, but now performed 
in d dimensions, i.e., 

dVld^iini] 
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d-l 



(4.21) 
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Here dQ,d-iis the solid angle element around the direction rii, and 1 = 2712 /T (I) is the 
volume of the unit sphere with d — 1 dimensions (see Appendix El for more discussion). To 
derive (j4.2Up we used the following relation linking k and ^Id-i, 

*^(^^- 

Concerning the other local integral, over things little bit more involved because 

we need to perform a multipolar re-expansion of the factor r^""^ present in that integral 
around the point y2. Writing down this multipole expansion presents no problem, and in 
symmetric-trace-free (STF) form it reads 

-r' = f,^-^dJ^)4u^. (4.23) 

The multipole expansion being then correctly taken into account, we obtain 

(4.24) 

As we can see, simple poles ~ 1/e will occur into our two local integrals, as determined by 
(I4.2(J|) and (j4.24p . only for the "critical" values p = —2 and p = —i — 3 respectively. 

Next we replace the explicit expressions ()4.20|) and ()4.24j) into the formula ()4.19|) we had 
for the "difference". As expected we find that the divergencies when s — 0, some value 
6^0 being given, cancel out between Eqs. ()4.20|) - ()4.24|) and the remaining terms in ()4.19p . 
so that the result is finite for any e ^ 0. Furthermore, we find that if we neglect terms of 
order 0{e), the only contributions which remain are the ones coming from the poles (and 
their associated finite part), i.e., for the latter critical values p = —2 in the case of singularity 
1 and p = — £ — 3 in the case of singularity 2. The other contributions in (j4.2U|) and (|4.24p 
have a finite limit when e ^ which is therefore cancelled by the corresponding terms in 
Hadamard's regularization. As a result we obtain the following closed-form expression for 
the difference, which will constitute the basis of all the practical calculations of the present 
paper, 

+ 0{6). (4.25) 

Notice that ()4.25|) depends on the two "constants" Inr'j^ and lns2. As we shall check these 
Inr'j^ and lns2 will exactly cancel out the same constants present in the "pure-Hadamard" 



The expansion is STF because Ar^ = in d dimensions (in the sense of functions). See Appendix IbI 
for a compendium of d-dimensional formulae on STF expansions. See also Eq. (|C6|I in Appendix [CJ 
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calculation, so that the dimensionally regularized acceleration will be finally free of the 
constants r[ and S2. Note also that the coefficients i/pfg and in d dimensions depend 
on the length scale £q associated with dimensional regular ization [see Eq. ()2.4|) ]. Taking 
this dependence into account one can verify that r[ and S2 in (I4.25p appear only in the 
combinations ln(r^/£o) and ln(s2/^o)- 

Let us give also (without proof) the formula for the difference between the gradients of 
potentials, i.e., 

Vd,P{l) = diP^''\yi) - (a,P)i . (4.26) 
The formula is readily obtained by the same method as before, and we have 



^ <io<q<qi 



L As) 



+ 0{e) . (4.27) 

Formulae ()4.25|) and ()4.27|) correspond to the difference of Poisson integrals. But we have 
already discussed that we shall need also the difference of inverse d'Alembertian integrals 
at the IPN order. To express as simply as possible the IPN-accurate generalizations of 
Eqs. (I4.25P and (j4.27j) . let us define two functionals H and Hi which are such that their 
actions on any ci-dimensional function F^'^^ is given by the R.H.S.'s of Eqs. ()4.25p and ()4.27p . 
i.e., so that 

VP{1) = n [F^'^^ , (4.28a) 
VdiP{l) = Hi [F^'^)] . (4.28b) 
The difference of IPN-retarded potentials and gradients of potentials is denoted 

VR{1) = i?('^)(yi)-(P)i, (4.29a) 
VdiR{l) = 9,i?("Hyi) - (diR), , (4.29b) 

where in 3 dimensions the potential -R(x') is defined by Eq. ()3.25|) and the regularized values 
(i?)i and {diR)i follow f rom (IT^, (HT^, (ITTfll . and where in d dimensions R^'^^yi) and 
diR^'^\yi) are given by ()4.10j) . ()4.15|) . ()4.16p . With this notation we now have our result, 
which will be stated without proof, that the difference in the case of such IPN-expanded 
potentials reads in terms of the above defined functionals 7i and Tij as 

VR{1) = 



Vd^RiX) = 



p{d) , n a2 p{d) 
^ ^2c\A-dr^ 
2 

p{d) _ ^1 f)2 p{d) 

2c\d-2r^ 



-^(f-4> + 0(^), (4-30a) 
-^(ni^„3> + o(^). (4.30b) 



These formulae involve some "effective" functions which are to be inserted into the functional 
brackets of Ti and Tii. Beware of the fact that the effective functions are not the same in the 
cases of a potential and the gradient of that potential. Note the presence, besides the main 
terms 7i[- ■ ■] and 7ii[- ■ ■], of some extra terms, purely of order IPN, in Eqs. ()4.3()j) . These 
terms are made of the average of some coefficients ikp of the powers r\ in the expansion 
when ri — >■ of the second-time-derivative of F, namely d^F . They do not seem to admit 
a simple interpretation. They are important to get the final correct result. 
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C. Distributional derivatives in d dimensions 



Let us end this Section by explaining in more details how we dealt with distributional 
derivatives in d dimensions. First, it is clear that if we were dealing with d-dimensional 
integrals of the type 

1= [ d''xipij{x)dijUi, (4.31) 



where (piji'x.) is some (formally) everywhere smooth function of x e M"^, with fast enough de- 
cay at infinity, and where Ui = (— 47r(5|'^'') is the elementary Newtonian potential in d di- 
mensions [see Eqs. ()4.5ap above], we should, in a straightforward ci-continuation of Schwartz 
distributional derivatives, consider that dijUi contains, besides an "ordinary" singular func- 
tion dij{ui)\^^^ (treated as a pseudo-function in the sense of Schwartz), a distributional part 
proportional to S^'^\x — yi). In other words, we would write 

Ayr 

d,,{ui) = d,,{u,)\ _ 5^^.5W(x-yi), (4.32a) 

I ord " 

dijki{vi) = dijkiivi)^ - ^^^^ (SijSki + 5ik5ji + 5iiSjk^5^'^\x - yi) , (4.32b) 

where the indication "ord" refers to the "ordinary" (pseudo-function) part of the repeated 
derivative. We have also added the corresponding result for the fourth derivatives of the 
"less singular" kernel Vi = ui, Eqs. fl4.5b|) . Note that the decompositions above of 
dij Ui or dijki f 1 into "ordinary" and "distributional" pieces arise because of our working in 
((i-dimensional) x-space, and of explicitly computing some derivatives, say as dij (^^1)10,^ — 
[n 6ij+n{n—2) n\ n\] r"~^. If we were working in the (d-dimensional) Fourier-transform space 
k (which is where dimensional continuation is most clearly defined jl^), the corresponding 
decomposition would be simply algebraic: e.g. /k.^ = A;^*-^^/k^ -|- d~^Sij, where fc^*-^^ 
denotes the STF part of k^^ = k^kK 

The decompositions ()4.32|) are clearly needed when dealing with simple integrals of the 
type ()4.3H) (with a smooth 93 (x)) to ensure consistency with the requirement that one 
may integrate by parts (which is one of the defining properties of dim. reg. [l^l), and 
we shall therefore employ it, when applicable. On the other hand, most of the singular 
integrals that we have to deal with look like (j4.3H) but contain a singular function (/^(x), 
of the type of Eq. (j4.8p . It is, however, a very simplifying feature of dim. reg. that 
when considering integrals like 1)4.311) with some singular V5(x) we can simply ignore any 
distributional contributions oc 5'^'^)(x — yi) or its derivatives. Indeed, as long as the integer 
q in the powers rf"*"^^ present in ()4.8|) is different from zero (which is precisely the case of 
all delicate terms involving several propagators shrinking towards a particle world-line), the 
"singular" expansion ()4.8|) can be considered, in dim. reg., as defining a sufficiently smooth 
function [by taking both qe and large enough in ()4.8p ] which vanishes, as well as its 
derivatives, at x = yi. Therefore, all the "dangerous" terms of the form Sing^''(x)(5'''^''(x— yi) 
unambiguously vanish in dim. reg. 

Let us now consider the consequences of this fact for the time derivatives occurring in 
expansions such as Eqs. ()4.3(jp . The distributional time-derivatives, acting in our present 
example on ui or fi, i.e., on functions of r\ = x"^ — y[{t), can be treated in a simple way 
from the rule dt = —v\di applicable to the purely distributional part of the derivative. For 
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instance we can write 

dUui) = dKui)\ --v?5W(x-yi), (4.33a) 

I ord 

d^d,,{v,) = d^d,,{v,)^ ^ - + - yi) • (4-33b) 

We have checked using these formulae that all the ci- dimensional terms coming from second- 
order derivatives of potentials, taken in the distributional sense (for instance the term 
WijdijV in the source of the X-potential yield the same purely distributional contri- 
butions, in the limit e — > 0, as the ones that would be computed using what we called above 
a "pure Schwartz", three-dimensional computation of such contributions [to "smooth" in- 
tegrals ()4.31| ]. On the other hand, the extended version of distributional derivatives intro- 
duced in j23| does yield some specific additional contributions, two of which were already 
mentioned in [i^l and are reported in the items (i) and (ii) of Section 1111 PI above, and a 
third one (also included in j2^) which comes in connection with the second time-derivatives 
in our formulae for the difference, Eqs. (j4.3U|) . 

Let us indicate here that the distributional second-time-derivatives in d dimensions have 



been obtained by using the following (generalizations of) Gel'fand-Shilov formulae |4J|, valid 
for general functions F^'^\'x) admitting some expansions of the type ()4.8|) : namely, for the 
spatial derivative, 

d,F^'^ = I + J2 f 3-2,-1)5^5^ + 1 - 2 , (4.34) 

lord ^- 1 

where diS^^ is the £-th partial derivative of the d- dimensional Dirac delta-function at the 
point 1 (L = iii2 ■ ■ - ig) and where the angular average is performed over the {d — 1)- 
dimensional sphere having total volume and, concerning the time derivative, 

dtF^'^ = dtF^'^ I - n^-i J2 -J-i^ii^iVi) f ^-1-2,-1)9 A'^ + 1 - 2 . (4.35) 

lord «• 1 

From the latter formula one deduces the second time- derivative in a way similar to 
Eqs. fl3.43|l . We have indicated in the item (vii) of Section Fill Dl the correction it leads to 
when comparing with the extended-Hadamard prescription for the second time-derivative, 
and we have subtracted it from a^^ to define the pure Hadamard-Schwartz result ()3.55p . 
Therefore, we consistently do not need to include such an effect into the differences VP{1) 
discussed here. 

Finally we are now in position to obtain the supplement of acceleration Dai induced by 
dimensional regularization, which is composed of the sum of all the differences of potentials 
and their gradients computed by means of the generals formulae of (j4.25p . (j4.27|) and (|4.3(J|) . 
The term Dai when added to the "pure-Hadamard-Schwartz" acceleration defined by (|3.55|) . 
gives our result for the dimensionally regularized ( "dr" ) acceleration 

af = af^ + Dai and 1^2. (4.36) 

More details on the practical computation of Da^ (which parts of the potentials contribute; 
what is the diagrammatic picture) will be given in Section FV CI 



Since this term is to be computed at the IPN order, not only does it contain second-order derivatives of 
ui , but also fourth-order derivatives acting on vi . 
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V. DIMENSIONAL REGULARIZATION OF THE EQUATIONS OF MOTION 



A. Structure of the dimensionally regularized equations of motion 

The preceding Section has explained the method we used to compute the dimensionally 
regularized equations of motion as the sum (a = 1,2; considered modulo 2) 

s-t'lejo] = Sif'^iCsa+ij+Vsi^lCsa+ueJo] , (5.1) 

where the label "pHS" refers to the "pure Hadamard-Schwartz" definition of the accelera- 
tion {i.e., the "raw" result of ji^l, after subtraction of the additional contributions quoted 
in Section Qll PI above. Eq. (j3.55|) ). and where Pa^ is the difference induced when using di- 
mensional continuation as regularization method, instead of Hadamard's one. A first check 
on our results will be that, as indicated in (j5.1|) . the four regularization parameters (with 
dimension of length), r'^, rg, Si, S2, that enter the Hadamard method must cancel between 
aP^^ and Pa^ to leave a result for the dimensionally regularized accelerations a^'' which 
depends only on the two regularization parameters of dimensional continuation: e = d — 3 
and the basic length scale io entering Newton's constant in d dimensions, G = Gn^q-, where 
we recall that Gn denotes the usual three-dimensional Newton constant. 
The dimensionally regularized acceleration 1)5.11) has the structure 

af [yi2, Vi, V2] = aNa[yi2] + aiPNa[yi2, Vi, V2] 

+ a2PNa[yi2, Vi, V2] + a2.5PNa[yi2, Vi, V2] + a3PNa[yi2, Vi, V2] , (5.2) 

where we denote yi2 = yi — y2- The 3PN term (which is the only one to have a pole at 
£ = 0) has a tensor structure of the form (say for the first particle, a = 1) 

Ani2 + 5'vi -5"v2, (5.3) 

where, as usual, ni2 = yi2/''"i2 denotes the unit vector directed from particle 2 to particle 1. 
The scalar coefficients A, B', B" entering the equation of motion of yi can be decomposed 
in powers of the masses, say 

A = 2^ c„,„2(vi,V2,ni2,lnri2) ^ m+n^+i ' (5-4a) 

l<ni+n2<4 '"is 

B = 2^ c„^„2(vi,V2,ni2,lnri2) ^ m+n^+i ' (5-4b) 

l<ni+n2<3 ^ ^12 

= Z2 c;;in2(vi,V2,ni2,lnri2) ^ m+^+i ^ ' (^-4^) 



l<7ii+n2<3 



c" r 



where ni and n2 are natural integers, with the restrictions indicated. Note that, in Eqs. ()5.4p . 
we have conventionally factored out an integer power of the "full" (rf-dimensional) gravi- 
tational constant G, and a corresponding integer power of ri2. This creates a mismatch 
between the usual 3-dimensional dimension of, say, ch^f and the dimension of c^^rij . Using 
G = Gm^o one sees that it is the combination ^o"^'''"^''^c„j„2 which has the same dimension 
as CnrnJ- Alternatively said, the ensuing fact that r[2^'^'^^''^c„^„2 has the same dimension as 
Cnrnf implies, as indicated in Eqs. ()5.4p . a dependence of c„j„2 lnri2 when e — * 0. Notice 
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also that in Eq. ()5.3p we have introduced separate notations for the coefficient of vi and that 
of V2. Actually, the Poincare invariance of the equations of motion imposes the restriction 
B' = B" so that the last two terms in Eq. ()5.3|) are proportional to the relative velocity 
Vi2 = Vi — V2. [Note, however, that B' is not a function of Vi2 only; it depends both on Vi 
and V2.] Because the calculation of the separate contributions aP^^ and Pa^ to the equations 
of notion breaks the over-all Poincare invariance of the formalism, our computation of the 
separate pieces a^^^ and Pa^ will involve partial contributions to B' and B" that do not 
coincide. It is only at the end of the calculation that the equality B = B' = B" will be 
satisfied, so that finally 

asPNi = Ani2 + B Vi2 . (5.5) 

Most of the coefficients c„^„2 5 c'^^^^; ^nin2 entering the 3PN acceleration are well behaved 
when e —>■ 0, in the sense that their evaluation never involves any poles oc By this we 
mean that whatever be the (reasonable) way of decomposing the integral giving a coefficient 
in separate contributions, the latter contributions do not involve poles oc 1/e. The subset of 
coefficients whose evaluation involves poles coincides with the set of "delicate" coefficients 
in the Hadamard regularization, namely the nine coefficients contributing to terms of the 
following form in the acceleration of the ffist particle: 



Q4 



.6 /^iS 



12 



C31 ml m2 + C22 m\ m\ + C13 mi ml 



m\ m2 



,6 ».4 



12 



ni2 



C2l(Vi, V2)ni2 + C2i(Vi, V2)Vi " Cgi ( Vi , V2) V2 



+ 



ml 



12 



C03(vi, V2) ni2 + Co3(vi, V2) Vi - Co3(vi, V2) V2 



(5.6) 



The ffist three terms in Eq. (j5.6j) do not depend on velocities and will be referred to as the 
static delicate contributions, by contrast to the kinetic delicate contributions involving the 
velocity-dependent coefficients C21, Cgi, Cg^, C03, C03, and Cq3 (they depend on vi, V2 and also 
on 1112). 



B. Pure-Hadamard-Schwartz static contributions to the equations of motion 

Correspondingly to the decomposition ()5.ip of the equations of motion, the dimensionally 
regularized static contributions^^ cf^, C22, 0^3 to the acceleration of the ffist particle can be 
written as the sum (m + n = 4, m > 1, n > 1) 

ctn[e] = ^2] + Vc^n[r[, ^2, s] . (5.7) 

In this subsection, we discuss the explicit evaluation of the pure Hadamard-Schwartz static 
coefficients c^^. 

As explained in the previous Section, the pHS static contributions c^^^fr'^^, S2] are obtained 
from the results reported in ^23\ by undoing two things. First, the "BP" results reported 
there for ai (Eq. (7.16) of [22j) were expressed in terms of the three parameters r[, r'^ and 



As explained above, we consider only the "delicate" ones. In the present case, this means that we do not 
consider the contribution = 16 + 0{e), unambiguously obtained from the test-mass limit m\ ^ m^- 
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A, instead of the two pure Hadamard parameters r[ and S2 more relevant for the present 
purpose. The introduction of the parameter A was motivated by requiring that the full set of 
equations of motion (which a priori depended on four independent regularizing parameters 
r[, Si, S2) admit a conserved energy. This led to the link, Eqs. (7.9) in jj^. 

159 , nil + ^2 



In = _ + A ^ " . (5.8) 

\S2j 308 7712 ^ ^ 

When inserting ()5.8j) in the expression of af^[r']^, A] we find, as it should be, that the 
result simplifies to an expression depending only on the two pure Hadamard parameters r[ 
and S2- This leads to the following net results from [j^ . 

RP, , , 34763 41 2 / , X 

^ 2 gg^ 

22 L 1. 2J 210 16 ' ^ ^ 
BFr , 1 1565 41 2 44, fr^ 



Second, Ref. [l^l obtained their results for the equations of motion by adding to the 
pure Hadamard-Schwartz contributions 7 additional corrections, imposed by their extended- 
Hadamard regularization and explained in Section Qll Dl above: see the items (i)-(vii) there. 
Note that these various corrections affect the "delicate" contributions to ai, in general both 
the static and kinetic ones, but only five of them contribute to the static part. These are the 
self term ()3.48p . the Leibniz term ()3.49p . the ^-correction given by ()3.50p . the EOM non- 
distributivity ()3.52|) . and the distributional time- derivative one ()3.54j) . Following Eq. ()3.55p . 
and focussing on the static contributions, we now subtract these static terms from the result 
()5.9|) in order to get the looked-for pure Hadamard contributions: 

4rKs2] = cffK,.^]-^, (5.10a) 



210 
97 

cir[r[,S2] = c^![r[,S2] + — , (5.10b) 



i.e., exphcitly. 



4"^[rl, .2] = cfiK ^2] - 5 + f - 1^ + A , (5.10c) 



pHSr / 1 1123 44 fru\ ^ ^ . 

4rK,S2] = l66-^n\ (5.11b) 

16 



We use here the link corresponding to the "particular" improved distributional derivative Di[-F] given 
by (|3.38(l above. Another derivative, the "correct" one, was also considered in and shown to yield 
equivalent equations of motion. The pure Hadamard result does not depend on this choice because we 
shall subtract below the specific additional contributions coming from the distributional derivative [F] . 
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Note in passing that though the coefficient C22 does not contain regularization logarithms, 
its evaluation involves many intermediate logarithmic divergencies that cancel in the final 
result. Such "cancelled logs" lead to as much ambiguity in the final result than uncancelled 
ones that explicitly depend on an arbitrary regularization scale such as r[ or S2 in C31 or C13. 



C. Dimensionally regularized static contributions 

We now turn to the evaluation of the "dim.-reg. minus pure-Hadamard" differences 'Dcmn 
in Eq. ()5.7|) . coming from the differences Vsli in Eq. ()5.1|) . We start from the (i-dimensional 
expression for the acceleration ai [see Eq. (j2.23p for a short-hand form], which is itself 
expressed in terms of the (i-dimensional elementary potentials V , Vi, K, Wij, Ri, X, Ztj, 
Yi and T defined in Section El Each elementary potential can be naturally decomposed in 
a "compact" (or, equivalently, "contact") piece (whose source is compact, i.e., involves the 
basic delta-function sources a, (Tj, aij) and a "non-compact" one (whose source is non-linearly 
generated and extends all over space). The potentials V, Vi and K are purely "compact", 
V = , Vi = Vf, K = , while all the other potentials admit a decomposition of the 
form Wij = + Wj^*", etc. For instance, the "compact" part of Wij is defined by 

□ W^c- = _47rG (^a,, - , (5.12) 

while its "non-compact" part is defined by 



□W^r = -^^5.^5.^- (5.13) 

C I vNC 



A more complicated example is the potential X = X + X with 



DX^ = -47tG 

and 



1 d-3 _ fd-SV /I 



d-2 " d-1 \d-l \2 



Van + 2 — — aiVi+{ — — 1 a { ^V^ + K 



(5.14) 



□ X^^ = W,jd,,V -2d,V,d,Vi + 2VdtdiV 

The NC contribution can be further decomposed into the piece of X^*" whose source is 
quadratic in compact potentials, namely 

□ X^v ^ y^c Q,, V -2d, Vj dj Vi + other VV terms , (5.16) 
and its "cubically non-compact" piece given by 

□ XCNC ^ ^Nc y ^ 1 d_zl vd, V^ d., V . (5.17) 
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To get a feeling of the actual evaluation of the difference Va.i let us consider a specific 
contribution to ai, say the term 

a\[X]^^^ida).=y,. (5.18) 

It can be decomposed into: (i) its "compact" piece ai[X'^], (ii) its "quadratically non- 
compact" one ai[X^^], (iii) its "cubically non-compact" part ai[X'^^'-^]. 

It is sometimes convenient to think of the various contributions to ai in terms of space- 
time diagrams. If we represent the basic delta-function sources [proportional to mi 5(x — yi) 
and 777,2 (5(x — y2)] as two world-lines and each propagator D"^ as a dotted line, a "compact" 
contribution to ai will be represented by one of the diagrams in Fig. ^ For instance, Fig. 




FIG. 1: Diagrams representing "compact" contributions to acceleration ai. The dotted line repre- 
sents D"^, the cross represents the field point x (here taken on the first worldline), and the bullet 
represents either a source point or (in the Figures below) an intermediate nonlinear vertex. 

can represent a term {diV)i in a\ in which the (compact) source a of is proportional 
to 77i2'5(x — y2) and involves no further powers of the masses, while Fig. represents a 
self-action term^^ in (9j V)i with source proportional to mi 5(x — yi). By contrast. Fig. 
might correspond to another term in [di V)i where the compact source cr is concentrated at 



While in usual regularization schemes using dimensionful cut-ofFs {e.g. small length scales si, S2) the 
self-action diagrams, such as Fig. ^ or Fig. are the first divergencies that one encounters and must 
then renormalize away, dimensional regularization has the technically useful property of setting all of these 
diagrams to zero. Indeed, when using a time-symmetric propagator D"^ = + + • • • these 

diagrams are seen to involve the coinciding-point limits of |x — yip"''^^", which vanish when x ^ yi by 
dimensional continuation in d. 
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y2, o"2 = (^2 ^(x — y2), and where a part of the "effective mass", 

2 maC (c/-2)+v|/c2 



Ai2 



(5.19) 



contains, besides the overall factor m2, another factor mi. As all the sources of X'^ contain, 
besides some "basic" a^^^, a potential (V, Vi, or K), the diagrams contained in a^^fX*-^] 
will be at least of the form of Fig. or will involve a more complicated mass- 

dependence. 

The quadratically non-compact terms a\ [X^^] will then contain diagrams of the type of 
Fig. 121 while the cubically non-compact term ai[X*"^'"] contains many subdiagrams of the 
type sketched in Fig. IHl 



1 2 
a 



1 2 
b 



1 2 
c 



etc. 



1 2 
d 



FIG. 2: Quadratically non-compact contributions to acceleration ai. 



etc. 



1 2 
a 



1 2 
b 



1 2 
c 



1 2 
d 



FIG. 3: Cubically non-compact contributions to acceleration ai. 

The particular term (j5.18j) that we considered contains only diagrams of the general type of 
Figs. or 01 Note, however, that there are also more non-linear contributions to a\, such 
as some terms in 

aim = ^(9,f)i, (5.20) 

corresponding to diagrams of the type sketched in Fig. |3] Similarly to the diagrams Fig. ^ 
or Fig. |21i, all the diagrams above can be modified by the presence of additional lines propa- 
gating directly between the two world-lines and corresponding to "potential" modifications 
of compact-support sources. 

As underlined in Section flV AI above, the 3PN equations of motion do not involve "quar- 
tically non-linear" contributions corresponding to diagrams such as those of Fig. El Terms 
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1 2 
a 



1 2 
b 



1 2 
c 



1 2 
d 



FIG. 4: Other non- linear contributions to acceleration ai. 



like A ^ (V'^dkVdkV) or A ^ ydkVdkXj are of this form, and they do occur in the 3PN 

acceleration ai, but since they involve double contracted gradients, it was possible to inte- 
grate them away thanks to rule (ii) of Section InJ see Eq. ()A12|) in Appendix ^ below. On 
the other hand, terms of the form diVdjVWij or diWjkdjWki do not occur at the 3PN order 
0(l/c^), although they are of the third post- Minkowskian order 0(G'^). 



•••1 1 etc. 



1 2 

a 



1 2 
b 



1 2 
c 



1 2 
d 



FIG. 5: Quartically non-compact contributions which do not occur in our calculation of acceleration 
ai at the 3PN order. 



Drawing diagrams often helps to highlight the nature of the UV singularities contained in 
the integrals they represent. As a rule of thumb, the "delicate" diagrams, that might involve 
poles, or cancelled poles, when e — > (corresponding to logarithms, or cancelled logarithms, 
in d = 3) are characterized by the presence of a subdiagram containing three propagator 
lines that can simultaneously shrink to zero size, as a subset of vertices coalesce together on 
one of the two world-lines. Examples of such UV dangerous diagrams are Fig. Eli and Fig. [Hb 
[for vertices coalescing towards (t,yi{t))] or Fig. Ell and Fig. (Hi [for vertices coalescing on 
the second world-line] . The former diagrams can give poles proportional to mf m2 (with 
some velocity dependence, or some extra mass dependence due to an extra line propagating 
between the two world-lines), while the latter can give poles proportional to ml (possi- 
bly with some extra velocity or mass dependence). The reason why three simultaneously 
shrinking propagators can yield poles as e — ^ is easy to see in the approximation where the 
relativistic propagators D"^ are replaced by non-relativistic ones A~^, = — — x'|~^~^. 
Indeed, when three such propagators shrink simultaneously, the overall integral contains a 
subintegral of the form J (P'^'^x ~ J^^drr"^"^^ ~ a~^'^/(— 2^). 

On the other hand, beyond our obtaining a heuristic feeling of what are the origins of 
the poles in ai, we did not use a diagrammatic technique for evaluating the equations of 
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TABLE I: Static contributions of 'DdiX(l). All the results are presented modulo some neglected 
terms 0(e). The "principal" part of a term corresponds to the term in the IPN symmetric 
propagator D^pj^ = + c~^5^A~^, while the "retarded" part corresponds to the purely IPN 
piece c~^9|A~^. The "extra term" refers to the last term in the R.H.S. of Eq. ()4.30b|) . Note that, 
in view of Eqs. H5.21() . one must multiply the results by a factor 4 in order to get the contributions 
to the coefficients 'Dcmn in the equations of motion. 
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13 
60 


_A _ J_ + Iln 

18 12£ ^ 6 

_53 _ J_ + Iln 

90 12e ^ 6 « 







Total 


^ + J- - i In 

15 ^ 6e 3 


7 

60 


_13 _ J_ 1 Iln 

15 6e ^ 3 



motion. [Note, however, that a generalization of the (2PN level) work j45|] would lead to 
a diagrammatic technique for evaluating the Fokker Lagrangian of two point masses.] Our 
actual computations used the techniques elaborated in the previous Sections. 

We evaluated the contributions to the difference Pai coming from all the terms in the 
expression for ai deduced from ()2.19|) together with the complete expanded forms flAll|) - 
()A12|) . However, as expected from various arguments — diagrammatic analysis, existence of 
(possibly cancelled) logarithms in the corresponding c? = 3 evaluation — most of the terms 
lead to a vanishing difference Pai . The only terms that give non- vanishing contributions to 
Pai are the four terms given in Eq. (j2.23p . 

Gi[X] = ^(9,X)i, ai[r] = i^(9,f)i, 

<m = ^ I (4)i , «iM = ^ I {y^)l ■ (5.21) 

Note that, for the contributions associated to X and -Rj, one needs a IPN-accurate treatment 
of both their respective sources and the propagator . Apart from the compact support 
terms in the sources for X, T, Ri and Yi which lead to zero difference, most of the non- 
compact terms do lead to some non-vanishing contributions to the difference of acceleration 
Pai. We give in Tables Ul lIVI the contributions to T>Cmn associated to the various individual 
source terms of the "delicate" potentials X, T, Ri and Yi, which were displayed in Section 
im Eqs. ()2.12|) [of course, we limit ourselves to non-compact source terms]. In these tables, 
we use the simplifying notation 

In^ = In (gr^ri2) , In^ = In (g S2ri2) , q = Aire^ , (5.22) 

where C = 0.577. . . denotes the Euler constant. 

Summing up the separate non- vanishing contributions displayed in Tables Ul lIVl we get 
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TABLE II: Static contributions of P5if (1). 







ll^C22 




ZijdijV 

1 / (7—1 \ ^ T ?-9 o9t r 

-2 (^t^O 


_119 _ J_ 1 Xin 
900 30e ^ 15 ^^^^ 

1 Q 1 1 1 

36 6e ^ 3 '■'■'■^ 







1429 , JJ_ _ 11 , 
900 ^ 30e 15 ^^^'^ 

450 15e ^ 15 ^^^^ 







28 , J. _ 4 , 

9 ^ 3e 3 

1 Q 1 In 

— 12. i kiln 

36 6e ^ 3 







Total 


_33 _ J_ 1 2 1 

50 5£ ^ 5 


397 1 3 3 1 
300 lOe 5 ^''^r 


31 I J_ _ In 
12 ^ 2e ^'^s 



TABLE III: Static contributions of The "principal" part, "retarded" part and the "extra 

term" have the same meaning as in Table Q] 







\VC22 




d—l fi T//) 1/ 1 

d—2 ^ * fc 1 principal 

~ 1^ '^^i^k\ retarded 
d{d-l) p. 

i{d—2)'^ * * Iprmcipal 
4(d-2)^ f^t^l^i^ retarded 

extra term 


_31 _ J_ I 2i 
18 3e ^ 3 ^'^'^ 

43 , J_ _ 2 1 
18 ^ 3e 3 ^^^^ 

5 I J_ _ 1 In 
4 4e 2 ^^^-r 

_I _ J_ 1 1 In 
4 ie^ 2 ^^^^ 

1 

4 


_31 _ J_ I 2 1 
18 3e ^ 3 ^'^^ 

43 , J_ _ 2 1 
18 ^ 3e 3 ^'^^ 

5 I J_ _ 1 In 
4 4e 2 ^^^J- 

_7 _ J_ , 1 l„ 
4 4e 2 ^^^f' 

1 

4 









Total 


1 

12 


1 

12 






the following total differences 

22 44 1 02 

2^C3i = -TT + ^ Hqr[ ri2) - + , (5.23a) 
6e 6 5 

Pc22 = 9 + 0(e) (5.23b) 

22 44 "ifiS 
-Dcr, = - \n{qs2n2) + + 0{e) . (5.23c) 
36 3 15 

Finally, adding Eqs. ()5.23|) to the pure Hadamard-Schwartz result fl5.11|) . we get the dimen- 
sionally regularized static contributions to ai: 

22 44 q5Q 
4 = + y l^(^-?2) - ^ + 0{e) , (5.24a) 

41 

C22 = 175 -- 7r2 + 0(£), (5.24b) 
lo 



dr 22 44, , 2 X 1799 41 . ^ x 

^13 = ^ - y In(grL) + - ^ + ^(^) • (5-24c) 
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TABLE IV: Static contributions of 





1 -T-> 


1 -T-> 


T6^Ci3 


WkidkiVi 











-I (s^j dtWikdkV 


Hii _L _ — In 

18 ^ 3e 3 


107 , 5 5 1 
45 12e 6 ^'^^'^ 


71 , 1 1 I 
_LL _L _I_ _ In 

450 ^ 60e 30 


diWki dkVi 


257 1 7 7 1 
450 ^ 60e 30 ^^^^ 


149 2 1 4 I 
— — — =— 4- — n 

225 15£ ^ 15 '■'■'■^ 


71 1 1 1 ] 

450 60£ ~ 30 


-dkWiidiVk 


450 60£ ^ 30 '-'■'■r 


14Q 9 zL 1 

J £ ^ In 

225 ^ 15£ 15 ^^^^ 


71 1 In 

Lt i ^ In 

450 60e ^ 30 


d—1 O T/^l D 


2681 1 19 19 1 
900 30£ 15 


3791 1 53 53 i 
900 60£ 30 ^"■'T 


71 1 1 1 In 
450 60e 30 ^'^^ 


-^^Vkd.vdkV 


100 lOe ^ 5 ^^^^ 


-| - i + 21n^ 


100 lOe ^ 5 « 




- ^ + ln 
4 2e ^^^^^ 


43 1 J. _ 4 1 
25 ^ 5e 5 '■^'■r 


^ 1 J_ _ 1 In 
100 ^ lOe 5 ^^^■^ 


-lilE^yVdkVdkV. 













-| - i +21n^ 


1 ^1 

-| - i + 21n^ 





9Vi 81 fl^V- 


n 


n 


n 


^\d-2r'^ 9kmvk 

d{d-l){d-3) a T/Q TV- 



















Total 


-M. _ _3_ _L 3 J 
100 lOe ^ 5 ^ 


100 lOe ^ 5 '^'■'■r 






As expected the two Hadamard regularization length scales r[ and S2 have cancelled be- 
tween and Vcmn to leave a result which depends only on the dim. reg. regularization 
parameter e = d — 3. One might be surprised by the presence in cjj^^ of terms ± y ln(r^2) 
compared to corresponding terms ±y ln(ri2) in cf^, and by the absence of any adimension- 
alizing length scale in these logarithms of r^2- These two properties can be understood when 
one remembers from the discussion above that the coefficients which have the same physical 
dimension the combinations (>^^^'"''>'^ q<^^_ in the present case, this means that 

are dimensionless. It is easy to see, thanks to the pole terms =f22/(3£:) and the 
expansion = exp(45 ln£o) = 1 + 25 In^Q + 0{e'^), that the combinations cf-^ and cfg 
do indeed depend only on the dimensionless quantities e and ln(r^2/^o)- 



VI. RENORMALIZATION OF THE EQUATIONS OF MOTION 
A. Poles in the dimensionally regularized bulk metric 

The first computation of the dimensional continuation of the 3PN gravitational interac- 
tion of point masses was done in ADM coordinates and resulted into a finite (z. e. , without 1/e 
poles) answer js^ . Our task in analyzing the physical meaning of the harmonic-coordinates 
result ()5.24|) is to interpret the presence of 1/e poles in it. For this we have to remem- 
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ber that, as in Quantum Field Theory (QFT), dimensional continuation is a regularization 
method which, like all regularization methods, transforms truly infinite results, say contain- 
ing Jq^^ (P:s./rf, into finite, but "large" ones, which depend on some cut-off parameter, e.g. 

J^^^ d?^/r\ = 47rln(ri2/si) or JJ^^ d^'^^:s./rf"^~^'^'^ oc 1/e. Any regularization must be followed 
by a renormalization process which allows one to absorb the cut-off dependent terms in 
some of the basic hare parameters of the theory. 

In order to have a clearer understanding of the poles in the (static) equations of motion 
(|5.24p [we shall prove below that our discussion extends to the full, velocity-dependent 
equations of motion], we need to analyze the presence of poles in the "bulk" metric, i.e., 
the metric g^v{^\ yi, y2) evaluated at a generic field point x, away from the two world-lines. 
Indeed, if we were considering the gravitational field generated by regular (z.e., non point- 
like) sources, a complete physical description of their gravitational effects would necessitate 
the simultaneous consideration of the bulk metric and of the equations of motion of the 
(extended) sources. Similarly, in the present formal study of two-point-like sources, we 
need to consider both the equations of motion = aa(yi, y2, vi, V2) and the bulk metric 
Qi^v (x; yi (t) , y2 (t) , vi (t) , V2 (t)) . 

It is here that the diagrammatic representation introduced above plays a useful role in 
highlighting the structure of divergencies in the equations of motion and in the bulk metric. 
Indeed, it is clear that the divergent diagrams of the equations of motion of the first particle, 
where the 1/e pole is due to the presence of a subdivergence induced by three propagators 
shrinking onto the second world-line (such as in Fig. Eli or Fig.^Ji), will correspond to similar 
1/e poles in the bulk metric, for the corresponding "bulk diagrams" where the special point 
marked by a cross in the diagrams above (denoting the coincidence x = yi) is detached from 
the first world-line to end at an arbitrary point in the bulk, as indicated in Fig. [HI 



X 




X 




X 




X 




































1 1 


i; 












1 1 


i 







12 12 12 12 

abed 



FIG. 6: Some divergent diagrams for the bulk metric. Here, contrary to the previous Figures, the 
field point, labelled by an x, is detached from the first worldline. 

Evidently, in addition to such diagrams as Fig. IHK and which will contain (at least) 
a factor m^, there will exist "mirror diagrams", containing a factor mf, and obtained by 
exchanging the labels 1 and 2. On the other hand, note that the bulk poles oc of the 
type of Fig. ^ and IHli do not (necessarily) correspond to poles in the equations of motion 
of yi because their coincidence limits x — yi induce diagrams of the type of Fig. or 
EJi containing four shrinking propagators instead of three. [Though such diagrams would 
exhibit worse divergences in a dimensionful cut-off regularization schemes, they are generally 
non dangerous in dimensional regularization because the integral J rf^+'^x/r^^^'*'^'' has no pole 
as e ^ 0.] 

A careful analysis of the possible presence of poles in the various potentials V,Vi, - ■ ■ ,T 
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we use to parametrize the bulk metric (aided by the structure of potentially dangerous terms 
sketched in Fig. IH)) shows that, at the 3PN approximation,^^ such poles can only be present 
in the IPN-level exp ansion of X and in the Newtonian-level approximation of T. Drawing 
on the results of [22] and 37 1 we can also see that all velocity-dependent terms in the poles 
present in Xipn and T [i.e., the terms proportional to m\vl or exactly cancel in 

the combination 4X/c^ -[- 16 T/c^ that matters for the bulk metric. [This shows up, for 
instance, in Eq. (7.1) of [i^ which implies that the divergencies linked to the second world- 
line, characterized by the presence of lns2, do not depend on velocities. This shows up also 
in the absence of poles in the velocity-dependent contributions to ai proportional to m\, see 
Eqs. fl6.47c|) - ()6.47d|) below.] We are therefore left with evaluating the poles present in the 
static limit (vi, V2 — > 0) of Xipn and T. Clearly, from Fig. IHl the poles in Xipn and T will 
come only from cubically non-compact (CNC) sources. Finally, as we are only interested in 
the pole part we can neglect the d-dependence of the coefficients in the sources of X and T 
(which we indicate by using a symbol ~). Thus, these poles can only come from 



X, 



CNC ^ r-|-l 
static 

CNC _ Q-l 



static 



d,,vn-\-diVd,v) 
1 

2 

'1 

2 



d'f V + dij V n-^ {-2di Vdt Vj) 



(6.1a) 



(6.1b) 



The static poles (involving factors m\ or m|) in Eqs. 1)6.11) are then obtained by: (1) consid- 
ering sources involving three times Vi or three times V2 (where Va denotes the piece oc iria 
in V)^ (2) evaluating the time derivatives in the static limit, using for instance 



(6.2) 



and (3) expanding up to the required accuracy the (time-symmetric) propagators according 
to = A-i + c-2 92 A-2 + 0{c-^). 

As an example among the simplest terms, let us consider the m\ contribution coming 
from the first term on the R.H.S. of Eq. ()6.1b|) . 



static(l) 



\ A-i[V;^ d', l^ilstatic = -\ ^-'M < V,] = a{ d, ^-'[Vl^] 



(6.3) 



Using Vi = 2^Gkmirl-'^ + 0[c'^) ~ Gm^rl^-' + 0(c-2) and r\ = r^^y[{\ + 
2) (A -|- d)], one finds that the pole part of ()6.3p reads 



static(l) 



I2e 



(6.4) 



Similarly an analysis of the second source term in Eq. ()6.1b|) yields 



G^ ml a{ dj r[ 



static{2) 2^ 



-l-3e 



(6.5) 



I.e., at order c 



in goo, c ^ in go*, and c 
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so that the full (static) contribution of T is 



f::L^^G'mlaid,r^'-'^. (6.6) 

The analysis of the pole part in the static limit of X, Eq. ()6.1a|) . is more intricate because one 
must expand to IPN accuracy both Va — Gnia r^^^^ + Gnia r^~^ and the propagator 
\Jr^. This yields 

^sTL--Y^^-?«i5.-r^-^^ (6.7) 

Let us now consider the improved \^-potential ()2.10|) that makes up the essential part of 

2 /rf-3\ _ 4X 16f _ 4X 16t 



V = V -^\——]K + — + — -r:^V + — + —-, (6.8) 

such that ^00 = - exp(-2V/c2)[l - %ViVil& - ?>2RiVi/c^ + 0(1/0^°)], see Eq. flTm . 
Combining the results above, we find that the only 1/e poles in the bulk metric g^ui'^, Yi, y2) 
show up in goo at the 3PN level and are (when expressed in terms of the improved potential 
V, and after cancellation of m?^v'^/e terms between X and T of the following static form 



V(x, yi, ys) = V + + \Vi + \ 



C2 c4 c6 



(6.9) 



where V2, V4,Vg's are finite when e — 0. To understand better the structure of result ()6.9|) 
let us introduce the notation 

11 Glml 
3e c6 



C = +T;-^V^«a, (6.10) 



where Gat is the 3-dimensional Newton constant and a\ the d-dimensional acceleration of 
y^. [This definition ensures that has the physical dimension of a length.] In terms of the 
definition ()6.1U|) . the result ()6.9|) can be equivalently written as 

V(x, yi, ys) = [k(x - y.) + Ci K(x - y,)] + ^V^ + ^V, + ^V, , (6.11) 

a 

where the pole part is entirely contained in the terms proportional to and C2 [^6 here 
differs from Vg in Eq. ()6.9p by some finite corrections when e ^ 0]. The fact that poles 
appear only in V, at order c~^, implies that there are no divergencies in the harmonic gauge 
conditions (j2.13|) in the bulk. Indeed, (j2.13ap needs V at order c~'^ only, and ()2.13b|) at 
Newtonian order only. 



B. Renormalization of poles by shifts of the world-lines 

Result (j6.1ip indicates a simple way of renormalizing away the poles present in the bulk 
metric. Indeed, the logic up to now has been to describe in the simplest possible manner 
a gravitationally interacting two-particle system, parametrized by the following bare pa- 
rameters: G^^^^, m}l^^^, 1712^^^, y\^^^, y2^^'^, considered in everywhere-harmonic coordinates, 
-pA ^ ^a/3pA _ g j^-^ particular, the internal structure of each particle has been, up to now. 
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entirely described by a monopolar stress-energy distribution, i.e., T^^ oc m}^^'^ — y^'^'^*^). 
In other words, we have set to zero any higher multipolar structure. Eq. ()6.1H) is most 
simply interpreted by saying that the non-linear interactions (see Fig. ^ dress each particle 
by a cloud of gravitational energy which generates, at the 3PN order, a divergent dipole 
in the Newtonian-like potential. Therefore, to get a net, finite bulk gravitational field we 
must endow each initial particle by an infinite, bare dipole, corresponding to a counterterm 
AT^^ oc —Tn]!^^^^Cl 9j ~ Ya), which will cancel the non-linearly generated one (jfj.llj) . An 
equivalent, but technically simpler way of endowing each particle by a bare structure able 
to cancel the dipolar pole terms in ()6.11|) is simply to say that the central bare world-lines 
used in our derivations up to now, henceforth denoted as y^^'''^, can be decomposed in a 
finite renormalized part yjj™ and a formally infinite shift involving a pole oc 

Ya =ya +?a, V„ = V„ + • (6.12) 

The gravitational potential of two point particles [oc 5(x — ya^'^'^)] is then 

v(x, y^-, y^<=) = J2 [^«(^ - yr - + a k(x - yr - ^a)] (e.is) 

a 

+ i V2(x, yr + \ V4(x, rr + Ca) + 4 ^^^^ yr + ^a) . 

c c c° 

Assuming that the vector is of 3PN order [i.e., = C'(l/c^)]) we can rewrite Eq. 



as 



V(x, y^^ y^^^^) = J2 [Vai^ - yD + id - a) K(x - yr)] (e.M) 

a 

+ 3 V2(x, yr ) + i V4(x, yr) + ^ V6(x, yr) +0^^ 



which makes it clear that the potential will be finite (at 3PN accuracy) when e — if we 
choose 

Ca = Ca + 0(^), (6.15) 

where by 0{e^/c^) we mean a term finite when e and of the 3PN order. We shall 
henceforth refer to in Eq. ()6.15|) as a shift of the a^^ world-line. The reasoning above 
shows that the introduction of such shifts, at the 3PN order and having the pole structure 
fl6.10|) . is necessary to renormalize away the poles present in the bulk metric. It remains 
to show that these shifts are also sufficient to renormalize away the poles present in the 
equations of motion. 

The effect of 3PN-level shifts on the equations of motion is easy to obtain. Indeed, 
the equations of motion we computed above concern the original, bare world-lines y^^'^'^. For 
the first particle, they had the structure (in dimensional regularization) 

"bare dr/ bare bare bareN 

yi — ^1 iyi2 '^1 '^2 ) 

— q'^'' -L Q<^'' f^r^^^'^ bare ,^bare\ 

— ^Nliyi2 J"r^lPNliyi2 5^1 ) 

I dr / bare bare bare\ i dr / bare bare bare\ 
"T ^2PNliyi2 ^^l 1^2 J "T ^2.5PNliyi2 i^l 1^2 ) 

+ a^^Nl(yl2^^v^^v^^), (6.16) 
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where yif^*^ = y^'^^'^ — yg'^''*^. Here denotes the dimensionally continued Newtonian-level 
acceleration, 

a^f ^ = di V2iyu) = fGm2~kd, rl^" , (6.17) 

where by a slight abuse of notation we pose di = d/dy[2, where G = GnC-o denotes the 
(i-dimensional gravitational constant, and where the d-dependent correcting factors 

d-2 l + e ~ V{^\ V{^\ 

tend to 1 as e ^ 0, but will play a significant role below. 

When inserting the redefinitions ()6.12p into ()6.16p one easily finds that the renormalized 
equations of motion, i.e., the equations for yj,""^, read [using only = 0{l/c^) at this stage] 

yr = ar(yir,vr,vr), (6.19) 

where 

ar(yi2, vi, V2) = af (yi2, vi, V2) + 6^ a^iy^2, Vi, V2) + O (^^^ , (6.20) 

with 

ai(yi2, vi, V2) = {^1 - ^i) dj 4\iyi2) - ii ■ (6.21) 

Let us note that the effect on the equations of motion of a (3PN) shift of the world-lines, 
Eq. ()6.21|) . is technically identical to the effect on the equations of motion of the restriction 
to the world-lines of a 3PN-level coordinate transformation, say x^g^ = x* — e*(x, t) and 
^new = t — c^^e°(x, t). Indeed, a coordinate transformation has two effects: (i) it changes the 
bulk metric by 6^ g^vix) = g^iui.^), where denotes the Lie derivative along e^, and (ii) it 
induces a shift of the world-lines y^new — ya~^\'^^~^^^aya (plus non-linear terms in e^), where 
we denote the coordinate change at the location of the a-th particle by ejj(t) = [e*(x, t)]x=y<j. 
Because of the diffeomorphism invariance of the total action, the effect (i) does not change 
the action, so that the net effect of a coordinate transformation on the equations of motion 
reduces to the effect fl6.21|) of the following shift induced on the world-lines: 

C = 4 - c"^e° + non-linear terms . (6.22) 



The coordinate transformations considered in j22], see Eq. (6.11) there, were of the type 
^iii.^) = i9^(X]a '^aAa) ' whcrc the Ka's are some coefficients, so that we see that the 
latter induced shift reduces at the 3PN order to the (regularization of the) purely spatial 
coordinate transformation evaluated on the world-line: 

ea = < + (^] . (6.23) 



We have checked that the formula given by Eq. (6.15) in 2^1 for the coordinate transforma- 
tion of the acceleration of the particle 1 gives exactly the same result as the one computed 



Actually, one should consider, as e.g. in 45], a non-generally-covariant gauge-fixed action. But the 
"double zero" nature of the gauge-fixing term, say cx ^J—gg^^i, , ensures that it does not contribute 
to first order in t^^. 
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from the effect of tlie sliift ()6.21|) . [Tlie agreement extends to d dimensions if we consider 
tlie straightforward extension of the latter coordinate transformation to d dimensions, 

namely e^{x) = df,{J2a '^'^^ ■] — 

Note that the coordinate transformations e^{x) were considered in 22] only in terms of 



their effects, Eqs. (|6.22|) - (|6.23|) . on the equations of motion. This was sufficient to prove, for 
instance, that the two constants r[ and r!^ in the 3PN equations of motion are not physical, 
because they can be gauged away in 3 dimensions and therefore will never appear in any 
physical result. However, we remark that the extension to d dimensions of the coordinate 
transformation e^{x) of the bulk metric, say e^(a;) = 9^ Ka A; (with coefficients 
oc s , as needed to remove the poles in the equations of motion), does not lead to a 
bulk metric free of poles. Indeed, assuming oc e~^, we see that the pole in the spatial 
coordinate transformation ej(x) would then induce a pole in the spatial components of the 
metric, 6^gij = d^ej + djei + ■ ■ ■, but this is inadmissible because we have proved above 
Eqs. that, at 3PN order, only the time-time component of the bulk metric contained 
a pole. A bulk coordinate transformation of the type above can then remove the poles 
in the time-time component of the bulk metric only at the price of creating a pole in the 
initially pole-free spatial metric. We shall leave to future work a complete clarification of the 
possibility of using, within our dim. reg. context, a coordinate transformation to induce the 
shifts ()6.15p . For the time being, what is important is that our introduction above of shifts of 
the world-lines (a priori unconnected to any coordinate transformation) is a consistent way 
of renormalizing away the poles in the metric, and that its effect on the equations of motion, 
Eq. (|6.2ip . is identical to the transformations of the acceleration obtained in Ref. 2^ . 



It remains now to show that the same world-line shifts (j6.15p that renormalize away 
the poles in the bulk metric, Eq. (|6.14p . do renormalize away also the poles present in the 
original bare equations of motion [see ()5.24j) for the static contributions and ()6.47|) below 
for the kinetic ones]. For this purpose let us consider a shift of the more general form 

e = ^^^<. (6.24) 

where e{d) represents a certain numerical coefficient depending on d, and where denotes 
the d-dimensional acceleration of given by its Newtonian approximation ()6.17p [but, for 
notational simplicity, we henceforth drop the label on such accelerations entering 3PN 
effects]. Inserting (j6.24j) into (j6.2H) yields (for the index a = 1) 



{ml — m'2) 02 dj a\ — ml v-^-^ djk a\ , (6.25) 



where = v{ — v{ and v-'^ = fi2't^i2 [and also, as before, dj = d/dy{2]. Before further 
evaluating (j6.25p by inserting the explicit expression (j6.17|) for the acceleration, we shall 
consider some simple but important consequences of the structure (j6.25p . 



C. Link to the general class of harmonic equations of motion 

As recalled in the Introduction, previous work on the 3PN equations of motion in har- 
monic coordinates has shown that these equations necessarily belonged to a three-parameter 
class of equations of motion, say 

yi'='^ = af (yi2, VI, V2; A, r;,r^). (6.26) 
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The dimensionless parameter A could not be determined by the previous work in harmonic 
coordinates. However, comparison with the work in ADM coordinates, has shown [i^, 113 
that, if there were consistency between the two calculations one should have the following 
link between A and the corresponding ADM "static ambiguity" parameter a;^, 

3 1987 , 

A = ujs . 6.27 

11 3080 ^ ^ 

If dimensional regularization is a fully consistent regularization scheme for classical pertur- 
bative gravity, we then expect that the dim. reg. determination of in ADM coordinates 
[35l |. namely o;^"' ^dm _ should lead to a dim. reg. direct determination of A (in harmonic 
coordinates) of harmomc _ —1937/3080. We will turn to this verification in a moment. 

The two other parameters, denoted above r[, rg, entering the general "parametric" har- 
monic equations of motion ()6.26|) have the dimension of length and have the character of 
gauge parameters. Indeed, they can be chosen at will (except that one cannot set them 
to zero) by the effect of shifts of the world-line, induced for instance [but not necessarily, 
cf. a discussion in Section IVIBI above] by some gauge transformations. In the way they 
were originally introduced (2^], the two parameters r[ and r'^ can be interpreted as some 
infinitesimal radial distances used as cut-offs when the field point tends towards the two 
singularities yi and y2. Therefore in principle Inr^ and Inrg should initially be thought of 
as being (formally) infinite. However, it is trivial to show that by a (formally infinite) gauge 
transformation, involving the logarithmic ratios ln(r^'/r^) and In(r2/r2), where r'l and r!^ 
denote any two finite length scales, one can replace everywhere r[, r'2 by the finite scales 
r'/, By this process it is therefore possible to identify the two sets of scales and thereby 
to think of the scales r^, r'2 as being in fact finite, as implicitly done in Ref. i^]. In the 



language of renormalization theory, the original (infinitesimal) scales and r'2 would be 
referred to as Hadamard-regularization scales entering the computation of divergent Poisson 
integrals [see Section llll Bl above] . while the (finite) scales r'/ and r'2 would be referred to as 
the arbitrary renormalization scales entering the final, renormalized harmonic-coordinates 
equations of motion. In the present paper, in order to remain close to the notation used in 
"22 1 ■ we shall keep the notation r'^ and r'2, but interpret them as arbitrary finite constants. 



which means that we shall identify them with the finite renormalization length scales r'/ and 
r'2. In other words, the scales r'^^, r'2 used in the present Section should in principle be distin- 
guished from the scales r'^^, r'2 used in Section llll Bl above. [Remember, in this respect, that 
the regularization scales r^, r'2 have disappeared when computing the dim. reg. equations 
of motion.] 

With our notation, and still focussing on the static contributions to the equations of 
motion, the "parametric" equations of motion ()6.26|) imply the following structure for the 
static coefficients Cmn'- 

li,.(rH)_||, 

Rp,, , 34763 44, 41 2 ^01 n 

4'^(A,ri,r^) = ^YT'Y^^Te ' ^^^^ ^ 

RF/^ , ,x 44, fri2\ 10478 44 , 41 „ / x 

It will be convenient to replace the parameter A by the parameter cu^, using ()6.27|) as a defin- 
ing one-to-one map between A and 0;^. With this change of notation the static coefficients 
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become 



c?fK,r;,r;) = fln(l^)-5|Z. (6.29a) 

41 

c^i{u;s,r[,r',) = 175 + AlUs-—7t\ (6.29b) 

lb 

RF/ / /X 44, /ri2\ 110741 , 41 o / x 

cfs'^KK,^;) = -yln(^^j +-^ + 4.;,--7r^ (6.29c) 



Note that there are two combinations of the three coefficients which do not depend on 
^22^, and the combination cff + cfl 



lnri2, namely C22, and the combination cW + cW, or even better the combination 



C3Y + 4^-4^ = yln(^)-^, (6.30) 

which depends neither on lnri2 nor on ujg (or A), and which contains, hke for C22, simpler 
looking rational numbers. 

We now come back to the effect of the general shift ()6.24|) on the dim. reg. equations of 
motion. Let us ffist focus on the static terms. We recall that the (dim. reg.) renormalized 
equations of motion had necessarily the form ()6.20p . By projecting the latter equation along 
the static terms Cmn, with m + n = 4 [recalling Eq. (|5.6|) ]. it will induce a result for the 
renormalized static coefficients of the form 

O = Cl(e) + %)C^n , (6.31) 

where the 6^(^^)CmnS are the static coefficients corresponding to S^a{, Eq. ()6.2H) . When 
choosing Cai^) of the form ()6.24p . we see from Eq. ()6.25|) that S^Cmn is simply obtained 
by projecting the ffist term on the R.H.S. of ()6.25p . Remembering that a2 oc mi and 
ai oc m2, we see that the latter term contains the factor {ml — m^) rrii m2 = m2 — mi m\. 
Therefore, without doing any further calculation, we see that the shifts 5^ Cmn have the 
special properties: 6^ C22 = and 6^ C31 + 5^ C13 = 0. In other words, a shift of the world-lines 
of the type (|6.24p leaves invariant both C22 and the combination C31 + C13 (as well therefore 
as the combination C31 + C13 — C22 considered above). As a consequence, we can compute 
without effort from our previous regularized (but unrenormalized) dim. reg. results ()5.24|) 
the following two combinations of the C,^{e)-renormalized static coefficients: 

clf = 4 = 175-^71', (6.32a) 

d d d 

^31 + ^13 ~ ^22 ~ ^31 + '^13 ~ '^22 ~ ~^ • (6.32b) 

By comparing ()6.32a|) with Eq. ()6.29b|) we discover that our present calculation using 
dimensional regularization in harmonic coordinates necessarily implies that 

1987 

Us = A = . 6.33) 

3080 ^ ' 

This nicely conffims the previous determination of Ug by a dim. reg. calculation in ADM- 
type coordinates [3^. We think that our present harmonic-coordinates dim. reg. result 
calculation is important in proving the consistency of dimensional regularization, and thereby 
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in confirming the physical significance of the result ()6.33p . A recent calculation 3^ 3^ has 
also confirmed independently the result ()6.33p by means of a completely different method 
based on surface integrals, and aimed at describing compact (strongly-gravitating) objects. 
By comparing ()6.32b|) with (|6.3(J|) . we further see that 

r[ = r'2 [in the case of the dim. reg. shift ()6.24j) ] . (6.34) 

Contrary to ()6.33|) which represents the determination of a physical parameter (having an 
invariant meaning), the result ()6.34j) has no invariant physical significance. Eq. ()6.34j) is 
simply a consequence of our particular choice for the shift vector (j6.24|) , in which we assumed 
that e{d) is a purely numerical coefficient, independent on any properties indexed by the 
particles' labels 1 and 2. In summary the particular shift (|6.24j) yields some equations of 
motion which are physically equivalent to a subclass of the general equations of motion 
considered in ji^], characterized by the constraint ()6.34|) . 

Next we relate the common length scale ()6.34|) to the basic length scale io entering 
dimensional regularization. For doing this we need to fully specify the value of the shift, 
i.e., to choose a specific coefficient e{d) = e(3) + ee'{3) + O^e"^) in Eq. ()6.24|) . We already 
know from Eq. (j6.1(J|) that the coefficient e{d) in Eq. (j6.24|) must tend to 11/3 when d ^ 3, 
if the ^-shift is to remove the poles in the bulk metric. As in quantum field theory (QFT) 
we could then define the Minimal Subtraction (MS) shift as 

Cms - ^ ■ (6.35) 

3e c° 

However, as is well-known in QFT, such a MS subtraction has the unpleasant feature of 
leaving some logarithms of vr and the Euler constant in the renormalized results. These 
numbers come from the expansion of the Gamma function and the associated dimension- 
dependent powers of tt entering the d-dimensional Green function. In our context, these 
numbers showed up in Eq. fl5.24j) in the guise of the combination 

ln(g) = ln(47re^) = C + ln(47r) . (6.36) 

Like in QFT, this leads us to consider the following modified minimal subtraction (MS) shift, 

_n G%~k^ml 

which differs from the MS shift ()6.35|) by the explicit factor of k"^ it contains, and by the use of 
the (i-dimensional (Newtonian) acceleration given by Eq. ()6.17p . The inclusion of two explicit 
powers of k in the coefficient e{d) entering Eq. ()6.24p . i.e. the definition ejj^{d) = y^^' 
means, when remembering that aj^, Eq. ()6.17|1 . contains one power of k, that the static terms 
in Eq. ()6.25|) will have four powers of k and the kinetic terms three. The overall factor k^ in 
the static terms is natural because these terms are of order and the x-space gravitational 
propagator in d dimensions always includes the combination G/c |x — x'^^'^. Finally, using 
the fact that the expansion of k{d) near d = 3 reads 



T(^) 1 

k{d)= ^,1^ =l--e\nq + 0{e'^), (6.38) 



TT 2 



2 
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it is easy to see that the MS shift defined by ()(i.37j) will cancel the Ing in the bare dim. reg. 
results of Eq. (|5.24|) . Finally, we find that the evaluation of Eq. (j6.25|) for the specific MS 
shift, given by e^jg((i) = y A;^, yields for the MS-renormalized static coefficients 

c^/ = 175 - ^ vr^ (6.39b) 

where the reader can note that the ln(r^2) entering the bare dim. reg. result ()5.24|) have 
been transformed into ln(ri2/^o) through the e-expansion of the factor r^2~^'^ present in 
5^ a\. 

We already discussed above the comparison of two simple combinations of the dim. reg. 
results ()6.39j) with the Hadamard-regularization results ()6.28j) . It is easy to see that the 
remaining independent combination, say (|6.39a|) . is fully consistent with its counterpart 
fl6.28ap . and allows one to relate the basic renormalization length scale io entering dim, reg. 
to the common length scale ()6.34j) entering the general equations of motion of |22l |: 

\ \ 327 

In I -i ) = In I ) = (for the MS renormalization) . (6.40) 

\to/ V^o/ 3080 

Evidently, the precise values one gets for r[ and depend on the precise choice of the 
compensating shift. 

Let us now remark that in fact one can recover exactly, provided of course that the 
crucial result ()6.33j) holds, the general "dissymmetric" class of equations of motion of |2^ . 
i.e., the general parametric result (j6.28j) or (j6.29j) with r[ ^ r'^. For this purpose it suffices 
to consider a slightly more general shift than the one assumed in the MS regularization. 
Namely, consider a shift of the same form as ()6.24|) . but in which one allows the (i-dependent 
coefficient e{d) to depend on the label of the particle in question, that is 

C = 5=f^<. (6.41) 

where now ei((i) and 62 (rf) are allowed to be different from each other. The most general 
way of parametrizing such dissymmetric ea((i) [however constrained by 6^(3) = 11/3] is 

eM) = yP [1 - 2epa + 0{e^)] , (6.42) 

with two independent numerical coefficients pa- It is then easily checked that the shift 1)6.4111 
defined by the particular choice 

= In f ^ V ^ , (6-43) 
^ \eoJ 3080' ^ ^ 

transforms the dim. reg. equations of motion into the general (r^-dependent) family of 



solutions obtained in |2^. If we suppose that the constraints f|6.40|) hold, then pa = and 
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we recover the shift assumed in the MS regularization. On the other hand, note than one 
can reach even more general classes of renormalized harmonic equations of motion in dim. 
reg. (as one could have also done in Hadamard regularization). Indeed, we could use the 
freedom indicated in Eq. ()6.15p of adding arbitrary finite parts to the shifts. Anyway, the 
result that the shift ()6.4H) - ()6.43|) gives equivalence between the dim. reg. and the extended- 
Hadamard 3PN accelerations [we check in Section FVI Dl below that the kinetic terms work 
also], constitutes the main result of the present paper (Theorem 2 in the Introduction). 

The two approaches we have discussed here are of course equivalent: choosing some dim. 
reg. basic length scale io and some specific, simplifying dim. reg. shift (such as the MS 
one), and then determining the values of the scales for which the dim. reg. results match 
with the Hadamard reg. ones; or arbitrarily choosing some Hadamard scales and then 
determining the corresponding general dissymmetric dim. reg. shift ()6.4H) - ()6.43p . in terms 
of the chosen r^'s. What is important is that we have checked that the three renormalized 
dim. reg. static coefficients ()5.24|) are fully compatible with the three extended- Hadamard 
reg. static coefficients ()6.28|) or ()6.29|) . and that their comparison yields one and only one 
physical result, namely: A = — ^kI^- 



D. Kinetic terms and check of the consistency of dimensional regularization 

Up to now we have verified the following aspects of the consistency of a dim. reg. 
treatment of the 3PN dynamics of two point particles: (1) consistency between the shift 
(I6.10|l needed to renormalize the bulk metric and the shift ()6.37|1 [or ()6.4H) - ()6.42|1 ] needed to 
renormalize the equations of motion; (2) consistency between the three finite, renormalized 
dim. reg. static coefficients ()(j.39j) and the general three-dimensional ones ()(i.28j) . [j^ : and 
(3) consistency between the present dim. re g. v alue of A and the previously derived dim. 
reg. value of tOg in the ADM-Hamiltonian [35|- It remains, however, to check that the 
velocity-dependent terms in the renormalized dim. reg. equations of motion do agree with 
their analogs in the harmonic-coordinates equations of motion of ^] . This will in particular 
prove that the dim. reg. equations of motion are Lorentz invariant. 

In the notation of Eq. ()5.6|) above, we need to consider the values of the velocity-dependent 
coefficients C2i(vi, Vs, riia), c'oifvi, Vg, riis), ^'^(vi, Vg, riis), Co3(vi, Vs, riis), Co3(vi, Vg, riis), 
and Cq3(vi, V2, 1112). In Ref. j2^, they were shown to take the following parametric forms, 
which actually depend only on the regularization scale r[ but not on r'^ nor A: 



BF, 00 [2 fn2\ , 48197 2 36227, , , 36227 

4 K) = -22 [vt2 - 5(ni2fi2)'] In (^—1 ^^^v\ - —^[viV2) + 



840 ' 420 ' ' 840 



45887, ,2 24025, , 10469, ,2 
[^(^i2fi) ^ ^^[nx2Vx)[nx2V2) -^[nx2V2) , (6.44a) 

4f (ri) = cr (K) = -44(n,2.i2) In + ^(-12^1) " ^(-12-2) , (6.44b) 

43 

(^l = l8iviV2) ~9vl- inx2Viy + 2(ni2i;i)(ni2i;2) + y(ni2t^2)' , (6.44c) 

c'of = = 4(ni2t;i) + 5{nuV2) • (6.44d) 

As explained in Eq. (j5.H) . the dim. reg. expressions of these coefficients can be computed 
as the sum of pure Hadamard-Schwartz (pHS) contributions, and c'^^, and the "dr — 
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pHS" differences, Vcmn, l^c'^n- The calculation of the pHS contributions has been explained 
in Section Fill Dl above, and we get the following results from Eqs. ()3.55p and ()6.44j) : 



.pHS 

--21 



JpKS 
-21 



."pHS 
-21 



„PHS 
•-03 

/pHS 

•-03 

//pHS 

-03 



-22 [vl, - 5(^12^^12)'] In 



^"12 



10639 



50885, ,2 1892, 

■(^i2^^i) H :^{ni2Vi){ni2V2 



168 ' 60 
3325 



5879, , 5843 . 
(^1^2) + ^2 



120 



168 



-44(ni2t;i2) In (^y-^ + 
-44(ni2t;i2) In (^y-^ + 



7279 
~8F 
5189 



1^12^^1. 



(^12^^1, 



12 

5879. 

5843 



("'12^^2)' 



1^12^2 J 
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(^12^^2 



311 



18(fii;2) - y'"2- (nuvi) + 2(ni2t;i)(ni2f2) + ^(?^i2^^2) 

4:{ni2Vi) + 5(r2i2t;2) , 
37 

^ni2Vi) + y (ni2f2) . 



(6.45a) 
(6.45b) 

(6.45c) 

(6.45d) 
(6.45e) 
(6.45f) 



Secondly, the method explained in Sections IIVBI and IV CI for computing the differences 
VsLa is found (after doing calculations similar to those reported in the Tables of "static" 
contributions above) to lead to 



21 



VCn 



21 



Vr" 



03 



11 33, ,_ 
— - y In (gri 



_ /4/3 2/3 



12 



r 2 r' / \2i 499 2 359 , ^ 



184 



15 
22 

e 

22 

e 



2957, ,2 425, , , , 217, ,2 a 
[ni2Vi) + —{ni2Vi){ni2V2) —{ni2V2) + 0{e) 
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33 In (^q r'^^^r^^ 
33 In (^r'i^^rl2 



, 359, , 368 

(^^12^12) + — (ni2t^i) - —I 
15 15 



ni2^^2) + 0{e) 



7 

0. 



-Xni2V2) 



-T^{ni2V2) ■ 



(6.46a) 
(6.46b) 

(6.46c) 

(6.46d) 
(6.46e) 
(6.46f) 



Together with Eqs. ()5.23|) above, these equations give the full difference between the dimen- 
sionally regularized and the pure Hadamard accelerations, and they constitute the main new 



input of the present work. The bare dim. reg. results, c! 



mn 



+ I^Cmn, read therefore 



-21 



-21 



J'dr 

-21 



-03 



11 33 2^ 



[vl2 - 5(ni2t;i2)'] + 



24 



1463 
T2~' 



^1^2) 



24 



1463 
y - 33 In {qr'^2) 



8959, ,2 2317, , 4193, 

('^12t^l) + (^12^^l)(«12^^2) 77^(^12^^2j 



12 



, 1325, 
(^12^12) + -^{ni2Vi} 



1463 
43 



U^12f2 



18(t;if2) - 9t;2 - (riuVi) + 2(ni2t;i)(ni2t;2) + y(^i2^^2)' 



(6.47a) 
(6.47b) 
(6.47c) 
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c'ot = = <nuvi) + 5{n,2V2) . (6.47d) 

Note that in the final result the equality between c^^ and c^„, i.e., between B' and B" 
in Eq. ()5.3|) . is recovered. The bare dim. reg. kinetic coefficients ()6.47|) contain poles 
(X 1/e but do not depend anymore of the arbitrary Hadamard regularization scale r[ which 
appeared in (j6.45|) . As in the case discussed above of the static coefficients the previous 
kinetic coefficients do not involve any adimensionalizing length scales in the logarithms of 
ri2 they contain. This is consistent with the fact that it is the combinations and 

c^^ which have the same physical dimension as their d = 3 counterparts. 

Finally, given a specific choice of shift, say the MS one, Eq. ()6.37|) . the renormalized 
kinetic coefficients are obtained by adding to ()6.47|) the velocity-dependent part of the effect 
of the shift, i.e., the second term on the R.H.S. of Eq. (IFT^ [with, say, e^{d) = (11/3) P]. 
Our final (MS) results for the renormalized kinetic coefficients are found to be 



cf^ = -22 [vl, - 5(n,,.,,)^] In + - ^K-2) 

979 2 6275, ,3 1646, , 2851, ,2 , 
+ -^^2 - ^^K2fi) + ^-{ni2Vi){ni2V2) - -^{ni2V2) , (6.48a) 

4f = cT = -44(ni2t;i2) In + ^{n,2V,) - ^-^{n,2V2) , (6.48b) 
— 43 

cfi^ = 18{viV2) ~9vl- {nuvif + 2{n^2Vi){ni2V2) + Y(rii2t^2)' , (6.48c) 

Co? = c'^ = 4(ni2t;i) + 5{n,2V2) • (6.48d) 



When comparing these results with the ones of Ref. [23, Eqs. (ICTll above, one remarkably 
finds that our previously derived link (j6.4Uj) is necessary and sufficient for ensuring the full 
compatibility between the renormalized dim. reg. results and the corresponding Hadamard 
reg. ones. Note that the rational coefficients entering the dim. reg. results are often simpler 
than the coefficients entering the equations of motion of 2^ . 

The results Eq. ()6.48|) complete our check of the full consistency of the dim. reg. eval- 
uation of the 3PN equations of motion, and the proof of Theorems 1 and 2 stated in the 
Introduction. 



VII. CONCLUSIONS 

We have used dimensional regularization [i.e. analytic continuation in the spatial dimen- 
sion d) to determine the spacetime metric and the equations of motion (EOM) in harmonic 
coordinates, of two, gravitationally interacting, point masses, at the third post-Newtonian 
(3PN) order of General Relativity. Our starting point consisted in writing the 3PN-accurate 
metric g^v{x) in terms of a certain number of "elementary potentials" V, Vi, Wy, ■ ■ ■ , satis- 
fying a hierarchy of inhomogeneous d'Alembert equations of the form □(potential) = source. 
The sources of the latter equations contain both "compact" terms, i.e., in the present case 
contact terms of the form J^[V{'x), V^(x), Wij{'x), ■ ■ ■]S^'^\x — yi), and nonlinearly generated 
"non compact" terms of the typical form, say, 5(potential)c?(potential). This representation 
of the 3PN metric, as well as the associated iterative way of solving for the potentials [using 
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the time-symmetric Green's function D"^ = + c~^A~^9j^ + (9(c~^)] is a direct general- 
ization of the one used in Ref. jl^. However, it has been crucial for our work to determine 
(in Section |n} the dependence upon the dimension d of the coefficients appearing in this 
representation, as well as the d-dependence of the kernels expressing the operators and 
A^^ in x-space. 

By studying the structure of the iterative solution for the metric, and that of the corre- 
sponding EOM (which are conveniently pictured by means of diagrams, see Figs. ^IHI); we 
determined, in the form of a Laurent expansion in e = c? — 3, the pole part of the metric 
Qfiuix), and the pole and finite parts of the EOM, namely = Aaiyt, ^b) where a,b = 1,2 
and Vb = dyb/dt. [See, however. Appendix O where the basic quadratically non-linear kernel 
^(x, yi,y2) is computed in any d dimensions, not necessarily close to 3.] Our calculations 
relied heavily on previous work in = 3 [22] , and were technically implemented in two steps 
(at least for the determination of the EOM, which are more delicate; the determination of 
the pole part of the metric uses only the second step). 

(i) A first step consisted of subtracting from the final, published results for the EOM 
ji^ . seven contributions that were specific consequences of the use of an extension of the 
Hadamard regularization method [ij, 0| (which included an extension of the Schwartz 
notion of distributional derivative). The result of this first step is referred to as the "pure 
Hadamard-Schwartz" evaluation of the EOM. 

(ii) The second step was the evaluation of the difference between the dimensional regu- 
larization of each contribution to the EOM (written in terms of the iterative solutions for 
the various potentials V, Vi, Wij, ■ ■ ■), and the corresponding "pure Hadamard-Schwartz" 
contribution obtained in the first step. This difference is obtained, similarly to the method 
used in [35^, by splitting the ci-dimensional integral into several pieces, and by carefully 
analyzing the terms due to the neighborhoods of the two singular points yi, y2 (including 
possible c?- dimensional distributional contributions). 

Concerning the "bulk" metric g^^{x), at a field point away from the singular particle 
world-lines ya{t), we derived only the pole part, that is the coefficient of 1/e in the Laurent 
expansion of gp,v{x] e). We found that at the 3PN order only the time-time component of the 
metric gm{.x) contained a pole [see Eq. ()6.9|) ]. For the EOM we derived both the pole part and 
the finite part, i.e. ~ + 0{e). The parts of the EOM for which the regularization 

was delicate are given by the nine coefficients C31, ■ ■ ■ , C03 defined in Eq. ()5.6|) . Our complete 
results for the dimensionally-regularized values of these nine "delicate" coefficients are given 
in Eqs. (IH:^ and (lOTfl . 

We proved that the pole parts of both the metric and the EOM can be "renormalized 
away" by suitable shifts of the world-lines of the form y\^^'^ = y^a'^ + ^ai^) , where yjj^'''' is the 
original world-line on which are initially concentrated the (5-function sources representing 
the point masses, where the shifts ^ai^) ~ + e° -|- 0{e) are of the 3PN order, and where 
the EOM of the renormalized world-line y^^^ is finite as e — >■ 0. The general form of the 
needed shifts is given by Eq. ()6.15|) with ()6.10p . The renormalized EOM corresponding to 
the "modified Minimal Subtraction" scheme ()6.37p are given by Eqs. ()6.39p and ()6.48|) . 

The finite renormalized 3PN-accurate EOM obtained by using the general (two- 
parameter) renormalization shift (j6.41|) - (j6.43p were shown to be equivalent to the final 
(three-parameter) EOM of 2^1 if and only if the Hadamard-undetermined dimensionless 
parameter A which entered the latter equations takes the unique value A = — This 
value is in agreement with the result of a previous dimensional-regularization determination 
of the Arnowitt-Deser-Misner Hamiltonian (in ADM-like coordinates) j3^, which led to the 
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unique determination of the ADM analogue of A, namely Ug = 0. The value for A or Ug is 
also in agreement with the recent work which derived the 3PN equations of motion in 

harmonic gauge using a surface-integral approach. Our result provides an important check 
of the consistency of dimensional regularization because our calculations are very different 
from the ones of jHI, notably we use a different coordinate system and a different method 
for iterating Einstein's field equations. However, the applicability of our general approach 
to higher post-Newtonian orders remains unexplored. 

Finally, the present work opens the way to a dimensional-regularization determination 
of the several unknown dimensionless parameters that were shown to enter the Hadamard- 
regularization of the 3PN binary's energy flux (in harmonic coordinates) ji^, 113 • The com- 
pletion of the 3PN energy flux is urgent in view of its importance in determining the gravita- 
tional waveforms emitted by inspiralling black hole binaries, which are primary targets for the 
international network of interferometric gravitational wave detectors LIGO/VIRGO/GEO. 
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APPENDIX A: THE d-DIMENSIONAL METRIC AND GEODESIC EQUATION 



We give in this Appendix several expanded expressions which are too lengthy to be 
included in the body of the article. The expanded form of the metric (j2.11|) is easier to 
compare with the literature, and notably with Eqs. (3.24) of Ref. [j^ : 
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The inverse metric is such that g^vd"'^ = 5^ in + 1 space-time dimensions. In terms of the 
modified Newtonian potential V defined in Eq. fl2.10|) above, it reads: 



,00 



10 



(A2a) 



= g (d-2)c^ 



1 + 



1 /rf- 1 1/ 



2 \d-2 c2 



+ 



8-R,; 



16 

77 



X 4 ^ 16 



(A2b) 
(A2c) 



2(d-2) 

Note the change of signs in the exponentials [with respect to the covariant metric ()2.11|) ] 
in front of W^Vj in Eq. (IX2E|) . as well as for the Oillc") and C(l/c^) terms in Eq. (IX2c|) 
Note also that the ^Vj- contribution to gij has disappeared in the inverse spatial metric g^^ 
The full post-Newtonian expansion of this inverse metric reads: 
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The determinant g = detg^^ of the metric is a useful quantity, notably to compute the 
"gothic" metric g^*^ = ^^—g g^'^ ■, which is the natural variable when using the harmonic- 
coordinate system. The simplest way to compute it is to use the exponential form ()2.11|) 
of the metric, and to perform a cofactor expansion across both the first line and the first 
column: 
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Since gij = exp 
gi^k j^i reads 
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Therefore, the determinant of the full metric is given by 
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where we have used the fact that goi = 0{l/c^). Note that this formula suffices to compute 
g up to order 0{l/c^) included if one knows the spatial metric gij up to this same order. At 
the 3PN order, we get easily 
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where we used the expansion det(]l + M) = 1 + TrM + (9(M^), valid for any matrix M 
whose entries are small with respect to 1. We can now compute the square root of this 
determinant, and give its full post-Newtonian expansion: 
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The gothic metric = y/—g g^'^ can now be written easily by combining Eqs. ()A2|) or ()A3|) 
with ()A8a|l or ()A8b|l . We shall not display here the explicit results, since they were not 
directly useful for the present article. Let us however quote the expression of the Ricci 
tensor in terms of the gothic metric, in c? + 1 space-time dimensions and in any gauge: 
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As usual, a comma denotes partial derivation, and g^u = gpLv/yf—g is the inverse of . In 
terms of the gothic metric, the harmonic gauge condition ()2.H) takes a particularly simple 
form: 

fl^% = -v^^/°^r^/5 = o. (Aio) 

This is the reason why this gothic metric can be useful to write the field equations. Note 
that several terms of Eq. ()A9j) vanish in this gauge, namely the first two (involving second 
derivatives) and those proportional to q"'^ ^ in the second line. Nevertheless, this expression 
for Rfj_y is slightly more complicated than the one we used in Section |Tll above, Eq. ()2.2j] . 
which does not depend explicitly on the spatial dimension d. It should be noted that many 
equations given in the book ji^ are erroneous when d ^ ?> (z.e., when + 1 = 7^ 4, in this 
book's notation), including Eq. (I, 14, 30) in ^4^ which gives the Ricci tensor in terms of 
the gothic metric. 

Let us end this Appendix by displaying the full expansion of the geodesic equation ()2.19p . 
or more precisely of the vectors P* and F*, quickly illustrated in Eqs. (I2.22j) . The following 
expressions are d- dimensional generalizations of Eqs. (3.35) of Ref. [22], and we keep the 
same writing and order of the terms to ease the comparison. The "linear momentum" 
reads 
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This c?- dimensional expression actually allows us to understand better some of the numerical 
coefficients found for c? = 3 in Ref. j22|. For instance, we find that a factor 33 comes from 
the expression 3{5d — A), and that a factor 49 comes from {3d — 2)"^. The full post-Newtonian 
expansion of the "force" P* is given by an even longer formula: 
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APPENDIX B: USEFUL FORMULAE IN d DIMENSIONS 

This appendix is intended to provide a compendium of (mostly well-known) formulae for 
working in a space with d dimensions. As usual, though we shall motivate some formulae 
below by writing some intermediate expressions which make complete sense only when d 
is a strictly positive integer, our final formulae are to be interpreted, by complex analytic 
continuation, for a general complex dimension, G C. Actually one of the main sources 
of the power of dimensional regularization is its ability to prove many results by invoking 
complex analytic continuation in d. 

We discuss first the volume of the sphere having d — 1 dimensions {i.e., embedded into 
Euclidean ci-dimensional space). We separate out the infinitesimal volume element in d 
dimensions into radial and angular parts, 

d'^x = r'^-Mrdnrf_i, (Bl) 

where r = |x| denotes the radial variable {i.e., the Euclidean norm of x G M.'^) and dQ^-i 
is the infinitesimal solid angle sustained by the unit sphere with d — 1 dimensional surface. 
To compute the volume of the sphere, fld-i = J dfld-i, one notices that the following ci- 
dimensional integral can be computed both in Cartesian coordinates, where it reduces simply 
to a Gaussian integral, and also, using ()B1|) . in spherical coordinates: 
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where T in the last equation denotes the Eulerian function. This leads to the well known 
result ^ 

For instance one recovers the standard results Q2 = 4:1t and Qi = 27r, but also Qo = 2, which 
can be interpreted by remarking that the sphere with dimension is actually made of two 
points. If we parametrize the sphere fld-i in — 1 dimensions by means of d — 1 spherical 
coordinates 6d-i, ' ' , which are such that the sphere in d — 2 dimensions is then 
parametrized by 6 d-2, ^d-s, ■ ■ ■ , and so on for the lower-dimensional spheres, then we find 
that the differential volume elements on each of the successive spheres obey the recursive 
relation 

dnd-i = (sin Od-i)''-'' d9d-idnd-2 . (B4) 

Note that this implies 

^= rded-i{smed-,f-'= r dx{l-x^Y-^ , (B5) 

Jo J -I 

which can also be checked directly by using the explicit expression flB3|l . 

Next we consider the Dirac delta-function ^'•'^•'(x) in d dimensions, which is formally 
defined, as in ordinary distribution theory [s^, by the following linear form acting on the 
set V of smooth functions G C°°(R°') with compact support: 'iip G V, 

< 5^'^\ ^>= J d'^x5^'^\x)ip{x) = ^{0) , (B6) 

where the brackets refer to the action of a distribution on ip E V. Let us now check that 
the function defined by 

u = fcr^-^ (B7a) 

k = , (B7b) 

TT 2 

[where r is the radial coordinate in d dimensions, such that = 'Y2i=ii^''Y] is the "Green's 
function" of the Poisson operator, namely that it obeys the distributional equation 

Au = -An6'-'^\^) . (B8) 

For any a G C we have Ar" = a{a + d — 2) r"~^, thus we see that Au = in the sense of 
functions. Let us formally compute its value in the sense of distributions in x-space. [The 
usual verification of ()B8|) is done in Fourier space.] We apply the distribution Au on some 
test function G "D, use the definition of the distributional derivative to shift the Laplace 
operator from u to (p, compute the value of the d-dimensional integral by removing a ball of 
small radius s surrounding the origin [say B{s)], apply the fact that Au = in the exterior 
of B{s), use the Gauss theorem to transform the result into a surface integral, and finally 
compute that integral by inserting the Taylor expansion of ip around the origin. The proof 
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of Eq. ()B8|) is thus summarized in the following steps: 
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= nd^Jz{2-d)^{o) 

= -AnifiO). (B9) 
In the last step we used the relation between k and the volume of the sphere, which is 

k^d-i = ^. (BIO) 

From u = kr'^~'^ one can next find the solution v satisfying the equation Av = m (in a 
distributional sense), namely 

v= „ . (Bll) 

2(4 -d) ^ ^ 



From (jBlljl we can then define a whole "hierarchy" of higher-order functions w, ■ ■ ■ satisfying 
the Poisson equations Aw = v, ■ ■ ■ in a distributional sense. 

However, the latter hierarchy of functions u,v,---is better displayed using some different, 
more systematic notation. This leads to the famous Riesz kernels, here denoted in d- 
dimensional Euclidean space [i^ . [These Euclidean kernels differ from the Minkowski kernels 
also introduced by Riesz, and alluded to in the Introduction.] These kernels depend 
on a complex parameter a G C They are defined by 

S^f{x) = K^r''-^ (B12a) 

= (B12b) 

2-vr.r (f ) 

For any a G C, and also for any c? G C, the Riesz kernels satisfy the recursive relations 

= -Si'^ ■ (B13) 

Furthermore, they obey also an interesting convolution relation, which reads simply, with 
the chosen normalization of the coefficients Ka, as 

^i'^ * 5f = &s ■ (B14) 



When a = we recover the Dirac distribution in d dimensions, 6q = Kqt = 6^"^^ (the 
coefficient vanishes in this case, = 0), and we have u = Air v = —Arc ■ ■ ■ . 
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The convolution relation ()B14|) is nothing but an elegant formulation of the Riesz formula 
in d dimensions. To check it let us consider the Fourier transform of in d dimensions, 



/«(k) = y rf^x|x|"e-*''-^. (B15) 

Using ()B1|) we can rewrite it as 

Uk) = / drr°+^-i / d^]rf_le-^'^•^ (B16) 







in which the angular integration can be performed as an application of Eq. ()B4|) . This yields 
an expression depending on the usual Bessel function, 

j ^^^]d_le-^'^•^ = fid_2 d9,,, (sin e'^'^^^^"^^^- = (27r)MA:r)^-^ J._,(A;r) , 

(B17) 

where k = |k|. The radial integration in Eq. ()B16|) is then readily done from using the 
previous expression, and we obtain 

/„(k) = 2°+'^7^i^,i^r"-^ (B18) 

where the factor in front of the power k~°'~'^, say A^, is checked from the Parseval theorem 
for the inverse Fourier transform, which necessitates that AaA^a-d = (2vr)'^. Finally we can 
check the Riesz formula by going to the Fourier domain, using the previous relations. The 
result, 

d'^ry, = vr^ V X ^jJ ri^^^'^ (B19) 

' ' F (-f ) F (-f ) r (^^±|t2rf) 12 , K J 

is equivalent to Eq. ()B14|1 . 

A set of formulae concerning symmetric-trace-free (STF) multipole expansions in d di- 
mensions is presented next, without proofs. We use the multi-index notation L = ii ■ ■ ■ if, 
more generally the notation is the same as in Appendix A of |47|. In particular denotes 
the STF projection of ul = rii^ ■ ■ -rii^, [|] means the integer part of |, T{jj...j^} denotes the 
(unnormalized, minimal) sum of Ti^^^^^.i^^^^ where the ex's are permutations of the indices 
such that T{jj...j^} is fully symmetric in L (for convenience we do not normalize the latter 



We adopt for the Bessel function the defining expression 

■''-M- r(.!f)r(i) 

To obtain Eq. IjBf 8|l we employ the integration formula 



2 ^—tZX 
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sum, for instance S^ijUk} — SijUk + SikUj + SjkTii). 

[|] 

til 

k=0 

where the coefficients are 

. _ ir{l + £-2k) 



' 2^ r(f + £-i) ■ 

fn particular (the brackets ( ) surrounding the indices mean the STF projection) 



niUiL — nr . 

d + 2e-2 

Spherical averages: 
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STF decomposition of a scalar function: 
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Decomposition of a function F(n.n') in terms of Gegenbauer polynomials: 
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Integration formulae: 
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APPENDIX C: EXPLICIT FORM OF g = A ^ (^r^ IN d DIMENSIONS 

A very important technical fact which allowed one to compute analytically the d = 
3 equations of motion is the possibility to obtain explicitly the quadratically non linear 
potentials, i.e., to evaluate in closed form the integrals appearing in the PN expansion of 
the cubic- vertex diagram of Fig. [3 



X 



FIG. 7: Cubic- vertex diagram related to function g = A ^(r^ '^r2 



At the lowest approximation in the c ^ expansion, the diagram of Fig. [7| leads, in = 3, to 
the integral 



\ri ra/ 
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By definition, the Gegenbauer polynomial Cj{x) is the coefScient of in the expansion 



(l - 2xa 



2\ — 1 



(d) — — 1 

The particular polynomial {x) = C/ (a;) represents an appropriate generalization of the Legendre 
polynomial in d dimensions [indeed {x) = P£{x)]. 
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which was (probably) first evaluated by Fock in 1939 ( "Sur le mouvement des masses finies 
d'apres la theorie de la gravitation einsteinienne" |4S]), with the simple result 



y{d=3). 



x) = ln(ri + ra + ^2) • (C2) 



Remembering that ri = |x— yi|, r2 = |x— ya] and ri2 = |yi— y2| the combination ri+r2+ri2 
entering the logarithm in Eq. ()C2|) is simply seen to be the perimeter of the triangle joining 
the three spatial points x, yi and y2 entering the (Newtonian approximation of the) diagram 
of Fig. [7| At the (9(c~^) level of the PN expansion of Fig. [71 there enter several new integrals 
which can be reduced to 

fid=3) = 2A-'g^''='\ (C3) 
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/g=3) = A-M - ) and 1 ^ 2 . (C4) 



together with 



The explicit evaluation of the integrals ( )C3D . ()C4|) is also possible, as was shown in 
Refs. ji^,!!^] (drawing on earlier works [50|,|51|). In this Appendix we shall explicitly eval- 
uate the c?-dimensional generalization of (jClll . It will be clear, however, that our method 
can be rather straightforwardly generalized to the (9(c~^) diagrams contained in Fig. [7| i.e., 
to the c?- dimensional generalizations of ()C3|) - ()C4|) . 

For our present purpose it will be more convenient not to include the two factors of k 
that accompany the two propagators issued from 1 and 2 in Fig. |3 We shall therefore define 

g{^) ^ A-\rl-'^ rl~') . (C5) 

The method we present here consists of four basic steps: (i) expand the integrand in series 
and construct a corresponding series for a particular solution g^^^^ = (r^"'^ r2~'^)part, 
(ii) resum the series to get an explicit line- integral form of gpart) (iii) compute Ag^^^^ in 
a distributional sense to discover that it satisfies Ag^^^^ = r'l^'^ r\~'^ + S where 5 is a 
distributional source (localized along a line), and finally (iv) subtract A"^^ (which is given 
by another line-integral) from gp^^j.^ to get g as a sum of line-integrals (which are expressible 
in terms of one special function of one argument). What is crucial in the argument is 
the uniqueness of the global solution (decaying at infinity) of any (distributional) Poisson 
equation Aip = a when the (distributional) source decays fast enough (or, at least, does not 
grow too fast) at infinity. In our case, the sources a involved will have fast-enough decay at 
infinity if we analytically continue d toward large enough real parts (say 3? [d] > 3). 

There are several ways of implementing our method. For instance, we could start by 
expanding rj"*^ in the source of ()C5|) in powers of ri, such an expansion being valid only 
in a neighborhood of yi. Namely, we have the (i-dimensional generalization of the familiar 
d = 3 Legendre-polynomial expansion of |x — ya]"^ near x = yi (more precisely in the ball 
ri < ris) 



Here, we denoted for visual clarity pf\x) = Cf where CJix) is a Gegenbauer polyno- 

1/2 fS) 

mial such that (x) = PI \x) is the usual d = 3 Legendre polynomial [see also Appendix 
IbI above] . The quantity C\ in ()C6|) denotes the cosine of the angle d\ between x — y 1 and 
y2 — yi- The notation we shall use is summarized in Fig. |H1 



d 
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yi ri2 y2 

FIG. 8: Notation for various angles and distances, yi and y2 denoting the positions of the two 
particles, and x the field point. 

When inserting the local expansion jnni) into the source of ()(^5|1 we are led to solving (locally) 
an equation of the form A g^^^ = ^ aeri+'^-'^P^'^\ci). However, using the general formula 

A(r^ nL) = iX-£)iX + d-2 + i) r^-^ til , (C7) 

we know a particular solution oi Aip = r^fii, namely 

The formulae apply to any source with fixed multipolarity {£) and a power law 

dependence on a radius. In particular, they apply when r — > ri, A —>■ i + 2 — d and 
fiL = nf^'"^''^ — > P^'^\ci) (because a generahzed Legendre polynomial is just proportional 
to the contraction of an STF-projected multi unit vector onto a fixed "z" direction; see 
Appendix El- This leads to a corresponding expansion of a local solution gi^ci (near yi) of 
Agioci = rl'^^rj^'^ of the form 



S-.-|f^Ej^^Pf'fe). (C9) 



In order to proceed further, we now need to resum the expansion This is done by a trick 

introduced, in a similar context of resummation of multipolar expansions containing extra 
^-dependent denominators, by Ref. js^. One introduces some radial-integration operators 
Ra [0] (r) = dX X" 0(A r). For instance, in the context of ()C9|) . one replaces r{/{i + 1) by 

i?o[?"i] = Jq dX {XriY or equivalently rf^^/(£+ 1) by J^^ diii{. This transforms back the 
multipolar series appearing in ()C9|) into the original "Legendre" series entering Eq. ()C6|) . 
This allows one to write giod as a simple line-integral: 



2(4'-,)y/"l^"-^'l"- (CIO) 

Here, y^^ is a point on the segment joining yi to x, located a distance ii away from yi. It is 
more convenient to replace the line- integration over the dimensionful length ii (0 < £i < ri) 
into an integration over the dimensionless parameter a = ii/ri (0 < a < 1). This leads to 
the explicit expression 

y^i = ax+(l-a)yi. (Cll) 

The resummed line- integral expression ()C10|) allows one to define everywhere gi^d, including 
in the domain ri > ru where the original series ()C9|1 was not convergent. Having in hands 
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such a global definition of gi^d then allows one to compute its Laplacian, in the sense of 
distributions, and to see how it differs from r\~'^r\~'^. The calculation of Agj^^,^ is done by 
techniques similar to the ones used in Ref. One needs to rewrite some terms in the form 
of a-derivatives. For instance, several of the terms appearing in Agj^^^ can be rewritten as 
the line-integral 



da 



2-d 







oa 



da 

da 



rr'rt\ (C12) 



where the last line-integral gave only the end contribution a = 1 corresponding to = 
X. Besides the terms yielding ()C12|) . i.e., the looked-for "source" of the complete g, the 
calculation of Ag^^^^ yields also the distributional source (where k = k/in entered through 

^^2-d ^ _Si<i)/k) 

d-d ^1 

daa^S^^'^iy^-y^). (C13) 
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2(4 -d)k Jo 

This is conveniently transformed by introducing f3 = 1/a (with 1 < f3 < +00) and 

y^^(l-/3)yi+/?y2, (C14) 

which varies along a semi- infinite line going from y2 to infinity along the direction y2 — yi, 
i.e., away from yi. This transformation allows one to rewrite ()C13|) in the more transparent 
form 



r 



A-d 
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At this stage, we recognize in ()(n5jl a very simple source, namely a uniform distribution 
of "mass" along the half-line along which (3 runs. This allows one to easily compute the 
unique, global (decaying at infinity) solution of the Poisson equation with source ()C15|) and 
to subtract it from g\oci to get the unique, global g in the form of two line-integrals: 

3— d /"OO 
\2-d '12 / j/7 l„ \2-d 



where dCa = \dya\ = ri da and dijs = \dy/s\ = ri2 djS. In other words, ()C16|) expresses g as, 
essentially, the difference between the Newtonian potentials generated by two uniform line 
distributions: a segment joining yi to x and the half-line starting from y2 in the direction 
away from yi. It is easily seen (modulo the slight delicacy of the logarithmic divergence of 
the potential of a semi-infinite line when d —>■ 3, i.e., the occurrence of a pole; see below) 
that the result ()C16|) yields, when d —>■ 3, the well-known result ()C2|) . [Actually, this was 
the way one of us (T.D.) had derived long ago for himself ()C2|) . unaware of its derivations 
in the literature.] 

The expression (IC16|) has the advantage of being explicitly regular (except at the point 
X = yi) in the ball ri < ri2. However, it has the default of treating dissymmetrically the 
two points yi and y2 (in spite of the fact that the result ()C16|) for g is, actually, symmetric 
under 1 <-> 2). One can derive an exchange-symmetric expression for g by modifying the 
first step of our method. Instead of expanding the source r\~'^ r1~'^ in the neighborhood of 
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X = yi, i.e., in a series of positive powers of ri, we can expand it in the neighborhood of 
infinity, i.e., in a series of negative powers of ri. Such an expansion is directly related to 
the expansions used in j^, which led to the decomposition of g('^='^) in two pieces denoted 
k and h there, where the source of h was a uniform mass distribution along the segment 
joining yi and y2. Let us briefly indicate the successive steps of this new calculation of g. 
Instead of the "local" expansion ()C6|) (valid for ri < ri2), one expands rg"*^ near infinity 
(ri > ri2) as 

£>0 



r2 = r 



?-'E(77) ^'f'fe). (C17) 



Solving term by term Ag^ = r\~'^ r1~''' "near infinity" by means of ((DHl), and transforming 
away the appearing ^-dependent denominators by means of Mi = r\~'^~^ / 3), 

leads to the following resummation of g^: 

A-d pco 

Here, Yq, is still defined by ()C11|) . but the parameter a now varies in 1 < a < +cx3 so that 
()C18|) is the potential of a semi-infinite line. Computing the distributional Laplacian of 
the particular solution g^, Eq. ()C18|) . leads to the presence, besides the looked- for source 



2 d „2 d Qf g^j^ additional distributional source localized now along the segment joining yi 



' 1 '2 

to y2, namely 

rf/35(")(x-y^), (C19) 





ri2 



2(4 -d)k 

where P = 1/a varies between and 1 and where is again defined by ()C14|) . It is then 
easy to subtract from g^ (which tends, in = 3, to the function ^(x, yi,y2) of jHH) the 
Poisson integral of the source ()(yl9|l (which is a uniform distribution along the segment 
yi-y2 and which tends, in c? = 3, to the function |/i(x, yi,y2) of [IH) to get the following 
alternative expression for g, 

„3—d />oo »,3— ci /•ri2 

|2-d , '12 / J/) l„ |2- 



« = -2(4^ /, - y-r' + ^ i |x - y,r' . (C2o) 

or, equivalently. 

The form ()C20|) - ()C21|) still does not look quite symmetric between 1 and 2 but a moment 
of reflection will show that it is. 

The two methods above have expressed g in terms of the Newtonian potentials generated 
by half-lines or segments, i.e., integrals of the type J |x' — y^^"'^ where y^ varies along a 
straight line (but where x' might be x or y2). Clearly, any such potential can be reduced 
(through linear decompositions) to the Newtonian potential generated by a half-line. Let us 
then consider a generic half-line starting at the point fq and going to infinity in the direction 
n, and let us consider the Newtonian potential generated by this half-line at the origin of 
the coordinate system (not located on the half-line). Let us denote Vi = yq + In, = \y(\, 
To = ro no, c = no ■ n (cosine of the angle 6 between the radius vector from the origin, 
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i.e., the "field point", towards the beginning point of the half-hne and the direction of the 
half-line, away from its beginning). Then it is easy to find that 



dirj-'' = ip{c)rl-'', (C22) 

where the function ip{c) is given by the integral 

'Pic) ^ / , , • (C23) 

Jo (l + 2cA + A2) — 

The integral ()C23|) converges for d > 3, has a pole oc l/(c?— 3) as c? ^ 3, and can be expressed 
in terms of hypergeometric functions, e.g. F2^i [|, |; ^] • It is, however, simpler to keep 
the form KJ^ .^° 

Finally, using the half-line potentials ()C22|) as building blocks one can write our result 
]20|) in the final, 1 <-> 2 symmetric, form 



„3-d 
12 



g(x,yi,y2) = ^^^^ [2rl-V(0)-rrV(ci)-rrV(c2 

3— d 3—d 

The quantities entering ()C24|) are those defined in Fig. |H1 notably Ci = cos^i, C2 = cos ^2, 
C12 = cos ^12, with r_|_ being the orthogonal distance between the field point x and the 
segment joining yi to y2 [with associated argument c± = cos | = in (p{c)]. Note the 
following properties of the function ip{c), 

ip{c) + ^(-c) = -^^^ , (C25a) 
(1 — c^j 2 

irfilrf^^) 1 

m = ^ ' ^(^) = ^- ^^25b) 

The simplest way to prove ()C25a|) is to notice that the Newtonian potential of an infinite line 
can be written either as twice that of two half-lines beginning at the orthogonal projection 
of the field point on the original line, so that 

7^ = 27=?' (C26) 



30 



The multipolar expansion of the function (f{c) reads 



On this expression one sees clearly the occurrence of the simple pole of Lp{c) when e = d — 3 ^ 0, which 
is given by the "monopolar" term = as 

^(c) = :^o^+O(eO) = l + 0(eO). 
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or as that of two other half-hnes obtained by a more arbitrary cut (under an angle 9 ^ ti/2 
and c = cos 6) . 

We can verify the c? — 3 hmit of Eqs. ()C16|) and ()C24|) by using the following e — > 
expansion of the elementary function '^{c), namely 

^(c) = --\n +0{e). (C27) 



e V 2 , 

To obtain ()C27|) we notice that the finite part of (p{c) when e 0, which is 

1 



V?o(c) = lim 

£^0 



(C28) 



can be re-expressed in the form of the following sum of two convergent integrals, 

;^ dX d\ f 1 



VI + 2cA + A2 io ^ VVl + 2cA + A2 
= -ln(^). (C29) 

Combining the expansion ()C27|) with the basic relations associated with the triangle of 
Fig. IHl our (i-dimensional expressions ()C16|) or ()C24|) are found to admit the expansion 

' Vlnf!l±^|±^V^(-), (C30) 



2e 2 \ 2 

which indeed reduces to the three-dimensional result ()C2|) modulo an additive constant 
linked to the 1/e pole. 

Nice as it is to have in hand an analytic expression for the d-dimensional basic non linear 
potential g, its practical utility in explicit computations of the (i-dimensional equations of 
motion is not evident because, contrary to the 3-dimensional expression ()C2|) . the expression 
is not explicitly regular along the yi — y2 segment. [The regularity of Eq. ()C24|) as r_|_ — >• 
comes by compensations between the three terms in the bracket, using ()C25a|) .] It would 
need some transforming [using ()C25a|) . and/or using the other expressions derived from the 
previous form ()C16|) . which are regular along the yi — y2 segment, but singular somewhere 
else] to write an explicit expression which is regular everywhere, except at the two isolated 
points yi and y2. 

Finally, let us just mention that the method explained above can, in principle, be straight- 
forwardly generalized to the computation of the higher post-Newtonian potentials contained 



Denoting S — ri + r2 + ri2 we have 

1 + ci _ ri + ri2 - r2 1 + C2 _ r2 + ri2 - ri 1 + ci2 _ ri + r2 - ri2 

S 2riri2 ' S 2r2ri2 ' S 2rir2 

The perpendicular distance r±_ is given by 

ri_ = ^V2(l + ci)(l + c2)(l + ci2) . 
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in the diagram of Fig. [7| For instance, in computing the d-dimensional analog of ()C3|) . say 
/ = 2A-ig = 2 A-2(ri"'^r2"'^), it is easy [by iterating JCHI)] to get the analog of JCi). 
Then a more complicated radial- integration operator (see, e.g., [52]) will allow one to resum 
the series to get a line-integral expression for /lod or foQ. We anticipate that a somewhat 
more delicate application of either A (to go back to g) or A^ (to go back to 2 rf^'^ r'^"'^) will 
yield additional line-distributed sources. It should then be a simple matter to compute the 
Poisson, or iterated Poisson, integral of these line-distributed sources. We leave an explicit 
study of these details to future work. 
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